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1.0 INTRODUCTION 

This report summarizes the results of several study tasks pertaining 
to various aspects of the Space Shuttle communication system. These tasks 
performed can be logically divided into the following six categories: 

1. Phase multiplexing for two- and three-channel data transmission. 

2. Effects of phase noise on the performance of coherent commun- 
ication links. 

3. Analysis of command system performance, 

4. Error correcting code tradeoffs. 

5. Signal detection and angular search procedure for the Shuttle 
Ku-band communication system* 

6. False lock performance of Costas loop receivers, 

A sununary of the phase multiplexing techniques for the Ku-band 
return link is presented in Section 2. 0 with detailed analyses in Appendices 
Af B, and C. The effects of phase noise on communication system perform- 
ance is summarized in Section 3. 0 with details in Appendix D. Analysis of 
the Range Safety command system is analyzed in Section 4. 0, Techniques 
for providing command authentication for Air Force missions are summar- 
ized in Section 4. 0 with details in Appendix E. Error correcting coding 
tradeoffs for the 50 Mbps Ku-band return link are discussed in Section 5. 0. 
Detailed analysis for predicting performance of the Shuttle convolutional 
coded communication links is presented in Appendix F, with the tradeoffs 
that were considered in the final recommendation of the 50 Mbps coded 
system identified in Appendix G. The signal detection analysis and angular 
search procedure required to point the Orbiter Ku-band antenna at the TDRS 
are discussed in Section 6. 0 with details in Appendices H and I* Finally, 
the potential false carrier lock conditions for the S -band and Ku-band link 
are summarized in Section 7. 0 with the detailed analysis in Appendix T. 


2. 0 PHASE MULTIPLEXING FOR TWO- AND THREE-CHANNEL 

DAT A TRANSMISSION 

Phase multiplexing offers a technique of accommodating multiple 
data channels with independent data rates and clocks as required on the 
Ku-band return link. Alternatives to phase multiplexing are time -division 
multiplexing and frequency-division multiplexing. Time division multi- 
plexing at the 100 Mbps data rate encountered on the Ku-band return link is 
extremely complex because of the requirement of reclocking the data. i<’re- 
quency-division multiplexing requires large bandwidths to avoid intermod- 
ulation products. The available bandwidth through the TDRS for the Ku-band 
return link is 225 MHz, which is just large enough to pass the 100 Mbps 
data rate, let alone frequency multiplexing additional data rates into this 
bandwidth. 

Appendices A, B, and C investigate the phase multiplexing modulation 
and demodulation. Appendix A investigates the general phase multiplexing 
techniques including conventional phase modulation with subcarriers, inter- 
plex, and hybrid interplex-quadrature modulation. Phase multiplexing using 
quadrature carriers is shown to have potential implementation advantages 
over both a conventional PM subcarrier approach and the interplex approach 
which has been utilized for deep space communications. To illustrate this 
potential implementation advantage, the particular cases of two- and three - 
channel phase multiplexing using various techniques are compared in terms 
of implementation and performance. It is demonstrated in Appendix A that . 
either the interplex or the hybrid phase multiplexing technique for transmis- 
sion of three independent and asynchronous digital data channels offers a 
significant performance advantage over a conventional PM subcarrier design. 
The hybrid phase multiplexing technique differs from interplex primarily 
in terms of implementation. Whereas interplex involves data channel com- 
bining at baseband and utilizes a linear phase modulator at IF, the hybrid 
technique involves data channel combining at IF and utilizes quadrature 
carrier modulation. Since the performance characteristics of these two 
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approaches are essentially identical, it appears that the selection o£ a tech- 
nique for any particular application should be dictated only by equipment 
design considerations. 

Appendix B investigates candidate receivers for a two-channel phase 
multiplexed Ku-band communication signal. It is shown that the fourth- 
power loop makes very ineffective use of the total available tracking power 
when the power division is near 4:1 . An additional disadvantage is that the 
stable lock points occur at multiples of tt/ 4, thereby allowing the possibility 
that the two data channels become interchanged. This could be detected and 
corrected in the system design, but only at the expense of additional complexity. 

The Costas loop which tracks on the high data rate signal of the 
unbalanced QPSK waveform is shown to perform satisfactorily. Appropri- 
ate error rate computations show that the Costas loop considered performs 
within a few tenths of a dB of the ideal receiver. Also, it is shown that the 
Costas loop has stable lock points at multiples of ir which guarantee that the 
two data channels cannot be sent to the incorrect data demodulator. Although 
it is usually expected that the power in both channels will always be present 
even if the data rates are reduced to aero, should the power in the high data 
rate channel fail, then the sign of the tracking signal is reversed and the 
stable lock angles are at tr/Z + kir, where k is 1, 2, 3, . . . . If this occurs, 
then the low data rate link will be processed into the high data rate demodu- 
lator. If the power in the low data rate link goes to zero, then the stable 
lock angles remain unchanged, so that the high data rate link continues to 
be processed into the high data rate demodulator. 

Appendix C presents the analysis of the recommended modulator and 
demodulator for the Ku-band return link. The exact analytical expression 
for the hard-limited three -channel hybrid phase multiplexed signal is derived 
From this analytical representation of the modulation, the performance of 
the Costas loop demodulator is analyzed. Analysis is presented that allows 
the optimization of arm filters of the Costas loop in terms of the high data 
rate and filter type. This optimization approach is new, since the usual 
analysis of the Costas loop considers that the arm filters are sufficiently 
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wide that the data modulation is undistortcd. It is shown that the optimum 
choice of the arm filters is narrow enough to actually distort the signal but 
the loop signal -to -noise ratio is maximized by limiting the noise that is 
passed through the filters. 

The result of carrier tracking phase error in the receiver introduces 
crosstalk from one data channel into another. In particular, the low data 
rate channels introduce crosstalk into the high rate channel in proportion 
to their relative received powers and the carrier tracking loop phase error. 
While Appendix A presents the crosstalk power as a function of fixed carrier 
tracking loop phase error, Appendix C presents the average crosstalk as 
evaluated by considering the loop phase error probability density function. 
The average crosstalk is important because it indicates the actual system 
performance degradation that can be expected. 

An area for further study of the phase multiplexing technique is 
numerical integration of the integrals in Appendix C that allow optimization 
of Costas loop tracking and performance degradation due to channel cross- 
talk. In addition, a similar analysis to Appendix C needs to be performed 
on the subcarrier tracking loop. Finally, an analysis of the false lock char- 
acteristics of the Costas loop with three -channel modulation needs to be 
performed similar to the analysis for bi-phase data in Appendix J. 
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3. 0 EFFECTS OF PHASE NOISE ON PERFORMANCE OF COHERENT 

COMMUNICATION LINKS 

Tran stricter phase noise limits the minimum probability of error 
achieved by a coherent phase -shift -keyed (PSK) communication system. 
Appendix D systematically considers the effects of transmitter phase noise 
onihe performance of two-way coherent communication links and indicates 
design approaches that may be used to reduce or minimize such effects. 

As a specific example of phase noise» the frequency synthesizer necessary 
for the S-band payload interrogator in the' DSN mode is considered* For 
this transmitter phase noise spectrum, the resulting degradation in bit error 
rate performance for baseband and PSK subcarrier channels is determined. 
For the forward portion of the two-way link, the transmitted phase noise is 
due primarily to jitter introduced by the oscillators and phase detectors 
used in the frequency synthesizer to generate the carrier frequency. For 
the return link, the transmitted phase noise is due primarily to the turn- 
around of the phase noise which is not tracked out by the forward link 
receiver carrier tracking loop. 

The analysis of the two -link with transmitter phase noise presented 

in Appendix D shows that there is a trade-off between the effect of phase 

noise on the forward link and the phase noise turned around for the return 

-4 

link* As a final result, it is foimd that at bit error probability of 10 , 

the performance degradation of the S-band payload interrogator in the DSN 
mode communicating with the DSN standard transponder is less than 0* 1 dB 
and 0.2 dB for the forward link and the total two-way link, respectively. 
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4, 0 COMMAND SYSTEM CONSIDERATIONS 

For Air Force missions, modifications are proposed in Appendix E 
to the Shuttle commeind system to provide command authentication. Three 
possible modifications are proposed. I*i each case, an authenticator vari- 
able up to 50 bits long is used to modify i;he parity portion of the command 
codeword. Of the three modifications proposed, the simplest of modulo -2 
adding the authentication variable to the parity is considered too easy to 
jam to be of much utility. The other two modifications were considered 
to be hard to jam. The simplest of these two techniques is presented in 
Figure 1, where the information word and the variable are encoded by com- 
pletely different (127, 50) BCH encoders. The two parities produced are 
modulo -2 added to form The codeword x = (i_,£) is transmitted. The 

same process is performed onthe received word in the error detector; the 
resultant parity is compared with the received parity £' . The remainder 
of Appendix E analyzes the susceptibility of jamming of this recommended 
modification and is found relatively Impossible to jam. 

An additional analysis that was performed on command systems was 
the performance analysis of the Range Safety Command System. A typical 
message format for the command system is shown in Figure 2. The sub- 
carrier frequencies (tones) are from 7.35 kHz to 13. 65 kHz, Each char- 
acter is made up of two tones transmitted simultaneously. The secure 
address is composed of 9 characters and 2 characters comprising the func- 
tional selection. 

To analyze the possibility of jamming the command system, an 
assumption must be made about the demodulator. It is assumed that, when 
the secure address is locked Into the command system, thresholds are set 
such that, at any given time if any tones are received that are greater than 
the thresholds besides the expected two tones making up the proper charac- 
ter, then the command will be rejected. If these thresholds are not set, 
then a jannmer could cause a command to be incorrectly accepted by placing 
jamming tones at each of the seven possible subcarriers. Using these 



(a) Encoder 


£’ 



(b) Detector 


Figure 1 . Command Authentication Technique 
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thresholds, the command system will accept a command if it received only 
two tones per character and only the correct characters corresponding to 
the secure address. Similarly, the command is accepted only if there are 
only two tones per character for the function selection. It should be noted 
that, with the use of these thresholds, it is easy to deny access to the com- 
mand system by the jammer simply placing jamming tones at each of the 
seven possible subcarriers causing all commands to be rejected. This type 
of jamming is not considered a threat, however, because of the small utility 
to the jammer. 

With the demodulator assumption of the use of thresholds, it is 
straightforward to calculate the possibility of a jammer causing an incor- 
rect command to be accepted. Assume a technically competent jammer 
(i, e, , possession of the encoder, transmitter, receiver, and decoder). 
There are seven possible tones taken two at a time ( only two) corresponding 

to Zl possible characters. For a nine -character secure address, there are 
9 11 

21 or 7.9x10 possibilities. To try all possible secure addresses, it 

7 

would require 2,4x10 hours of continuous transmission. In 30 minutes, 

4 . 

1.6x10 possibilities could be tried or the probability of trying the correct 

-8 

secure address is 2x10 . 
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5. 0 ERROR CORRECTING CODE TRADEOFFS 

To design a convolutional coded communication system, techniques 
are required to predict probability of bit error (P^^) versus energy per bit 
to single-sided noise spectral density (E^/N^) for the modulation/demodu- 
lation scheme employed. In addition, it would be desirable to predict the" 
nature of the errors out of the decoder. That is, do the errors occur in 
bursts and, if so, what is the probability of an error burst of a given length. 

There are three basic analytical techniques available for predicting 
versus E^/N^. The first technique uses Gallager’s exponential bound 
for random codes. This boimd allows quick comparisons of various modu- 
lation/demodulation and encoding /decoding schemes in terms of the basic 
coding parameters of code rate (R) and encoding constraint length (K). 
Therefore, it can be used as an initial optimizatiou technique. 

The second technique is to analyze the distance between codewords 
using a generating function (transfer function) developed by Viterbi. This 
analysis technique has been used to obtain tight upper bounds on the error 
probability for the binary symmetric channel (BSC), the additive white Gaus- 
sian noise (AWGN) channel, and general memoryless channels. These bounds 
can be used to optimize the communication system parameters for a particu- 
lar convolutional code and Viterbi decoding. While the bounds on error 
probability using the transfer function approach are very precise, the bounds 
are difficult to calculate for long constraint codes. Also, for even moderate 
constraint lengths, the transfer function bounds are evaluated in matrix form 
and the individual codeword distances are not available to determine the 
effect of decoder memory truncation or to determine the statistics of the 
errors at the output of the decoder. 

The third technique uses an approximation of the transfer function 
developed by Huth and Weber, which is easy to calculate even for long con- 
straint length codes. The resulting bounds are as tight as Viterbi transfer 
function bounds and the individual codeword distances appear in the bound 


I 
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so that It is possible to determine decoder memory truncation effects and 
output error statistics. 

Appendix F compares the three approaches for predicting perform- 
ance and presents results for Viterbi decoding bit error probability, burst 
error probability and decoder memory length. The predicted performance 
with the Huth-Weber approximations is compared with actual system tests 
and found to be very accurate. 

Error correcting coding is required for 50 Mbps data link from the 
Shuttle Orbiter through the Tracking and Data Relay Satellite System (TDRSS) 
to the ground because of severe power limitations. Convolutional coding has 

been chosen because the decoding algorithms (sequenticd and Viterbi) provide 

" 6 

significant coding gains at the required bit error probability of 10 and can 
be implemented at 50 Mbps with moderate hardware. While a 50 Mbps 
sequential decoder has been built, the highest data rate achieved for a 
Viterbi decoder is 10 Mbps. Thus, five multiplexed 10 Mbps Viterbi 
decoders must be used to provide a 50Mbps data rate. Appendix G dis- 
cusses the tradeoffs which were considered in selecting the multiplexed 
Viterbi decoder approach for this application. 

It is felt that the selected 50 Mbps encoding /decoding approach offers 
a reasonable compromise between performance and system complexity. 

The coding gain provided by this approach is relatively high (5.3 dB at a 

-6 

probability of error of 10 ), the system degrades gracefully (a coding gain 

-2 

of 2.7 dB is available at a probability of error of 10 ), the spaceborne 

encoder is extremely simple to implement, and the ground-based decoders 
are available commercially emd can be combined in a straightforward 
manner. Such an approach is not limited to 50 Mbps data rates. It appears 
reasonable to extend the parallel 10 Mbps Viterbi decoder concept to sys- 
tems requiring operation at data rates well in excess of 100 Mbps, although 
reliability becomes of increasing concern as the ntimber of parallel decoders 
becomes large. 
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6. 0 SIGNAL. DETECTION AND ANGULAR SEARCH PROCEDURE 
FOR SHUTTLE KU-BAND COMMUNICATION SYSTEM 

To establish the Husband communication link between the Shuttle 
Orbiter and the TDRSi the antennas of these two satellites must be pointed 
at each other. Various techniques for pointing the TDRS antenna at the 
Orbiter have been considered, and the latest information available to 
Axiomatix indicates that initial pointing can be accomplished without 
requiring the Shuttle antenna axis to be in perfect alignment with the TDRS 
antenna axis. During this initial pointing, however, the effective EIRP 
available from the TDRS is only 40 dBw, instead of the nominal 50 dBw. 

But. before the TDRS makes its final pointing correction, the Shuttle 
antenna must search out the random designation errors which exist between 
its computer estimated line-of-sight (LOS) to the TDRS and its actual point- 
ing direction. Appendix H analyzes several procedures for performing this 
search and estimates the time requirements for each. 

Since the object of any search is to Increase the energy of the 
received signal, there must be a signal detector to stop the search when 
the signal energy has exceeded a preset value or threshold. A simple and 
fast circuit for determining this is a square -law device followed by an inte- 
grating (LPF) filter and a threshold crossing indicator. Ideally, the pre- 
detection IF filter is matched to the RF spread spectrum; however, a prac- 
tical Butterworth or Chebishev filter will degrade the predetection SNR 

only slightly. For B__ = 20 MHz, this loss is 0, 9 dB, so for P /N. = 

ix* rec 0 

64,8 dB-Hz (assumes EIRP = 40 dBw), the predetection SNR is about 

-7. 6 dB. Appendix H computes the required integration times for various 

values of probability of detection. P , and for a false alarm probability, 

-9 ^ 

p£^ = 10 , Assuming that the design is to detect the signal within a 1 dB 

beamwidth with = 0. 999, then the corresponding integration time is 
220 microseconds, which is short compared to the dwell times of mechan- 
ically scanned antennas. 
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Having established the criteria for detecting the presence of the 
signal within the 1 dB antenna beamwidth, Appendix H considers the angular 
search procedure which insures that the designation uncertainty volume is 
scanned to provide a predetermined probability of true LiOS intercept. The 
Shuttle antenna beam is spiraled out from its designated LOS to search for 
the true LOS. The spiral scan procedure is analyzed in Appendix I in terms 
of dwell timef total scan time, and scan rate. Several numerical examples 
pertinent to the Shuttle Ku-band radar /communication system are presented. 

When the energy detector commands the Shuttle antenna to stop 
scanning cind to start tracking, there may be a pull-in transient. A con- 
servative estimate is that the pull-in time should be approximately one 
third of the inverse of the servo loop beuidwidth. For a 1 Hz loop, this 
implies that about 0, 3 second will be required for the servo loop to go from 
scan to track. Consequently, the time that the true LOS stays within the 
moving antenna beam, called the dwell time, must be at least 0. 3 second. 

The analysis presented in Appendixes H and I indicate that the times 
required to search the uncertainty angles for TDRS antenna acquisition are 
short compared to radar search times. This is due to the fact that the SNR 
available from TDRS for signal detection during acquisition is much higher 
than the SNR of radar targets. Consequently, the scan time limitation for 
TDRS acquisition is determined by servo considerations rather than signal 
detection. The approach for radar search is the same as that for TDRS 
acquisition, and only the dwell time requirements are different. We there- 
fore recommend that the radar search pattern be adopted for TDRS 
acquisition. 
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7. 0 FALSE LOCK PERFORMANCE OF COSTAS LOOP RECEIVERS 

False lock can occur in Costas loop receivers when the frequency 
search requirement due to Doppler and oscillator instabilities is large 
enough to include false lock points. Appendix J demonstrates under what 
conditions the Costas loop can exhibit a false lock condition in the sense 
that the loop can lock up on a data sideband. In particular, for the case of 
random data, the loop can theoretically lock at a frequency which is any 
integer multiple of half the data rate. The required frequency search for 
both the S-band and Ku-band is larger than half of the data rate* Therefore, 
false lock conditions are expected on both sides. 

Appendix J derives an expression for the DC output of the lock detector 
in terms of the Fourier transform of the symbol pulse shape, the squared 
magnitude of the lowpass arm filter transfer function, and the loop phase 
error. This expression is then evaluated at the true lock and false lock 
conditions and numerical results are given for the case of Manchester coded 
data and a single pole (RC) arm filter. It is found that, by widening the arm 
filter bandwidths beyond the optimum value which yields minimum squaring 
loss for the Costas loop, the separation ' eiween the true lock and false lock 
voltages at the lock detector increases. Therefore, the true lock case can 
be more readily distinguished from the false lock condition. However, 
widening the arm filter bandwidth also allows more noise to pass through 
to the lock detector, thus degrading the signal -to -noise ratio at the input 
to the threshold detector. Therefore, widening the arm filter bandwidth 
is potentially a solution to the false lock detection problem up to the point 
where the signal -to -noise ratio which distinguishes between true lock and 
out-of-lock (as opposed to false lock) becomes too low for a specified level 
of lock detection performance. 
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APPENDIX A 

PHASE MULTIPLEXING FOR TWO- AND THREE-CHANNEL 

DATA TRANSMISSION 

Sergei Udalov 

This appendix considers various techniques for simultaneous trans- 
mission of two or three independent and asynchronous digital data channels. 

A phase -multiplexing technique, using quadrature carriers, is shown to 
have potential implementation advantages over both a conventional PM sub- 
carrier approach and the interplex approach which has been utilized for 
deep space communications. The phase -multiplexing technique is a hybrid 
approach which has some of the features of both quadrature modulation and 
interplex methods. Performance and implementation considerations are 
discussed for each of the various approaches. 

1.0 INTRODUCTION 

Phase -multiplexing offers a technique for accommodating multiple 
channels with independent data rates and clocks. This technique is especially 
important when high data rates are involved. Time-division multiplexing at 
high data rates is complex because of the requirement for reclocking the 
data. Frequency-division multiplexing requires large bandwidths to avoid 
intermodulation products. This appendix investigates the general phase 
multiplexing techniques, including conventiotial phase modulation with sub- 
carriers, interplex [1-4], and hybrid interplex-quadrature modulation. 

The particular cases of two- and three -channel phase multiplexing using 
the various techniques are compared in terms of implementation and 
performance. 

2.0 PHASE MULTIPLEXING PRINCIPLES 

2. 1 Conventional Phase Multiplexing 

The general form of a phase -modulated RF signal of carrier fre- 
quency o> is 
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z(t) = ^2P sin [co^t + 0 (t)] . (1) 

In the conventional modulation system, adapting the nomenclature used in 
[l], the combined phase modulation signal resulting from N channels is 
defined by the term: 

N 

0(t) = y\ (t) , (2) 

n=l 

where 6 is the modulation angle associated with the nth channel and 
n — 

S„(t) := d„(t) sq (e>t) == + 1 , ( 3 ) 

ntin 

which defines the nth data stream biphase modulated on a squarewave sub- 
carrier of frequency * 

To show how the total RF power of the conventional PSK/PM system 
gets distributed among the various components, we shall consider the rela- 
tively simple cases of a three-channel and two-channel system. 

The modulator for the three-channel conventional PSK/PM system 
is shown in Figure 1. The simplified block diagram for the corresponding 
receiver is also shown here. For the case of three channels, the received 
phase modulated RF signal is 

z(t) = ‘^JzP sin [tu^t + 9 jSj^{t) + 62^2^^^ 

where n(t) is assumed to be white Gaussian process of spectral density Nq 
and zero mean. 

Equation ( 4 ) can be rewritten in the form which shows the two orthog- 
onal components of the RF signal: 

z{t) = ^2F sinw^t cos [0jSj(t) + e2S2{t) + 0 jS^(t)] 

+ ^ZP coBCO^t sin[ 9 ^S^(t) + ^2^2^^^ 

= ^ 2 P sin cu^t cos 0 (t) + ^ 2 P cos w^t sin 0 (t) + n(t) . ( 5 ) 



Figure I. Conventional Three - Channel PSK /PM System 
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The trigonometric manipulations of the cos 6(t) and sin 6 (t) terms 
results in detailed expansion of ( 5 ) into its components. Thus* 

z(t) = ^|^ cos Sj(t) + ^ S^(t) + ^ Sj(t) - ^^^( 1 , 2 . 3 ) S^(t)S2(t)Sj(t)] 

+ yjz sin<n^t - ViPcm*2,3) S2(t)S3(t) - >^^(1,2) S^{t)S^{t) 


- Vp^^( 1 , 3 ) Sj^(t)S3<t)] + n(t) , 


( 6 ) 


2 2 2 

where Pj^ = P sin 0j^ cos cos 83 

2 2 2 

P^ = P cos e j sin 0^ cos 0 3 

2 2 2 
P3 = P cos 0 sin ©2 sin 83 

2 2 2 

P^ = P cos 0 cos ©2 COB 0 3 

P^„{ 1 » 2 , 3 ) = P sin^0, sin^G sin^0 

cm 1 2 3 

P^^(l, 2 ) = P sin 0^ sin 02 cos 83 

P ( 2 , 3 ) = P cos\ sin^0 sin\ 

cm 12 2 

P ^( 1 * 3 ) s P sin^0 cos^G sin\ 
cm 12 3 


~ Power in channel 1 

= Power in channel 2 

= Power in channel 3 

= Unmodulated carrier power 

= Third-order cross - 
modulation power 


Second -order cross* 
modulation powers • 

( 7 ) 


Besides the obvious significance of wasted RF power, the role of the 
cross -modulation terms can be examined further by considering the signal 
recovery process in the receiver. As shown in Figure 1, the receiver con- 
sists of a phase-locked loop which tracks the phase and frequency of the 
unmodulated carrier component and, at the same time, coherently recovers 
the multiplexed subcarriers. The subcarriers are then applied to data 
demodulators for baseband data recovery. Thus, as shown, the incoming 
signal z(t) of (6) is multiplied by the local reference cos ( t + 0 ), where 
^ is the residual RF carrier phase error. The result is: 
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y(t) = cos 0[-^^S^(t) + S^(t) - ^^^{1,2,3) Sj(t)S2(t)S3<t)3 
+ sin0[-^ - ^P^JiZ,3) S2(t)S3(t) - yjp^^^JT^ S^{t)S^{t) 

- yip (i^) SAt)S (t)] + n»(t) . (8) 

From (8), it becomes evident that,, if the frequency components of 
the second-order terms fall within the bandwidth of the lowpass filter, car- 
rier tracking capability of the receiver will be impaired. Furthermore, the 
third-order term's spectrum, depending on its nature, can degrade either 
the carrier tracking or data recovery, or both. The selection of the sub- 
carrier frequencies and the corresponding data rates thus becomes an 
important design consideration. 

Equations (7) indicate that, for the case of the conventional PSK/PM, 
the unmodulated carrier power cannot be fully suppressed, such as for 
Costas loop reception, without suppressing all channels except channel 1. 
Furthermore, as pointed out in [1], the conventional PSK/FM system does 
not make use of the cross -modulation terms which do contain information 
which could be utilized if properly recovered. In [l], it is also shown that 
for less than five channels, a system called Interplex (see also [2-5]) is 
more efficient. Since, in this appendix, we are concerned primarily with 
two- and three -channel multiplexing, we will therefore concentrate on the 
interplex and related systems. 

2. 2 Interplex Modulation 

With an interplex implementation of a PSK/PM system, the major 
cross-products of the conventional system are converted into a different 
set of modulation components by multiplying all of the secondary channels 
by the signal generated by the primary channel. The primary channel is 
generally the channel which carries the highest data rate and thus has the 
major portion of the RF power assigned to it. Thus, for an N-channel 
interplex system, the phase modulation term 0 (t) becomes 
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N 

0(t) = 0^ S^(t) + 5^ 0 ^Sj(t) S^(t) . (9) 

n =2 

• Figure 2 shows the implementation of a three -channel interplex 

modulator. The corresponding equation for the received RF signal is 

z(t) := ^/zF 8 in[&)^t+ 0 jSj(t)+ + 0 jS^(t)Sj(t)] + n(t) ( 10 ) 

where all of the terms are as defined earlier. 

Rewriting (10) in terms of its orthogonal RF components, one obtains 

z(t) = i^2P sinw^t cos [Oj^Sj(t) + ®3 

+ cosw^t sin[ 0 jSj^(t) + 02 S^(t)S 2 (t) + © 3 S^(t)S 2 {t)] + n(t) 

» <^ 2 P sin w^t cos 0 *®‘(t) + ^ 2 P coscu^t sin 0 *(t) + n(t) , ( 11 ) 

where 0 ^(t) indicates interplex time -varying modulation angle. 

Following the same procedure as for conventional PSK/PM, we 
expand ( 11 ) into its power term components. This results in 

z(t) = JT cosw t[‘^S,(t) -^jp (1,2,3) S,(t)S_(t)S,(t) 0T^S,(t)S (t) 

^ c ’ 1 i ' cin 1 6 3 ^ cm 1 2 


+ Sj(t)S3(t)] + sinw^t[^ - S2(t)S3(t) 

- ^|P^ S^[t) S 3 (t )3 + n(t) ( 12 ) 


where 


2 2 2 
= P sin cos 0^ cos 0 j 

- Power in channel 1 


*>2 

2 2 2 
= P sin 0 j^ sin 0 ^ cos 0 ^ 

= Power in channel 2 



2 2 2 
P sin 0 ^ COB 0 2 sin 0 ^ 

= Power in channel 3 

* 


P 

c 

2 2 2 
= P cos 0 ^ COS 0 2 cos 0 3 

= Unmodulated carrier power 
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P (1,2,3) 

= p 

. 2 

sm 0 , 

. 2 „ 
sm 0 _ 

. 2 ^ 
sm 0 _ 

cm 


1 

2 

3 

P (1,2) 

2 

= P cos 0 . 

. 2 

sm 0 _ 

2„ 

cos 0 _ 

cm 


i 

2 

3 

P (2,3) 

= p 

2 

cos 0 , 

* 

sm 0 _ 

• 

sm 0 - 

cm 


1 

2 

3 

P (1,3) 

= p 

2 

cos 0 , 

2 

cos 0 _ 

2 

sin 0 _ 

cm 


1 

2 

3 


Third-order cross- 
modulation power 


Second-order cross- 
modulation powers 

(13) 


The obvious significance of (13) is that letting 90° provides the 

following advantages; 

(a) power in all three channels is maximized, 

(b) unmodulated carrier is completely suppressed (Costas loop 
tracking can be utilized), and 

(c) all second-order cross -modulation terms are suppressed. 

For the case of 0^^ = jh 90 °, the signal portion of (12) can be written 
as 

z(t) = “^2P cos oi^t [Sj(t) COS0 2 cose^ - Sj^(t)S 2 (t)S^(t) sin 0^ sin 0^] 


- ^Jzp 


sino) t 
c 


[S^(t) 


sin 0^ cos 9 ^ + 


S3(t) 


cos 0 . 


sinOg] 


(14) 


Equation (14) is useful for providing a geometrical representation of the 
three-channel interplex PSK/PM signal. 

Figure 3 shows such representation for two of the eight possible 
vector states of the composite signal. The salient feature of the vector 
diagram of Figure 3 is that the cross -modulation term, although interfering 
with the amplitude of the main channel signal, keeps the envelope of the 
composite Signal constant. This can also be proved by showing that the 
sum of the squared bracketed terms of (14) is constant, regardless of the 
polarities of S^(t). 

The interference to channel 1 caused by the cross -modulation term 
depends on the relationship between the signal formats and the data rates 
of the three channels. If there is no spectral overlap between the three 


Constant envelope locus 



Figure 3. Geometrical Representation of a Three-Channel interplex 

PSK /PM Signal 
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channels, there will be no interference regardless of the relationship shown 
in Figure 3, 

On the other hand, if the bandwidth of S^(t) is such that it extends 
into the frequency regions occupied by S 2 (t) and interference will 

show up in channel 1 but it will not exceed the bounds set by power alloca- 
tions indicated by (13). Also, depending on the demodulation implementation, 
channel 1 can interfere with channels 2 and 3 if there is a spectrum overlap 
and carrier tracking error is not zero. This type of interference is dis- 
cussed in other sections of this appendix. 

The two-channel interplex system power allocations can be derived 
from (13) by setting 6, = 0, This results in 

Power in channel 1 

Power in channel 2 

Unmodulated carrier power 

Cross -modulation power ( 15 ) 

which, of course, is the same as the equation set (15) in [1]. Furthermore, 
setting “ +90^ for maximum efficiency reduces the two-channel inter- 
plex system to a quadriphase system with xinequal powers in the two orthog- 
onal components. Thus, (14) becomes simply 

z (t) = ^2P cos <y t [S. (t) cos 0 _] ~ -dzP sin<u t[S (t)sin0,] (16) 

where the subscript u indicates that we are now dealing with an unequal 
power division quadriphase modulation. The corresponding polar form of 
this equation is 




*’1 

2 

= P sin 0 ^ 

2 

cos 0 2 

^2 

« . 2 
= P sin 

• 2„ 
sm 02 

P 

c 

r, 2 

- P cos 0J^ 

2 

cos 0, 

cm 

2 

= P cos 0J^ 

2 

sin 0^ 


= cos [<o^t + S^(t) S^(t) 0 2 ] 


(17) 
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2,3 Equivalence of Two -Channel Interplex (e^ = +90°) and 
Quadriphase Modulation (QPSK) 

Consider now a signal which results from addition of two quadrature 
carriers* each of which is amplitude-modulated (DSB-SC) by a separate 
bipolar data sequence. The block diagram for generation of such a com- 
posite signal is shown in Figure 4. The resultant QPSK signal is given by 

f(t) = A a(t) cos at t - B b(t) sino> t . (18) 

c c 

where a(t) = binary sequence (+1) of channel 1 
A = amplitude of chaiuiel 1 
b(t) = binary sequence (j_l) of channel 2 
B “ amplitude of channel 2 , 

(The inversion of b(t) polarity has been introduced only for sign comparison 
consistency with interplex, ) It is customary to reduce (18) to an equivalent 
polar form, 

f(t) = Rcos(ot + 9). (19) 

c 

For this, we use trigonometric identity: 

cos (oj t + 0) = cos Qi t COB 8 - sinn) t sin 0 . (20) 

o o c 

Thus, 

f(t) = A a(t) COS a> t - B b(t) sinru t 

c c 

= R cos &) t [a(t) cos 6 ] - R sin6j t [b(t) sin 9 ] , (21) 

c c 

Note that (20) is of the same form as (16). Furthermore, for (20) to be true 

for all values of t, the coefficients of cos<u t and sinw t must be such that 

c c 

A a(t) = R a(t) cos 9 
B b(t) ~ R b(t) sin 0 




Figure 4. Two -Channel Quadriphase Modulator 
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Hence, (18) can be written as 

f(t) = V 


Z 2 
A + B COB 


<0 t + tan 
c 


■1 [Bb(t)1 l 
[Aa(t)]j * 


(23) 


The equivalence with two-channel (0^ = jt90 ) interplex is estallished by 
conversion coefficients 


R cos 6 2 

= ^ZP cos 0 2 

R sin e 2 

= ■^2P sine 2 

tane 2 



a(t) = Sj(t) 


b(t) = S (t) . 

Lr 


(24) 


Thus, we have shown that, for the two-channel case and antipodal modula- 
tion of the primary (high power) channel, the results of interplex and quad- 
riphase modulation are identical. 

From the standpoint of implementation, there are differences which 
must be taken into account. The interplex transmitter performs all the 
data combining at baseband and thus requires a phase modulator whose 
linearity and gain must be carefully controlled. This requirement is par- 
ticularly important for suppressing the cross -modulation and carrier terms. 
Quadriphase transmitter implementation, however, is based on combining 
individually modulated orthogonal components of the same carrier, a task 
which is relatively easy to implement with off-the-shelf components. Addi- 
tional similarities and differences are discussed in other sections of this 
appendix. 


2.4 Three -Channel PM Multiplexer /Modulator Implementation 

Based on Quadrature Carrier Summation 


A three-channel multiplexer /modulator which does not require a 
linear phase modulator, such as required by an interplex transmitter, can 
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be synthesized from the basic configuration of an unbalanced two -channel 
quadriphase modulator (6-8], Such a three-channel multiplexer /modulator 
is shown in Figure 5, This configuration can be considered as a hybrid 
between the interplex and quadrature carrier multiplexing. The similarity 
with interplex arises from the fact that angular modulation of the two secon- 
dary channels is proportional to the fractional power allocated to each of 
these channels. On the other hands the similarity with quadriphase modu- 
lation is due to the fact that two quadrature components of the same RF" 
carrier are used to synthesize the composite signal. 

As shown in Figure 5, the modulator consists of two balanced mixers, 
BM^ and which are supplied, respectively, by the two orthQgonal com- 

ponents of an RF carrier, BMj accepts the primary (wideband) data stream 
a(t) and thus provides a constant envelope biphase modulated carrier to the 
summing network, BM^ accepts a linear sum of baseband secondary' data 
stream b(t) and a square wave subcarrier biphase modulated by the third 
data stream c(t). The output of BM^ is then a double sideband, suppressed 
carrier whose amplitude and phase vary according to the instantaneous sum 
of the two signals at its modulation input terminal. 

The outputs of BM^ and BM^ are summed and applied to an amplitude 
limiter. In this limiter, the amplitude modulation of BM output is converted 

Ca. 

into a phase modulation which adds to the antipodal modulation of the main 
carrier vector. The addition is such that it provides an orthogonal carrier 
whose amplitude and phase vary in accordance with the output of BM_. The 
amplitude limiter, however, keeps the amplitude of the composite signal 
constant. This process can be understood by referring to the vector dia- 
gram shown in Figure 6. Before limiting,the composite signal consists of 
an orthogonal addition of (A'| and IB' + C*) . where A', B' and CV are the 
amplitudes of the three data streams prior to limiting. After limiting, the 
amplitude of the A vector gets modulated, the nature and the degree of mod- 
ulation being determined by the instantaneous polarities of the B and C 
vectors. This modulation is similar to the third-order cross -modulation 



Subcarrier 


Figure 5. Three-Channel PM Multiplexer/Modulator 

Implementation Based on Quadrature Carrier 
Summation 
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generated by a three -channel interplex signal. The magnitude of this modu- 
lation can be computed by considering the terms: 


A = R cos 
max 


A . = R cos 

mm 


A = R cos 
ave 




(25) 

(26) 


(27) 


where R is the magnitude of the constant envelope of the RF signal out of the 

amplitude limiter and A » A . t and A are as indicated in Figure 6, 

max mm ave 

The cross-modulation term’s amplitude, normalized to the average 
amplitude of the A vector, is approximated by 


A - A . 
max mm 


ave 


2 A 


ave 


cos 


\A7 


(28) 


2 cos tan 


This term is e-ssentially equivalent to the corresponding term obtained from 
(13) for the interplex system. The latter term is 


A (Si) 



cm 


(1.2,3) 


tan 0 tan 0 „ , 


(29) 


where Si stands for the amplitude of the cross -modulation term superimposed 
on the primary data stream Si(t) of the three -channel interplex system with 
= i90° , Note that both (28) and (29) are third-order cross -modulation 
terms and, as shown below, these terms are approximately equal provided 
the primary/ secondary channel power division ratios are 6 dB (4 to 1) and 


up. 
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2.5 Demodulation 

For both the interplex and the hybrid phase multiplexing, it has been 
shown that the antipodal modulation should be used for the primary channel, 

1. e. , the channel which carries the highest data rate and thus contains the 
largest fraction of the total RF power. Ff r such power allocation, com- 
bined with antipodal modulation, a Costas loop demodulation can be used 
for carrier tracking and data demodulation. As described in [7] and in 
Appendix B, such a demodulator, although not optimal, provides good per- 
formance and yet is relatively simple to implement. Such a demodulator 

is shown in Figure 7. With the carrier frequency and phase lock established 
within residual phase tracking error (p , the demodulated primary data stream 
appears at the output of the balanced mixer BM^. The secondary data streams 
which are modulated on the orthogonal component of the carrier appear at the 
output of BM^* These signals are then processed in their respective data 
recovery circuits and the data streams recovered. 

2. 5. 1 Two -Channel Quadriphase Input Demodulation 

Assuming that the frequency lock is perfect, i. e. , (U ' = , and that 

c c 

only the residual carrier phase error is present ( <p «45°), the expres- 
sions for the outputs of BM^ and BM^ in Figure 7 can be obtained for the 
baseband terms by neglecting the terms of the appropriate trigonometric 
expansions* Thus, the baseband output of BM^ is 

y (t) = [A a(t) cos tu t + B b(t) sin t + n(t)] • -v/i” cos {to t + <^ ) 

X c c ’ c 

= — “ [A a{t) + n (t)] cos<^ + — ~ [B b{t) + n (t)l simi (30) 

yjZ ^ ® 

' V ' V ' 

Channel 1 data and noise Channel 2 crosstalk and 

noise 

Similarly, the baseband output of BM^ is 



Carrier 

Error 

Detector 


To data 
processors 


Figure 7. Costas Loop Demodulator for Interplex or 

Quadriphase Modulation 


•I 
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Y2<t) 


= — ^ [A a(t) + n (t)3 eim^ 

- - - J 


[B b(t) + n (t)] cos ift . 

VT 


Channel 1 crosstalk and 
noise 


Channel 2 data and noise 


(31) 


The terms n (t) and n (t) in both (30) and (31) are the projections of the band- 

C B 

pass noise in the cosine and sine directions of the composite modulation 
process. 

From (30) and (31), it is evident that the penalty paid for using only 
the major component of the quadriphase signal is the existence of the cross- 
talk between the channels when the tracking error <p is not zero. Also, as 
is shown below, the power available for carrier tracking is also reduced, 
but this reduction is not significant (less than 1 dB) if power division ratio 
between the two channels is 6 dB or more, 

Liet ^1^(0 l>o thi. ratio of the power in the kth channel to the crosstalk 
caused by the f channel data. Also let B^ and B^ be the bandwidths (to 
the first null) of the A a(t) and B b(t) signals, respectively. Now if we 
assume that Bj > B^, which is generally the case, we obtain, neglecting 
the noise terms in (30) and (31), 


D^(2) 


A^[a(t)l^ 

B^[b(t)f 



(32) 


D^d) 


B^fb(t)1^ / / -cos0^ 
A^[a(t)3^^®2'^ 



(33) 


where and are the powers of the respective carrier components and 
<f) is the residual carrier tracking error. 

Note that in (32), the bandwidth ratio does not appear because as 
long as is located within the spread of B^, all of the power P^ contributes 
to the crosstalk. However, because the rate of b(t) is less than that of a(t), 
the crosstalk appears as a bias which shifts the mean of the A a(t) signal. 


21 


On the other hand, the crosstalk has a general nature of filtered 

pseudorandom signal and thus is a noise -like process with statistics being 
determined by the properties of the a(t) and b(t) data streams. 

The carrier tracking error is provided by multiplying y^^ and in 
the carrier error detector M and lowpass filtering. The output of the 
carrier error detector is thus: 


e(t) = y,y, = ~[A^ +2Aa(t)n (t)-2Bb(t)n (t)+n^(t) -n^(t)] sin 2<6 

L Lt ^ C 6 C S 

- -7 [AB a(t)b(t) + B b(t)n (t) + A a(t)n (t) + n (t)n (t)] cos Z(f> . 

Z C 5 c 8 


Note that, in the absence of noise, (34) reduces to 

2 _2 


e(t) = 


A - B 


. „ , AB a(t)b(t) „ , 
sin 2<p ~ r cos 20 , 


(34) 


(35) 


The first term in (35) is the carrier tracking DC term, which is reduced 
from its full possible value by not utilizing the power in the second channel. 
This reduction is negligible, however, for channel power division ratios of 
6 dB or more.v The second term represents wideband crosstalk which enters 
the carrier tracking loop and is a function of spectral properties of a(t) and 
b(t). The effect of this term on the carrier tracking performance is, further- 
more, the function of the carrier tracking loop bandwidth. In general, the 
faster the data streams a(t) and b(t) and the narrower the loop bandwidth, 
the smaller is the interference. 

A typical situation may involve Manchester type data for both data 
streams, resulting in spectrums shown qualitatively in Figure 8 . The con- 
volution of the two spectra results in the cross-product spectrum shown in 
Figure 9. From this figure, it is evident that, in the region of the low fre- 
quencies, where the carrier tracking bandwidth lies, the spectral density of 
the cross-product signal is relatively flat. 

The actual magnitude of the spectral density in the low frequency 
region is given in Appendix B by 
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S(0) = 


2 T 


(36) 


B 


where S(0) = spectral density of a(t)b(t) at f = 0 
Ta = symbol time of a(t) 

= symbol time of b(t) . 

The e 3 q>ression for the signal-to-crosstalk power in the carrier tracking 
loop, for the noiseless case and both data streams Manchestered, is thus 


S(0) B, 


, (A^ - 
2 A^B^ 


R. 




for Ra^ 2R^ 


(37) 


where R. = Channel 1 data rate 
A 

R„ = Channel 2 data rate 

A = Channel 1 signal amplitude 
B = Channel 2 signal amplitude 
Bj^ = Carrier tracking loop noise bandwidth . 

For the case where a(t) is biphase modulated and only b(t) is Man- 
chestered, the spectral density in the low frequency region is expressed 
in Appendix B by 

T 

S(0) 




(38) 


and, correspondingly, the signal-to-crosstalk power in the carrier tracking 
loop is 


S(0)B, 


. (aL^. 


R 


4aV 


( 39 ) 


where all the symbols in (38) and (39) are the same as previously defined* 
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2, 5.2 Two -Channel Interplex Input Demodulation 

A two -channel, quadrature interplex signal can also be recovered 
by the Costas demodulator of Figure 7. The corresponding’ expressions 
for the y. and y signals are obtained by multiplying z (t) of (16) with the 

* w \1 

cosine and sine components of the locally generated replica of the RF car- 
rier. Thus, one obtains 

[Sj^(t) cos 6 2 cos^ + 2 ^ sin^ (40) 


^ V, 


Channel 1 data and' noise 


Channel 2 crosstalk and noise 


and 


= \jp [s^(t) cos 0 2 + ^ ^ ® 2 ^ COS0 (41) 




j ^ 


Channel 1 crosstalk and noise Channel 2 data and crosstalk 


For a noiseless case, it is a simple matter to derive from (40) and 
(41) the data crosstalk equations similar to (32) and (33). These are 

2 


V ■ = 


2 2 , 
cot 0 2 cot <p 


(42) 


V. ■ (S:)(^K^) -(St) 


2 2 

tan0 cot <f> . (43) 

& 


'W'e can also derive the equivalent expression for the carrier tracking error, 
which for the two-channel interplex is 


e(t) = yj^y^ = ^ ^®2 ^^2 

Note that (42), (43), and (44) can be derived directly from (32), (33) and (35) 
by use of the conversion coefficients of (24). 
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2, 5. 3 Three-Channel Interplex Demodulation 

When the input to the Costas demodulator in Figure 7 is a three- 
channel interplex signal defined by (14), the corresponding Yj^(t) and 
signals are as shown below. The Yj^(t) signal is 

y^(t) = [Sj^(t) cosQ^ cos ©2 " cos^ 

^ . , ( . 1 . — , — t 

Channel 1 data Cross -modulation Noise 

+ ^ ®2 ^ ®2 ®3 ^ 

V I I » > * ' 

Channel 2 data Channel 3 data Noise 

crosstalk crosstalk 


Equation (45) shows that, even if the carrier phase tracking is perfect (i, e, , 
(p- 0), a small third-order cross -modulation term contaminates data 
output. This is not the case for the output where data crosstalk occurs 
only for ^ 0; 

y 2 (t) = JP^ sine 2 cos - S^Ct) cos 0^ sin 0 ^ - "^(t)] cos<^ 


Channel 2 data Channel 3 data Noise 


^[S 


j(t) cos 0^ cos 


Si(t)S 2 ,(t)Sj(t) sin sin 0^ 


+ n (t)] sin^ 
c 


Channel 1 data Cross -modulation Noise (46) 

crosstalk 


Let D, ( f ,m) be the ratio of the power in the kth data channel to the 

lx 

crosstalk caused by either the f ^ or mth channel, or both. Also let 
be the bsindwidth (to the first null) of the nth channel. Equations (45) and (46) 
can now be used to define the D^( £,m) ratios. Thus, we have for the higher 
rate data channel: 
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2 2 

Dj^( 2 ) = cot 0 ^ cot 0 
Dj^(3) s: cot^G^ cot^0 
D (2,3) = cot^e cot^e, 

X b 

. and, for the secondary, lower rate data channels; 

2 2 

D^U) = ^ tan 6 ^ cot 0 
3 

2 2 

03 ( 1 ) = — tan 83 cot (p 

®13 2 2 

= " 5 ~“ cot 83 cot <p 
2 

B 

D2(1,2) = cot^e 2 Cot^^b . 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 


Note that in (50) through (53), the data bandwidth ratios are included. These 
ratios, however, provide only the first-order approximations, particularly 
in (52) and (53) where batndwidths determined by convolving 

B, with B_ and B with B , respectively. Therefore, to obtain an accurate 
bandwidth ratio improvement factor, each specific combination of data rates 
and formats must be examined in detail. This is also an area where a great 
deal of optimization can be obtained for any particular link by selecting the 
appropriate modulation formats for the component data streams of the com- 
posite multi-channel signal. 

The carrier tracking error signal which is applied to the loop filter 
is, prior to filtering, the product of (45) and (46). Thus, for a three-channel 
interplex signal, we have 


I 


i 
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e^(t) 


= yi(t) y,(t) 


= — cos 26 ^ cos 2 0^ sin 2 <f> 


^ fi 

111 'M 

14 o ™ 
M rt S, 

UHw 


i- 

O 


- ^ [S2(t)S3{t) sin 2 02 sin 2 0^] sin 20 

- ^ [Sj^(t)S2(t) sin 20^ cos 263 + S^(t)S3(t) cos 20^ sin 283] cos 20 

- P[S,(t) cos0„ cos 0_n - S,(t)S_(t)S,(t) sin 8, sin8, n 

I L 3 B L 5 C as 


I 


CO 


+ S-(t) sin0 „ cos 0 _ n + S„(t) cos 0 _ sin 0 _ n ] cos 20 

Ci a aCa C ac 


a 

to ** 


u er 


+ PfS, (t) cos 0 , cos 0 _ n - S_{t) sin 0„ COB 0^ n 

1 ^ D Q £, C aS 

- S_(t) cos 0_ sin e n - S,(t)S,(t)S (t) sin 0 sin0. n ] sin 20 
5 L asi^a ^ aC 


- P (n n ) cos 20 
c s 


.P, 2 2 ^ . 

+ - (n^ - n^) sin 20 . 


0) 

n 

•rl 

0 

a 


( 54 ) 


A similar eiqjression can be obtained for the hybrid multiplexer signal 
by substituting the appropriate magnitudes of the A, B and C components. 
Equation ( 54 ) and its equivalent for the hybrid multiplexer can then be used 
to obtain carrier tracking performance as a function of power division ratios, 
data formats, data rates, and the P/N^ ratio, 

2.6 Numerical Example for a Three-Channel System 

L»et us consider a numerical example for a three-channel system 
which may be typical of a space communication link. The nominal data 
rate bandwidths, formats auid power divisions may be as shown in Table 1 . 
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Table 1. Bandwidth, Data Formats, and Power Assignments 
for the Numerical Sample 


Channel 

F-ndwidth 

(MHz) 

Data Format 

Power 

Assignment 

{%L 

1 

100 

Biphase, baseband 

80 

2 

4 

Bi-0-L, baseband 

18 

3 

0.4 

Bi-<^-L, subcarrier 

2 


Note that, in this example, we arbitrarily have assigned more power 
to the secondary channels (2 and 3) than would be required if power division 
had been carried out based on equal error rates (same E^/Nq) in all three 
channels. Also, to minimize the cross-talk between channels 2 and 3, the 
baseband data of the latter channel is assumed to be riding on a squarewave 
subcarrier whose frequency is well above the high end of channel 2 data 
baseband spectrum. On the other hand, the frequency of the channel 3 sub- 
carrier should be low enough to allow at least the 3rd and 5th harmonics to 

« 

pass through the RF channel without significant attenuation. 

2. 6. 1 Interplex Implementation , 

Assuming that a three-channel interplex system is used, the modu- 
lation angle assignments are as follows: 

9 ^ = 90 ° (±) 

0^ = tan'^JP^TpJ" = tan‘^ = tan“^0.47) £ 25.4° (±) 

0j - tan"^ ^2/80 s tan‘^0.16) S 9*^ (±) (55) 

From these angle assignments, the data -to -crosstalk ratios Dj^(j?,m) as a 
function of carrier tracking error (f) can be computed using (47) through (53), 
These are plotted in Figure 10, From the plots, it is evident that Dj^(2), 

D«(l), and D.(l) are affected most by high values of 0 . But, because the. 
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channel 1 data rate is much higher than the channel 2 data rate, Dj^{2) in 
f this case is the ratio of the S^Jt) signal power to the time-varying bias due 

to data. The same is true of Dj^(3), Consequently, at least in principle, 

both of these ratios can be improved by the use of coarse estimates of S_(t) 

^ and S^(t}. 332 ( 1 ) and on the other hand, are determined by filtered 

Sj^(t) data stream and thus their crosstalk noise floor is a fluctuating signal 
. with statistics approaching that of a Gaussian process. Thus, these ratios 
can be improved by optimal filtering of S (t) and S>(t) data streams. Of 

Ct J 

course, once bit synchronization is established, such optimum filtering will 
be the integrate-and-dump detectors for the S_(t) and S (t) data streams. 

2, 6.2 Hybrid Modulator Implementation 

To achieve the same power division between the three channels as 
was assumed for the preceding example for interplex, the relative amplitudes 
of the respective data streams will be as follows: 

A' = 1,0 

B* = 0,474 ) normalized with 

C’ = 0,158 ) respect to A* (56) 

Since we are primarily interested in comparing the third-order terms (i. e,, 
the terms independent of carrier tracking error , we compute the A A 
term by substituting values of (56) into (28), This results in 

A A = 0,081 . (57) 

ave 

The equivalent D^(B, C) ratio for the hybrid implementation is therefore 

D^(B.C) = 20 1og^Q(l/AA^^J 

= 20 logjQ (1/0.081) = 21.8 dB. (58) 

^ ■ The corresponding term, i. e., Dj(2, 3) for the three-channel interplex sys- 

tem is 22. 5 dB, which is only 0,7 dB above the value shown in (58). Con- 
sidering the absolute values of the D (B, C) and D,(2,3) ratios (*^20 dB), 

A 1 
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this is a small difference in view of potential implementation simplicity of 
the hybrid modulator /multiplexer, 

3. 0 CONCLrUSIONS 

It has been demonstrated that either the interplex or the hybrid 
phase -multiplexing technique for transmission of three independent and 
asynchronous digital data channels offers a significant performance advan- 
tage over a conventional PM subcarrier design. The hybrid phase - 
multiplexing technique differs from interplex primarily in terms of 
implementation. Whereas interplex involves data channel combining at 
baseband and utilizes a linear phase modulator at IF, the hybrid technique 
involves data channel combining at IF and utilizes quadrature carrier mod- 
ulation. Since the performance characteristics of these two approaches 
are essentially identical, it appears that the selection of a teclinique for 
any particular application should be dictated only by equipment design 
considerations. 
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1.0 INTRODUCTION 

Signal waveforms for the Orbiter Ku-band relay communication 
links have been described in [3], along with the link power budgets. Data 
multiplexing, encoding, and modulation schemes for both the forward and 
return links have also been described in general terms. 

We consider various candidate receivers for the unbalanced QPSK 
signal format. The basic modulator configuration is shown in Figure 1 
for both the return and forward links. In all cases the transmitted signal 
is 

S(t) = A a(t) cos cj t + B b(t) sin uj t (1) 

c c 

where a(t) represents the encoded data of the high rate, or I, channel and 
b(t) represents the encoded data in the low rate, or Q, channel. 

Referring to Figure la, in the I channel, a(t) is a sequence of ±l's 
of NRZ data which has symbol rate up to 100 Msps. It represents rate 
R- 1/2 convoiutionally encoded data which has variable rate from zero up 
to 50 Mbps, In the Q channel, b(t) is either at the fixed rate of 192 kbps 
of Bi-<j>-L signals, or it is variable rate from zero to 2 Mbps of Bi-(|)-L 
signals. The data rates for the two links are completely independent and 
may be either synchronous or asynchronous. 

I 

The proposed power division for the return link is 80 percent and 
20 percent for the high data rate and low data rate links, respectively. 

With this power division ratio of 4:1, the power available to the high data 
rate link will be reduced by about 1 dB with respect to the total power. 

The purpose of unequal power ratio is to cause the symbol energy, and 
therefore the error rates, in the two channels to be closer to the same 
value. 


I 
















The present design calls for a fixed 4:1 power ratio in both the 
forward and return links. From an error rate point of view, this over- 
powers the low data rate channels. However, the maximum to be gained • 
by placing all the power in the high data rate channel is only 0. 97 dB. 

This power ratio of 4:1 is considered satisfactory, and we shall emphasize 
this power ratio throughout the present study. 

In Section 2.0, the fourth-power receiver tracking loop is con- 
sidered. When the power division in the unbalanced QPSK signal is 4:1, 
it is shown that this receiver tracking loop is an unattractive choice. 

In Section 3.0, the Costas loop is evaluated as a candidate 
receiver for both the forward and return links. In this appendix, we restrict 
attention to design of the Costas loop receiver which is tracking primarily 
on the high data rate signal. In sections 3. 1 and 3.2, it is shown that the 
Costas loop quadriphase demodulator in Figure 6 performs with an RMS 
phase error of less than 10 degrees over a variety of tracking loop band- 
widths and over a variety of rates in both the high and low data rate chan- 
nels. The performance curves are shown in Figure 7, where it is seen 
that this performance is attained for a power division of 70/3 0 up to 100/0. 

In sections 3.3 and 3.4, similar results are obtained for the for- 
ward link Costas loop shown in Figure 5. 

In section 3.5, the return link Costas loop is examined when the 
rate in the high data rate link goes to zero, but the power remains on. 

The tracking is then carried out primarily on a monochromatic tone and, 
as a result, performance improves with respect to that reported in the 
earlier sections of 3.0. 

It should also be noted that, if the power in the high data rate link 
is turned off, tracking is then carried out only on the low data rate link. 
Without any alterations, the Costas loops in Figures 5 and 6 will shift 
their reference by 17/4 radians, thereby sending the low data rate signal 
into the high data rate demodulator. This problem can be easily overcome, 
but a better solution may be to not remove the power in either of the data 
links without removing both. 


An additional candidate receiver is the four -phase Costas loop 
receiver [4], This receiver requires significantly more equipment for 
implementation with little anticipated improvement in performance. This 
receiver would also exhibit phase a^nbiguities of multiples of ir/4 radians, 
so that an algorithm would have to be developed to process the proper data 
links into the appropriate demodulator. For these reasons, this receiver 
was not considered in this feasibility study. 

In the last section, some approximate error rate computations 
are presented that show that performance is within a few tenths of a dB 
of ideal performance, 

2.0 FOURTH-POWER RECEIVER TRACKING LOOP 

A block diagram of the fourth -power tracking loop is shown in 
Figure 2. The signal at the input to the fourth -power device is 

y(t) = s(t) + n{t) , (2) 

when s(t) is given by (1) and the narrowband noise n(t) is given by 

n(t) = [n (t) cos a> t *i" n (t) sin w t] . (3) 

CCS c 

The independent baseband noises and n^(t) each have one-sided power 

spectral density equal to that of n(t), namely, watts/Hz. The bandwidth 

of n (t) and n (t) is B _/2, where is the bandwidth of the IF filter H(p), 
c s Xr II? 

as shown in Figure 3, 

The difficulties associated with the fourth-power tracking loop 
can be demonstrated without considering the additive noise when the power 
division is 4:1, 

With signal only, therefore, the output of the fourth-power device 
in Figure 2 is ’ 
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4 4 

y (t) = s (t) 

4 4 3 3 

= A cos u) t + 4A B a{t)b{t) cos u t sin w t 
c c c 

+ 6A^B^ cos^ cj t sin^ a> t + 4AB^ a(t)b(t) cos w t sin^ w t 
c c c c 

^4.4^ 

+ B sin (j t 
c 

, ^ ^ ^ J cos (Zco^t) + 3\^— j 

1 O ^ 

+ “a(t)b(t)jA B[ sin (4ca t) - sin (2 <jj t)] - AB [sin (4cij t) + sin (2w t]}- 

£i C- C C C 

+ ~ (A^B + AB^) sin (2u) t) + |* A^E^[l - ^(l + cos (4w t))] . (4) 

C C Z ' 4 \ c / 

After passing the BPF centered at 4f^ in Figure 2, the result is 


z(t) = 


A + B 


/ a2„2 

- 6A B 


cos (4.t)+2^^^(A^B 

c ^ 


AB ) sin (4io^t) . (5) 


The first term in (5) is the tracking term, and the second term 
is a self-noise term. In terms of the average power in each link, namely, 

= A^/2 , Pg = B^/2 , (6) 

2 (t) can be expressed as 


z(t) = i + p| - 6P^Pg 


] cos (4uj^t) + 2a(t)b(t) J. 


P P 
A B 




( 7 ) 


The tracking term contributes a spectral component at The 

self-noise term contributes a spectral density at 4 w^t. The level of this 
spectral density within the receiver tracking loop is dependent on whether 
the forward or return link is considered, since the waveforms are different. 
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in (7). 


The amplitude of the tracking term is the coefficient of cos (4(o t) 

c 

The power in that term is 


track 


= -[P^ + 

8 ^ A B 




( 8 ) 


The self-noise power in the tracking loop is the equal to the spectral density 
of the second term in (7) at 4f^ times the tracking loop bandwidth. This 
self -noise power can bo expressed as 

PSN = <’> 

where B^^ is the tracking loop noise bandwidth and S(0) is the two-sided 
spectral density of the signal a(t)b{t) at f = 0, In Addendum A, we show 
the spectral density of a(t) b(t) at f = 0 is 


S(0) 




( 10 ) 


when a(t) is NRZ and b(t) is Bi-(j)-Ij, and 


S(0) 



( 11 ) 


when a(t) and b{t) are both Bi-c|}-L. In (10) and (11), T . and T are the 

A B 

symbol times of a(t) and b{t), respectively. 

In order to obtain a preliminary evaluation of the performance of 
the fourth-power tracking loop, let the total power be designated as 
so that 


P 


A 


aP 


T 


Pg = (l-a)P^ . 


( 12 ) 


Upon substitution of (12) into (8) and (9), the ratio of tracking power to 
self -noise power within the tracking loop can be written as 


t 


I 
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Ptrack ^ J_ l^A " Pb ~ 

^SN. PA*“Bt^A- 

1_ N 

32 D[2S(0)] 

where 

N = [a + (1 -a) - 6o'(l -a)]^ 

and 

D = a'(l-Q’)[o'-(l-a)]^, 


(13) 

(14) 

(15) 


The parameter a represents the fraction of total power- in the high data rate 
link, which assumes values over 0.5 £ o' < 1, The numerator N in (13) and 
(14) represents the variation in tracking power as a function of the power 
division and is shown in Figure 4 as the dashed line. The denominator D 
represents the variation in self-noise spectral density as a function of the 
power division. The denominator is also shown in Figure 4, where it can 
be seen that D reaches a maximum and N has a minimum of 0 near f( - 0, 85, 
The signal power to self -noise spectral density, 


N 

32D 


= 2S(0) 


P 


track 



( 16 ) 


is also shown in Figure 4, where it is seen that the worst power division 
values are within [0, 75, 0, 95] . We see that the fourth power tracking loop 
makes very ineffective use of the total available tracking power when the 
power division is in the vicinity of 0,8, An additional disadvantage of the 
fourth-power tracking loop is that the stable lock points occur at multiples 
of tt/ 4, thereby allowing the possibility that the two data channels become 
interchanged. This could be detected and corrected in the system design, 
but only at the expense of additional complexity. 

With these results, we conclude that the fourth-power tracking loop 
should not be recommended for the demodulation of an unbalanced 4:1 QPSK 
signal. 
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Figure 4, Tracking Power, Self Noise Density and Ratio of Tracking Power to Self Noise 
Density for the Fourth Power Tracking Loop. 
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3.0 COSTAS RECEIVER 

In this section, we consider the Costas receiver for demodulating 
an unbalanced QPSK signal. Block diagrams of the receiver as shown in 
Figures 5 and 6 for the forward and return links, respectively. 

The received signal is as described in (1) and (2), where we elect 
to model the noise [6, 7] as 

n(t) = [n (t) cos (co t +<j)) + n (t) sin (w t + (j>)] {17) 

c c s c 

where n (t) an i (t) have the same statistics as described in (3). 
c s 

We initially develop a generalized description of the receiver 
operation, then consider several specific cases. 

Referring to Figures 5 and 6, the reference signals are given by 


r(t) = ''/z' cos (w^t + (ji) 

and 

r'(t) = ^/2' s in ( u> t + 4> ) , 

c 


(18) 

(19) 


Under the assumption that the IF filter preceding the Costas loop 

receiver is sufficiently wide to pass both channels of the received signal 

without distortion, then the baseband outputs, e (t) and e (t), of the mixers, 

c s 

neglecting the 2f^ terms, are given by 

ejt) * y(t) 

= (n/^/Z)[A a(t) cos 4> - B b(t) sin <j)] + n (t) ( 


''"’Ivideo 


= (<72 /2) [A a{t) sin <}i + B b(t) cos (j> ] + n (t) 


The Return Link Costas Loop Receiver 


In considering the i*eturn link, we first consider the bandwidth of 


the LPF to be 
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= 100 MHz . (22) 

thereby allowing the passage of the signal portions of e (t) and e (t) without 

c s 

significant distortion. This is equivalent to saying the tracking loop is 
going to function primarily on the high-rate, high-powered channel* The 
dynamic error signal e(t) is then given by 

e(t) = e (t) e (t) 

C S 

= ( ^ sin Z* 

+ a(t)b(t) cos 2(t) 1 Self-Noise 

+ n (t) n (t) } Noise 

+ a(t) sin 4> + JPg b(t) cos <j) ] ^ 

>S X N 

+ n^(t) a(t) cos 4> - JP^ b(t) sin <})] ' 

where, in addition to the desired tracking signal, self -noise, noise, and 

(signal X noise) terms are generated, which we designate as the equivalent 

loop noise, n (t). 

eq 

Observation of the tracking signal •• pt^o- 

vides the following: 

(a) The stable lock points are at ^ - kir, with the result that the 
two data links cannot be sent to the incorrect data demodulator# The high 
data rate link will always go to the high data rate demodulator, and the same 
for the low data rate link. This assumes that the power in both channels 
is always present, which is the most likely case. 

(b) Although it is expected that the power in both channels will 

always be *'on'‘ even if the data rates are reduced to zero, should the power 
in the high data rate channel fail, i*e. , then the sign of the tracking 

signal is reversed and the stable lock angles are at -n/Z ± kir* If this 



Tracking 

Signal 
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occurs, then the low data rate link will be processed into the high data rate 
demodulator. 

(c) If the power in the low data rate link goes to zer j, i. e. , 

Pg-* 0, the stable lock angles remain unchanged, so that the high data 

rate link continues to be processed into the high data rate demodulator. 

The statistics of the equivalent noise, n , (t), are obtained by 

eq ' 

assuming that the tracking loop bandwidth is much smaller than The 

phase error t|i{t) will then be slowly varying with respect to the other 
random processes present in e(t). The first and second conditional moments 
are given by 

Efn^^(t)|4)] = 0 (24) 

and 

= lyr) R^(t) cos^ (2 4>) 

+ B^(t) + Vt)[P^ R^(t) + P^R^(-r)] , (25) 


where R (t) and R, (t) are the autocorrelation functions of a(t) and b(t), 

S U 

respectively, as discussed in Addendum A. 

The spectral density of is essentially flat over the bandwidth 

of the tracking loop, so that the significant statistic when considering track- 
ing loop performance is the spectral density at f = 0. The spectral density 
of the terms in (25) can be shown [5] to be 


2 2 ^0 ®1 
Rj^(t) S(0) = (Nq/ 2)^ df = 

V^)[P^RJt) + P3R^(t)] - S(0) =. ^(P^ + Pb) 
P^Pr \(t) \(t) cos^ (2 4*) S(0) « P^P^ T^ (l - 



(26) 

(27) 

(28) 


I 


where, in (28), we have imposed the upper bound, cos (2<|)) < 1 . In addi> 

tion, we have used the fact that the spectral density of at f = 0 

is given by = 2T^[l - (T^/T^)] for the return link which consists of 
NRZ waveforms in the high data rate channel and Bi-cjj-Ij waveforms in the 
low data rate link. This result is derived in Part I of Addendum A. 

Combining these results, the two-sided spectral density of n^^(t) is 


No B, 

S ,0) * + 

»eq 2 2 




We shall base our performance analyses on a linearized assumption# 
We shall see that^ at the signal-to-noise ratios which are anticipated, the 
Costas loop is well within its linearized region* When such is the case, the 
variance of the phase error is given by [7] r 

n -2 1 z 1 , , ,2 

■’4’ = 4 "2<|. =47 

where p Is the signal-to-noise ratio in the Costas loop, which is given by 


Nq B 
*^eq L 

where P is the available power in watts for tracking, B is the one-sided 

C Lj 

equivalent noise bandwidth of the tracking loop in Hz, and Nn is the one- 

^eq 

sided spectral density of the equivalent additive noise in the loop in watts/Ha, 
In our receiver, Nq is given by (29), and from (23); 

Upon substitution of (29) and (32) into (31), the Costas loop signal- 
to noise ratio is given by 


<^a-Pb> 

I 2Pa^b 


(■ ■ k)\ 
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Using the power division parameter ce , as given by (12), 


P = 




(2 g - 1) 


N B 2P g(l-g) 

1 + — — - + = T U 

P N A\ T 

T 0 B‘ 


('-aJ 


(34) 


The performance of the return link receiver using a Costas tracking 
loop and tracking on the high-powered channel is given by (30) and (34), 

3*2 Performance Computation - Return Link 

The performance of the return link Co.stas loop using the high data 
rate channel is shown in Figure 7. The standard deviation of the phase error 
is shown versus the power division parameter a. Two sets of curves are 
shown* Set (1) shows performance for the symbol rates = 384 ksps and 
192 ksps* This is a worst set of rates from the point of view of maxi- 
mizing the standard deviation of the phase error. As shown, the two curves 
are for loop bandwidths of 5 kHz and 10 kHz, As cv ^ 0* 5, the Costas loop 
is unable to track and tr n) . As cv 1, the signal becomes single channel, 
the high data rate channel, and the self-noise goes to zero* The effect of 
the self -noise can be determined by comparing the standard deviation at 
some a to that at a =1. The increase is due to the effects of the self- 
noise, and the fact that the available tracking power is also decreasing. 

The second set of curves in Figure 7 is for the best set of rates 
in the return link. The improvement due to the resulting reduction in self- 
noise at the higher symbol rates is seen to be approximately an order of 
magnitude in (r . The points on the curves corresponding to a = 0* 8, a 
4-to-l power division, have been emphasized by dots, since this is the 
most likely choice of power division, 

A summary of the various parameters used in the forward link 
and return link computations is shown in Table I. 






TABLE I 


FORWARD AND RETURN LINK PARAMETERS 



F orward Link 

Return Link 

Channel A 

0 to 1 Mbps 

0 to 50 Mbps 


R = 1/3 

R = 1/2 


0 to 3 Msps 

0 to 100 Msps 


Bi-(j>-L 

NRZ 

Channel B 

72 kbps 

1 92 kbps 


F ixed 

Fixed 

or 


Bi-«|>-L 

2 Mbps 


No coding 

Bi-(j)-L 

One-Sided 

N' = T ^/T 
0 A^ B 

. / '^a\ 

Self-Noise 


Spectral Density 


\ B/ 

at f = 0 

Bi-c|)-L and 

Bi-(|>-L and NRZ 


73. 8 dB-Hz“^ 

86.8 dB-Hz 

When Tracking 
on Channel A 

B^ = 6 MHz 

B^ = 100 MHz 

(Curves 3 and 4) 

(Curves 1 and 2) 

When Tracking 
on Channel B 

= 72 kHz 

B^ = 4 MHz 


Includes all cases except symbol sync. 


Reference B. H, Batson power budgets. 
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Forward Ltinlc Costas Loop Receiver 


We now consider the forward link Costas loop receiver. The 
development is similar, with the fallowing exceptions: 

(a) The two data links are both Bi-<j>-L waveforms, with the result 
that the self-noise spectral density at f = 0 is given by 

No = 


as given in Table I and derived in Addendum A, Part II. 

(b) The IF bandwidths are set at = 6 MHz, which will pass the 
signal essentially undistorted at all anticipated rates, 

(c) The carrier power -to -noise spectral density is 73,8 dB-Hz, 
as shown in Table I. 

The performance is still given by (30) and (31), but the equivalent 
spectral density is now given by 


'n 


eq 


No 

(0) = = 




N 

+ — Pb> ^ 


P, P„ T 
ABA 


2T 


(36) 


B 


The resulting signal-to-noise ratio in the tracking loop bandwidth is then 
given by 


(2 a - 1)' 


4N B 
^0 L 


1 + + 


N„T„ 
0 B 


^(l - o) 


(37) 


The performance of the Costas loop for the forward link employing 
the high data rate link is then given by (37) and (3 0), 

2. 4 Performance Computation - Forward Link 

The performance of the forward link Costas loop using the high data 
rate link for tracking is shown in Figure 8. The RMS phase error in degrees 
is again shown versus the power division for the unbalanced QPSK signal. 
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The set of curves designated (3) gives performance for the best 
set of rates fromthe RMS phase error point of viev/. In both the forward 
and return links, this corresponds to the highest rate possible, since this 
minimizes the spectral density of the self-noise. The set of curves (4) 
gives tracking performance for the worst set of rates. 

The variation in performance over the possible sets of rates is 
not as large for the forward link as with the return link. Throughout, we 
have imposed the restriction that R 2 since this restriction was used 

in the development of the self -noise spectral densities in Addendum A. 

F or both the forward and return links at 0 = 0 , 8, it is seen that 
the Costas loop receiver performs satisfactorily over the entire range of 
anticipated rates. We shall consider the case when the rate in the high 
data rate link goes to zero separately, 

3,5 Return Link Costas Loop Pei'formance When R^ = 0 

The performance of the Costas loop receiver is now considered 
for the return link when the rate in the high data rate link is set equal to 
zero. That is, we set 

R^ = 0 . (38) 

The powers and P^ remain present, so that the predominant 
signal for tracking purposes is a monochromatic tone in the high data rate 
link. 

The dynamic error signal when R =0 and B =100 MHz can be 

A 1 

determined directly from (23) by setting a(t) = 1, so that 

P - P 

e(t) = ^ 

+ cos 24> 

+ n (t) n (t) 

c s 

+ n^(t)[^^^^ sin 4> + JPg b(t) cos rt>] 

+ cos ^ - /Pg b(t) sin 4>] 



(39) 
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By using the same approach as that used in Section 3. 1, the auto- 
correlation function of n ^(t) is now given by 


+ R.^(T)fP^-f Pg \(r)] . 

2 

The spectral density of ’’) is given by (26). From (27), the 
spectral density of Pg Rj^(t)R^{t) is (NQ/2)Pg, The spectral density of 
P^Rj^(t) is (Nq/ 2)P^. For b(t) a Bi-c|>-L waveform, the spectral density 
of P^ Pg R^(t ) is (see e. g. [4] ): 


Pa Pn T 
A B s 


sin^ (iTf Tj^/2) 
(vf Tg/2)^ 


(40) 


(41) 


In the vicinity of f = 0, this spectral density is negligibly small, and can 
the 
by 


therefore be neglected. The spectral density of n^^(t) at f = 0 is then given 


Nr 




ea 


N ^ B 
0 1 




(42) 


and the corresponding Costas loop signal -to-noise ratio is 


P = 


4 N B 
0 L 


(2g- -U 
" -T- 


(43) 


The phase error variance of the Costas loop under the assumption of linearity 
is then related to p by (31?). The result is plotted in Figure 9 for the return 
link for loop bandwidths of 5 kHz and 10 kHz. It can be seen that perform- 
ance is tiniformly better than when 


384 ksps < < 100 Msps , 


(44) 
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which is shown in Figure 7. Since this is the case, we henceforth consider 
only the case where the range of rates is as shown in Figure 7. 

4, 0 ERROR RATE PERFORMANCE 

The resulting error rate for the forward link has been approximated. 
The tracking loop is assumed to be operating at sufficiently high signal-to- 
noise ratio so that the linearized assumption that we have made throughout 
can be maintained. The phase error is then approximately constant through- 
out the symbol time of T^ seconds and is approximately a zero mean Gaus- 
sian random variable with variance 

2 1 


O' I 


> “ 4p 

where p is given by (3 7). The symbol energy-to-noise spectral density is 


(45) 


fi 

N, 


P T 
T A 

N, 


0 0 

In the computation presented in Figure 10, we have set 
1 


R, 


= 3 Mbps 


A 


B, = 0. 01 R,, = 30 Iddz 

L A 

= 6 MHz 

Q = 0. 8 . 

The RMS phase error is then 

P 

O', = 3.8°RMS, when = 73.8dB-H2. 

The expression for the error rate used in Figure 1 0 is 

2 

G(4> ; a ; ) d(j) 


-/2 _ 

i 


1 eric 

-•rr/2 

y 

o 

a 

0 


(46) 


(47) 


(48) 


( 49 ) 


where erfc( ) is the complementary error function 





f 


I 


* 
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erfc(x) 



- 1/2 

e 

1J2 T! 


du 


and G(4i; a ; or^jj ) is the Gaussian pdf of with mean a and variance . 
The result was shifted by 5. 3 dB to the left to represent the ideal gain 
obtained from Viterbi decoding of an R = 1/3, K = 7 convplutional code. 
The result is compared to 


^°ideai 


erfc 


i 


2 E, 

D 

^0 


(50) 


( 51 ) 


which was also shifted to the left by the same amount. These approximate 
results show that only a few tenths of a dB are lost from ideal when employ- 
ing the Costas loop which is tracking the high data rate channel. This rough 
analysis does not account for (1) differential encoding, (2) interchannel 
interference, and (3) symbol synchronization. 

It is felt, however, that all of these effects would not alter the 
fact that the Costas loop considered in this study should remain as a viable 
candidate for both the forward and return link receivers. 
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ADDENDUM A 


Part I. Spectrum of Product of an NRZ Data Stream and a Bi-tp-Xi 
Data Stream . 

In this first part of Addendum A, we determine the spectrum of the 
' product of an NRZ waveform and a Bi-cp-L (Manchester) waveform. In 
the computation we impose the following restrictions; 

i) The high rate signal has rate R = 1/T and is an NRZ 
signal. 

ii) The low rate signal has rate R_ = 1/T_, is a Bi-cp-L 

is B 

signal, and 


R 


A 




or equivalently 


(A-l) 


T 


A 



(A -2) 


ill) The two sequences are independent and are operating 
asynchronous ly * 

In Figure A-l, the autocorrelation functions R (r) and R„(t) are 

shown for the case T. = 'aT^, 

A B 

The autocorrelation of the product of NRZ and Bi-Cp-L under the 
above as j ^mptions is given by 



(A-3) 


NRZ 


R 3 (t) 


Tb/2 


Ra(t) 



'T. T . 
A A 


s T3/2 


Figure A- 1 . The Autocorrelation Functions R (t) and R (t ). 

■A, S 


The power spectral density is given by 


S(tu) = ^{R(t)] 


= 21 cos (UJT)R{T)dT 
*^0 


2(3 1 / 3 1 \ 

2 )T~ ■’'It” 

«) I ^A \^B ^a/ 


A B _ 


2)3 1 / 3 1 \ 6 sin UJT 

2|T **’It "tI *^^A “ ~T T ti) 

«) 1 B A \ B ^a/ ^ 


sin OUT 


A B / 

2 

sin (ujT^/ 2 ) 
^ (U)T,/ 2 )^ 


(A- 4 ) 
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This spectral density is essentially flat over the bandwidth of 
the tracking loop. The two-sided power spectral density at iu= 0, is 

S(0) = lim S(cu) = T (l - ^ 7 ^ ) watts/Ha (A-5) 

uhO \ ^ B / 

Part II . Spectrum of the Product of Two Bi-cp-L Data iStreams , 

In Figure A-2, the autocorrelation functions of two Bi-cp-L waveforms 

are shown where the assumption is made that T ^ , and the data streams 

A i3 

are statistically independent and asynchronous. The autocorrelation of 
the product of the two Bi-cp-Li signals is then equal to the product of the 


Rb(t) 



T 



T^n 





Figure A-2. The Autocorrelation Function of Two Bi-cp-L. Signals 
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individual autocorrelation functions. Hence 


R(t) = 




) 


1 - 3 ! T 




2 
r 

. T 
A B 


(A-6) 


and 


R(t) 


Bi.iVi.iiii 


B 


■) 




' G'a'' 4 ) ■ 


(A-7) 


where we still assume R aZR^,, The spectral density of the product 
waveform is obtainea from 


= 2 


I. 


S(uj) = 7{R(t)1 

0 , ■ 'l'» ■ 'a'ii 


COS OJ T dr 




1 

)■(■•*)!= 


COS U)TdT 


(A-8) 


These integrals can be evaluated by parts, which after some algebraic 


manipulation 


X 


I 
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S(oj) = T, 


( 


sin (ujT^/4) 


(ujT^/4) 


) 


3T^ /sin^(ujT /4) 

cos(ujT^/2)| 


2T 


B 


\(U)T . /41 
A 


) 


+ - 


12 


2 

U) T. 


B 


sin(U)T /2) sin(UJT ) "] 
1 - 2 . „ ^ 


(«JT^/2) 




{A-9) 


This spectral density is shown in Figure A-3 along with the original 
Bi-tp-L (Manchester) spectmms for the forward link. In this case 
also, it can be seen that the spectral density is very flat near the 
origin. The spectral density at f = 0 is given by 

"a 

Urn S(w) = (A-10) 

uhO B 


These spectral densities are used to determine the effect of self 
noise on the performance of various candidate QPSK tracking receivers. 
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APPENDIX C 


MODULATOR POWER ALLOCATION AND CARRIER TRACKING 
ANALYSES FOR THREE -CHANNEL ORBITER 
Ku-BAND RETURN LINK 


by Marvin K. Simon 


1. 0 THREE-CHANNEL PHASE MULTIPLEXING MODULATION 

Consider the modulator block diagram illustrated in Figure C-l> 
whose purpose is to generate a signal for simultaneous transmission of 
three channels of information on the Orbiter's return Ku-band link. Two 
of the channels represent independent data channels (one having a data rate 
up to 50 Mbps while the other has a rate up to 2 Mbps) and the third channel 
consists of operational data at a rate of 192 kbps. The structure of the mod- 
ulator is such as to form an unbalanced quadriphase signal wherein the high 
rate data signal, is bi-phase modulated on the in-phase carrier 

n/F sin Wpt, and the sum of the two lower rate signals, 

^ 3 ^^^ (after being modulated onto separate squarewave subcarriers 
Sq^it) and Sq^it) ), arc added and bi-phase modulated onto the quadrature 
carrier, ^/F cos w^t. The sum of the in-phase and quadrature modulated 
carriers, i. e., the unbalanced quadriphase signal s(t) is then bandpass, 
hard limited and power amplified. 

In view of the above, the signal generated by the modulator of 
Figure C-1 is 


where 


s(t) = ^^2 [C(t) cos + S(t) sin co^t] 


C(t) = Sq 2 (t) m 2 (t) + Sq^{t)m^(t) 

S(t) = m^(t) * .j(t) 


(C-1) 


(C-2) 









1 


1 






3 


and binary waveforms. Alternately, in polar 

coordinates (amplitude and phase), equation (C-1) can be rewritten as 

s(t) = Jz V(t) sin + m) (c-3) 

where 

V(t) = Jc^{t) + S^(t) 

4>(t) = tan ^ “1^ . (C-4) 

Passing s(t) of equation (C-3) through a bandpass, hard limiter 
preserves the phase 4>(t) . Thus, the resultant first zone output is given 
by 

zj(t) = JzF sin (w^t + 4>(t)j (C-5) 

where P is the total power in the first zone -after amplification. In terms 
of in-phase and quadrature components, equation (C-5) can be rewritten as 

z^(t) = sinu^t + cos w^jt] . (C-6) 

From (C-Z) and (C-4), we have that 

V(t) = + 2 .^(t) 

= ^2 + ^3 " =2<‘> '3<'> • I'"-"' 

Let P,_ - P + P + P denote the total input power. Then, 

1 I Cl i 

v(t) = (i + (c-8) 


or, equivalently, 


I 


I 


i 


1 


t 



t*- 

^1 

n 



1 

4. 

1 


2 

[Ji+2^P^/Pt7P3/Pt ^ 



- 3,(t) S3(t) - 




(C- 9 ) 


Substituting (C-E) and (C-9) into (C- 6 ) gives 


i<‘> = " K^/^ - V 

+ ^[^ B^(t) - C^Jp^/P^ s^(t) B^lt) S 3 (t) j sin (Upt 

= =i“* • 


/P^ 8^{t) S^ii) S3(t) 


sin w^t 


+ 


'^U^z 


s^(t) + /P 3 S 3 (t) 


cos u^t 


(C-10) 


whe 


re 


Pi = (P/^T>PlS 


P 2 = <P/PT>hjP 2 

P 3 = (P/P^) |c JP 7 - 


Pd = IP/Pt>Pi''2 


(C-11) 


If we arbitrarily set the output level of the power amplifier equal 

to the total input power level, i. e. , P =Ptp » then for fixed values of Pj^/P^, 

P_/P_, and P-/P_ , one can compute from (C-11) the corresponding values 
b X ^ X 

ev . 

of for the transmitted unbalanced quad- 

riphase signal of (C-10). For example, suppose the high rate channel con- 
tains 80% of the totai power while the remaining 20 % is split with 80% going 
to the next highest rate channel and 20% to the lowest rate channel. Equiv- 
alently, we have that P^ /P^ = 0.8, ^ 2 ^^^ " ~ Then, 









J 


^ 1 ^ 


6 


from (C"ll), the transmitted power ratios are P /P = 0.8157, P /P = 

X Jm. b X 

0. 1503, P,/P.„ = 0.0287, and P^/P^ = 0. 00529* 

3 1 CL 1 

2. 0 CARRIER TRACKING LOOP ANALYSIS 


The receiver for the three-channel quadrature multiplex signal 
employs a Costas loop for carrier recovery (Figure C-2). The received 
signal plus noise, x(t) = s[t, 0(t)] + n.(t), is obtained from (C-lO) and is 
of the form 


x(t) 



■ JK ' 3 “’ 


sin #(t) 


(t) + /P^ S2(t) 


cos $(t) + n^(t) 


(C-12) 


where ^'(t) = to^t + 6(t), with the radian carrier frequency and 
0(t) = 0Q + Opt the input phase to be estimated, and n^(t) is the additive 
channel noise which can be expressed in the'form of a narrowband process 
about the actual frequency of the observed data, i.e. , 


n.(t) 


sjZ |N^(t) cos ^(t) “ #(t)| 


(C-13) 


where N^(t) and Ng(t) approximately statistically independent, stationary, 
white Gaussian noise processes with single-sided noise spectral density, 

Nq w/Hz, and single-sided bandwidth, < u^/Ztc. 

The in-phase and quadrature -phase reference signals (see Figure 
C-2) are respectively given by 


r^(t) = JY' sin ^{t) , 

r (t) = Jz" K, cos <t(t) . 
c 1 


(C-14) 


Denoting the in-phase and quadrature -phase detector gains by K^, then 
the phase detector outputs are (ignoring second harmonic terms) 


I 


I 


1 


1 


e (t) ^ K x(t) r (t) 
s ms 


= ‘l<*> -/"d ^(‘)«2(‘)=3W -N^(‘)|<=o= *(« 


- K 




sin <j>(t) 


e (t) = K x{t) r (t) 
c me 


= ^l^rrlsFl =,'«.^2(‘)=3W -^(‘) 


sin 4 >(t) 


+ K 


Sj(t)+N^(t) 


COS 4*(^) 


(C-15) 

where <^(t) = <l> - ${t) is the loop phase error. After lowpass filtering 
with in-phase and quadrature -phase arm filters, G(s), the phase detector 
output signals become, respectively, 


2 (t) = K K 
s 1 m 


- K 


Si<‘) - ^/^ =d<‘> - 


cos <j>(t) 




z^(t) 




d ®d<‘> - 


sin <^(t) 


+ K 


K \ Jp~ 

1 m I.N/ 2 


s^(t) + ./P 3 ^ 3 (t) + N^(t)| cos 


(C- 16 ) 


where Sj^(t) = G(p) Sj^(t) ; i = 1,2,3 
yt) ^ G(p)[s^(t)s2(t)s3(t)] 
N^(t) = G(p) N^(t) 


A (t) = G(p) N^(t) 


{C-17) 


Multiplying these two lowpass filter outputs (assuming, for simplicity, that 
this multiplier has unit gain) gives the dynamic error signal 


I 


i 


I 
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z^{t) ^ zjt) zjt) 


K ^ 
1 m 


Kf K ^ 
1 m 


K ^ 
1 m 

2 


^^/^^(‘> -'7^ ®2<‘' ®3<*> 

l^/^ ^<*> - 7^ «d'‘'|17^ ^2'« ®3'‘> 

- - z[j¥; s,(t) -7?7 V‘>l ^s<‘> 

- ==I7^V«+,^V‘'1 V> 


sin 2^(t) 
cos 2 0(t) 


sin 2 0(t) 


K ^ K ^ 

+ M I?' §,W - ,/?7 



1 


AJh ®2<‘> - ^3<‘>1 ^s“’ - ™s<*> ^c<‘> 


COS 20(t) 


(C-18) 


The instantaneous frequency of the VCO output is related to z^Ct) by 

^ = K^[F(p) z„(t)] 1. 

and hence the stochastic integro -differential equation of loop operation 
becomes 


(C-19) 


2d0{t) 

dt 


= 2 n* - KF(p) 


J^Z «l(‘'®2<‘> ®.<‘>®3<‘> -7^ 


7P^ 


cos 2 0(t) + V [t, 2 0(t)] 


(C-20) 


'y O 

where K ^ K, K K„ and 
1 m V 


M 


1 


t 








f 


v,[t.2p(t)] = 


-N^(t) + N^(t) - 2 
c s 


^^(t) + ^3(t) 


+ 2 


=i<‘> - JK “d<*> 

4 - 1^2 =2<‘>^•^/^3 


=1<‘> - ^d'‘> 

sin 2 0(t) 


N (t) 
s 


N (t) 
c 


N (t) 
c 


N (t) 
s 


2N (t) N (t) 
s c 


cos 2 0(t) 


(C-21) 


As is customary in problems of this type, we now decompose each of the 
signal terms in (C-EO) into its mean value plus the variation about this 
mean, e. g. , 

(t) sin 2 0(t) = s^(t) sin 2 0(t) + ['s'^(t) - s^(t)] sin 2 0(t) (C-2E) 


where the overbar denotes statistical expectation. The bracketed term 

2 

in (C-22) represents a self -noise due to the modulation, s^^ (t), and such 
terms would result for each signal component in (C-20) when decomposed 
as in (C-22). Since all three data rates, = 1/T„. (i = 1,2,3), are large 

Si 

with respect to the single -sided loop bandwidth B , then as has been shown 
previously for the case of a Costas loop tracking a bi-phase modulated 
carrier, the self -noise terms can be neglected. Furthermore, 

a.^(t) ^ D. = f'”s.(f)lG(j2irf)pdf ; i = l,E, 3 

^ *^C0 


8^(t) = = r"s^(f)|G(j2TTf)l^ df (C-23) 

^CO 

1 2 

where S^(f) is the power spectral density of s^(t); i = 1, 2, 3 , |G(j2'iTf)| is 

the square of the magnitude of the arm filter transfer function, and 

Sj(f) = S,(f) * S_(f) * S-(f), where the asterisk denotes convolution. In 
d 1 2 3 






10 


addition, since m^(t) is independent of m^(t); i / j = 1,2,3, then 


S,(t) & (t) - 0 i.j = 1,2,3; i j 

J 


l(t) ^ ,(t) = 0 i = 1, 2, 3 


(C-24) 


Making use of (C-23) and (C-24) in (C-20), along with the assumption 
regarding the neglect of the modulation self-noise terms, we get the 
simplified stochastic equation of loop operation 


d<a(t) 

dt 


2flg -KF(p) -P 3 D 2 +^d^d^ sin2,*(t) 


+ [t, 2 0 (t)]| . 


(C-2S) 


We now proceed to evaluate the noise spectral density of the equivalent 
delta-correlated noise process, V-[t, 20(t)], i. e, , 


/ 


CO 


N = 2 / R„ (t) dr 

sq ' f 


(C-26) 


where 


Ry (r) = V [t, 20(t)] V [t + T, 20{t)] . 


(C-27) 


Substituting (C- 21 ) into (C-27), we obtain, after considerable algebraic 
manipulation, 

(t) = 4j[P R. (t)+P R^ (t)+P R^. (t) + P R;, (t)R (-)R^ (t)]Ra(t) 
2 i Si ^ as, N 




R^(t) 


(C-28) 


where 


R^ (t) = §.(t) ^.(t + T) : i = 1 . 2,3 

S, 11 

1 




(C- 29 ) 
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Integrating (C-28) between -co and co and noting that 

rN. 


/'“R^(T)dr = r[^lG(j2.ff 


i 2 


df 


.00 


f R^(t)Ra(t) dr 
/ s . iN 

fc'-CO 1 


/*”ra (t)Ra (t)Ra (t)Ra(t) dx 

J-w ®1 ®2 


N 

~Z 


- ^”s.(f)lG(j 2 Tr£)P df 

i = 1,2,3 


N 


° r"s^(f)|G(j2Tr£)l^ df (C-30) 
•'-00 


with 


S^(f) = S^{f)lG(j2TTf)l^ 82 (f) lG(i2TTf)l'^ ❖ S^(f)|G(j2Trf)!‘^ (C-31) 


we get an expression for given by 


N = 4N^ 
sq 0 


3 

ZP. y*“s.(f)|G(j2Trf)l^df + ^“s^(f)|G(j2Trf)|‘ 

2 /“[“lG(j 2 irf)l^] df 


df 


S.(f)lG(j 2 irf)l^df + ^*S^(f)|G(j 2 Trf)|‘ 


df 


(C-32) 


Finally, from (C-26), the equivalent loop signal-to -noise ratio p' is 


, 'Pl°l - " Vd> 6 _P , 

’’ N B 4'^L 

sq L 


(C-33) 


where p = P/N_B is the loop signal-to -noise ratio of a phase -locked 
0 L> 

(CW) loop operating on the total power P and 




0 


L N 


sq 




4 N„ fP,D, - P^D„ - P^D^ + P ,D ,1 


1 1 2 2 3 3 d d 


(C-34) 
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is the loop squaring loss. What remains is to express in terms of the 

basic system design parameters. Defining the modulation indices q . ; 

i = l,2,3, and q, by 

^ p. 


- ^ 

M P 


we obtain, after much simplification, 


(C-35) 



(C-36) 


where 


D = n, D + 1 == average filter distortion factor 

11 2 2 33 dd 


K, 


A •'-00 


r"|G(j2Trf)f 


df 


= filter constant 


y'*|G(j2irf)l‘ 


df 


B, = r |G(j2irf)|^ df = two-sided arm filter noise bandwidth 

^ “'-CO 


Ehi /■“s.{f)|G(j2i,f)l^df +ri^ /*“sA(f)|G(j2Trf)l^df 

^7” i=l •'-(a -m 

~ 3 CO 03 

E 1 j /“s.(f)|C(j2.Tf)|^ df + /■"s_|(f)lG(j2rf)|^ df 

i=l *^00 *'-<0 

= combined filter - signal spectrum constant (C-37) 

and = 2P/NqB. is the total power -to -noise ratio in the arm filter band- 
width. Alternately, since the bandwidth will be optimally chosen relative 




1 


t 


to the high data rate symbol stream at .’ate 01 , then 


- "zOa - ’l3°3 ’'d“d> 


4 = 


fiA 

k 2 Rd, / D 




where = Pj^T^/N^ is the data signal -to -noise ratio for the high data 

rate channel, r* • 

The simplest case, at least mathematically, to use as an illustra- 
tion is that of an ideal arm filter for which 


|G(j2irf)l' 


1 ; |fl SB./2 


0 |f| >B./2 

' ' 1 


In this situation, we see that K = K = 1 . Thus, (C-36) simplifies to 

JD 1j 

^ - '^2°2 - ’^3°3 ^ ’^d°d>' 

^ ’'d^d’ k 

To simplify (C-40) any further requires specification of the format of the 
three signal components. As an example, suppose that s^^(t) = m^(t) is 
an NRZ modulation of rate = l/T. . Then, 


/sin irf Tg \2 


For an ideal filter, 
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which can be expressed in terms of the tabulated sine integral Sj{x) function. 

Furthermore, suppose S 2 (t) = m 2 (t)Sq^(t) is an NRZ modulation of rate 

~ bi -phase modulated on a square wave subcarrier of 

frequency f .Then, 
s c 


where 



(C-43) 


(C-44) 


Finally, if ^^(t) = m^(t)Sq^(t) is a Manchester coded modulation of rate 

^Sg “ 1/^83 * bi-phase modulated on the quadrature square wave sub- 

carrier Sq_(t) = Sq_(t - l/4f ), then 
3 c sc 


s,(f) = Z J 

i , mr m„ sc 

n=l,3,5,... i 


(C-45) 


where 


S 

m 


(f) 

3 



(C-46) 


For an ideal filter, and can be evaluated by direct numerical inte- 
gration of (C-43) and (C-46), respectively. 


3. 0 PROBABILITY DENSITY FUNCTION OF THE LOOP PHASE ERROR 

Using Fokker -Planck techniques, the steady-state probability den- 
sity function (p.d.f.) p( 2 <^) of the modulo 2tt reduced phase error 24> can be 


1 
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determined from (C-25). Assuming a loop filter of the form 

i + ST_ 

= rTT^ > 

1 + ST, 1 2 1 


{C-47) 


then (Ref. 1, Ch. 2, and Ref. 2, Chs. 9-10), 

^2<p + Zv 

p(2^) = C' exp (p2<t> + a cos Z<t>) f exp (~p2x - a cos 2x)dx ; 

° 2<*> 


where C’ is a normalization constant 
0 ■ 




f''p{Z^)d4> = 2 r p(2(/.) 

-Tt/2 


d.^ = 1 


and 


« = (^)p' - 




2 


rtr. 


|(^| < Ti (C-48) 


(C-49) 


(C-50) 


2 

p = (^) sir 


+ a sin 2<p 

where r = second order loop damping parameter 
loop damping 


P' 


effective signal-to-noise ratio in the loop bandwidth 


(C-51) 


<T - sin 2(l> - (sin 20) 
Cl 


(C-52) 


t 


In. the linear region of operation, the variance of the loop phase 
error (tracking jitter) is given by 


where p' is found from (C-33) and (C-36); whereas, in the nonlinear 
2 

region, tr must be determined by numerical integration, viz. , 


tr =2 
2 ^ 


/ 

•'-•rr/; 


(20 - 20) p(20) d0 


4. 0 CROSSTALK CONSIDERATIONS 

The presence of tracking phase error in the receiver introduces 
crosstalk from one data channel into another. In particular, we observe 
from (C-15) that the low data rate channels introduce crosstalk into the 
high rate channel in proportion to their relative received powers and the 
carrier tracking loop phase error. Defining D^(il)j^ e<lual to the ratio of 
power in the kth channel to the crosstalk power due to the data in the jfth 
channel conditioned on a fixed carrier tracking loop phase error 0, we see 
from (C-15) that 







cos 0 

^1 

p 

2 

. 2 

sin 0 
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1 

2^ 
cos 0 

P 
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^3 
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cot^ 0 


2 
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where the ratios P^/P^ and P^/P^ may be determined from (C-11). The 
average crosstalks D^(2) and Dj^(3) may then be evaluated by averaging 
(C-55) over the p.d.f. of (C-48), i. e. , 
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(C-56) 


The compxftations required in (C-56) must be carried out by numerical 
integration. 
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APPENDIX D 

EFFECTS OF PHASE NOISE ON THE PERFORMANCE OF 
TWO-WAY COHERENT COMMUNICATION LINKS 

Gaylord K. Huth 

This appendix systematically considers the effects of transmitter 
phase noise on the performance of two-way coherent communications links 
and indicates design approaches that may be used to reduce or minimize 
such effects. For typical transmitter phase noise spectra, the resulting 
degradations in bit error rate performance for baseband and phase -shift- 
keyed subcarrier channels are determined. 

1.0 INTRODUCTION 

Figure 1 illustrates functionally the generalized two-way coherent 
communications link which is considered in this appendix. Each portion 
(forward and return) of this link is assumed to consist of a baseband digital 
data channel, sinnmed together with n phase -shift-keyed (PSK) digital data 
subcarrier channels prior to phase modulation of the transmitted carrier. 

Of concern here are the effects of transmitted phase noise on the bit error 
rate performance of these various channels. For the forward portion of the 
two-way link, the transmitted phase noise is due primarily to jitter intro- 
duced by the timing source which generates the carrier frequency. For the 
return portion of the link, the transmitted phase noise is due primarily to 
the turnaround of the phase noise which is not tracked out by the forward 
link receiver carrier reference loop. 

2. 0 TRANSMITTER PHASE NOISE SPECT RA 

As a first approximation, the effects of transmitter phase noise on 
bit error rate perforjnance of communications links are sometimes deter- 
mined assuming that the phase noise power spectrum is white, as in Figure 
2(a), In practice, the transmitted phase noise spectral characteristics are 














(rad /Hz) 



10 1000 
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Figure 2. Transmitter Phase Noise Spectra 
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very dependent on the configuration and characteristics of the circuits used 
to generate the transmit carrier frequency. Ht-nce, trcinsmitted phase noise 
spectra may vary considerably from one system to another. Figure 2(b) 
illustrates the phase noise spectrum typical of what might be encountered 
over a number of systems. The approach presented subsequently for eval- 
uating the overall effects of transmitted phase noise is general and applies 
for any spectral characteristics. Only the final numerical results are 
dependent on the assumed spectral characteristics. 

As a specific example of phase noise spectrumst consider the fre- 
quei..cy synthesizer necessary for the S-band payload interrogator in the DSN 
mode presented in Figure 3. Phase noise is generated from each oscillator 
(i. e. , the TCXO, the VCXO, and the VCO) and ' om each phase detector. 
These sources of phase noise are accentuated by the multipliers and attenu- 
ated by the loop filters. First consider the phase noise generated by the 
TCXO, the VCXO, and the associated phase detector as modified by the 
multiplier and loop filter. Figure 4 analyzes each of these sources and 
combines them to form the total phase noise at the output of the VCXO loop. 
The asymptotes for the loop filter are also shown. Note that the loop filter 
attenuates the low frequency response of the VCXO due to the feedback and 
the high frequency response of the TCXO and phase detector by direct 
filtering. By adding up the resultant curves of the filtered phase noise 
source shown in Figure 4, the total output phase noise spectrum is obtained. 
Note that the loop filter is actually 3 dB down at the 100 Hz cutoff frequency 
rather than the asymptotes shown and the total phase noise spectrum curve 
takes this into account. 

The next source of phase noise in the system of interest is the output 
of the VCXO multiplied by 36. If this multiplier is implemented with a 
phase locked oscillator (PLO), then the effect at the high end of the spec- 
trum can be improved. Figure 5 analyzes the resulting phase noise spec- 
trum using the X36 PLO. First, by multiplying the output VCXO phase 
noise spectrum obtained in Figure 4 by 36 will give the resulting phase noise 


4 ~ 

t 


445. 0694444 
to 



Figure 3, Frequency Synthesizer for S-Band Payload Literrogator in the DSN Mode 


U1 



dBc/Hz (single ' sided) 



Figure 4 


Estimated Phase Noise Density for VCXO Loop 
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spectrum at the low frequency. The high frequency spectrum will be attenu- 
ated by the loop filter. To this filtered VCXO spectrum must be added the 
effect of the VCO and the phase detected of the PLO as filtered by the loop 
filter similar to the procedure used in Figure 4,. 

Using Figures 4 and 5, the total phase noise spectrum for the fre- 
quency synthesizer is obtained in Figure 6. First, the VCXO output phase 
noise from Figure 4 is multiplied by 5 and added in at the data modulation 
input. Second, the VCXO multiplied by the X36 PLO from Figure 5 is added. 
Finally, the output of the VCXO is passed through the variable multiplying 
PLO in the upper channel of Figure 3. The worst case for the low frequency 
spectrum in this channel is when the multiplication of the VCXO is 1361/135. 
Again the loop filter attenuates the high frequency of the multiplied VCXO 
spectrum. The other contributions of the phase noise are the loop filtered 
phase detector and the VCO, with its low frequency attenuated by the loop 
filter. By adding all these contributions, the total phase noise spectrum 
at the output of the frequency synthesizer is shown in Figure 6. 

In the remainder of the appendix, the phase noise is referred to as 

2 2 

a random process with zero mean and a variance , To determine (T^ 
from the total phase noise spectrum in Figure 6, note that, for a zero mean 
random process, 

= f ( 1 ) 

where Sj^(u) is the power spectral density of the phase noise process. It 

remains to be shown that the oscillator output power below the carrier 

(i. e. , the center frequency) in Figure 6 is actually the power spectral 

density of the phase noise process. 

A waveform of frequency f and constant amplitude A can be 

c 

expressed as 

V = A sin[27rf t + 0(t)] (2) 

c 


t 
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where 0{t) represents the phase variations about the linearly progressing 
phase* The instantaneous frequency is defined by 


(3) 


Often the phase or frequency modulation of the carrier can be repre- 
sented as a sinusoidally varying quantity 


6(t) = ^<p sin Zirf t 


m 


(4) 


where /i<p is the peak phase deviation in radians* and (2) can be written as 
V = A sin (2irf^ + A^sin 2irf^t) , 


(5) 


Using a Bessel function series expansion* (5) can be written to empha- 
size the frequency components of a simple phase modulated signal. With 
Cl) “ 27rf, 

V 

T = sin Cl) t + J (A<^)[sin (ci) +0) )t - sin (o - w )tl 

Au Cl cm cm* 


+ J (AC;P)[sin (w + 2Ci) )t + sin (w - 2ci) )t] + . . . , 


m 


m 


(6) 


Now assuming that A(p «1, which is reasonable for a good oscillator* 

the Bessel fvinction Jq(A<P) = 1, , and J^* J^* etc.* are negligible. 

Thus* the waveform appears in the frequency domain as a carrier at f with 

' c 

sidebands spaced at f +f and f -f The ratio of one sideband amplitude 

c m ^ c m 

to the carrier amplitude is « Therefore* if the power of spectrum below 


the carrier is known, then 


2 : 


sb 


(7) 


where P . is the power in the sideband f +f and P is the power in the 
sb c m c ^ 

carrier. If more than one sideband exists* then the total phase deviation 
would be the sum of the square root of the ratios of power in each sideband 
divided by the carrier power. Therefore* the spectrum estimated in 


.. . 
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Figure 6 is the power spectral density of the phase noise process and the 
variance is determined by equation (1). 

It will be shown in the following section that the variance of the phase 
noise process outside the loop bauidwidth of transponder affects the perform- 
ance of the forward link and the phase noise inside the loop bandwidth affects 
the performance of the return link. The loop bandwidth (double sided) of 
the DSN transponder is 18 Hz. Therefore, for the forward link, the phase 
noise spectrum can be integrated from 9 Hz to, say, 1 MHz and then doubled 
for the two-sided bandwidth to determine the phase varicince. For the phase 
noise power spectral density in Figure 6, the rms phase noise is 2.7° 
for a double sided bandwidth from +9 Hz to +1 MHz. 

Two additional points must be made about the phase noise spectrum 
in Figure 6. First, with the final loop bandwidth equal to 5 kHz in the upper 
channel of Figure 3, spur suppression is difficult with a digital phase detec- 
tor even though it is fairly easy with a sample -and -ho Id detector. To alle- 
viate this difficulty, the final loop filter could be decreased in bandwidth 
to 2 kHz. This would make spur suppression easy with either type of phase 
detector but would increase the rms phase noise. Figure 7 presents the 
phase noise power spectral density with th j loop bandwidth at 2 kHz. The 
resulting rms phase noise is 3.4*^ for a double sided bsoidwidth from 
jf 9 Hz to ^1 MHz. Finally, to simplify the frequency synthesizer, a X36 
miiltiplier could be used instead of the X36 PLO. Figure 8 presents the 
phase noise spectral density for the X36 multiplier with no PLO. The rms 
phase noise for this simplified frequency synthesizer is = 5. 9° for a 
double sided bandwidth from jl9 Hz to _+ 1 MHz. 

3. 0 PHASE ERROR MODEL, 

As shown in Figure 9, the received signal Z(t) is 


Z(t) = s(t) + n(t) 


( 8 ) 


dBc/Hz 










s(t) = Jz Asia («^t + 9(t) + ^^p(t)) 

.(t) = >lz[n cos 0 ) t - n sin u) t] 
i c c s c 


r(t) = JZK^ cos I 


Figure 9. Carrier Tracking Loop 
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where s(t) = A ain(,w^t + 9(t) + ^^p(t)) 
= angular carrier frequency 
0 (t) =s modulating process 
'5',pp(t) = transmitted phase jitter 


A = total received i>ower = P 


rec 


and n(t) » channel white Gaussian noise with flat spectral density 

Nq (one-sided) over the frequency reinge of the receiver. 

The bandpass filter with bandwidth centered at carrier frequency 
is used to reject as much noise as possible to the phase -locked loop (PLIj) 
without distorting the signal. The output of the bandpass filter can then be 
given by 


Zj(t) = s(t) + n^(t) 


( 9 ) 


where the filtered noise n.(t) can be expressed in the form of a narrowband 
process [l] about the nominal center frequency of the bandpass filter, i. e. , 

n.(t) = fz [n (t) cos 0) t - n (t) sin u t] , (10) 

1 V c c s c 

Both n (t) and n (t) are statistically independent Gaussian noise processes 
c s 

with identical symmetrical spectra corresponding to the lowpass equivalent 
of n(t), i. e. , 

B. 


"c't> ' ”s<‘» = ”o(^) = I ■ 


( 11 ) 


Assuming that frequency lock has been achievedf the expression for the 
multiplier output can be written as 

X(t) = K Z.(t)r(t) 
m i 


( 12 ) 


where K 

m 

r(t) 


«(t) 


multiplier gain 

output of the voltage controlled oscillator (VCO) 

JT Kj cos ^fiiJ^t +$(t)) 

rms voltage of the VCO output 
FLL estimate of '9^,pp(t) , 


and 


^ i I I < i • * 
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Since the PXjL cannot respond to the double -^frequency terms« we have 
X(t) = AK^ Bin ( 0(t) + «^p(t) - ^(t) ) + Kj K^n^(t) cos ^^(t) 

+ K. K n (t) sm$(t) , (14) 

1 m s 

Let "^(t) - 'St^p(t) - ^(t) = phase error; we obtain 

X(t) = K, K [ (a cos e(t) - N (t)) sin*(t) 

1 m ' S' 

+ ( A sin e(t) + N^(t) )cos \p(t)] (15) 

where 

Njt) = n^(t) cos ^^p(t) +*ig(t) sin ^.j,p(t) (16) 

and N^(t) = n^(t) cos ^,j,p(t) - n^(t) sin ’®^,j,p(t) (17) 

The modulating process can be generalized as 

m 

e(t) = p^d(t) + 53 Pgi sl» ^ 

i=l 

baseband data channel modulation index 
+ 1 K baseband digital data, bi-«^-L format 
data rate of d(t) 

subcarrier channel i modulation index 
channel i angular subcarrier frequency 
subcarrier channel i modulating process . 

By substituting (18) into (15) and noticing that the loop is narrow* 
band (as compared to R^) and cannot respond to high frequency terms. 

We obtain 

X(t) = KjK^[A‘ sin nt) 4l(t,^(t)) +N(t,^(t))] 09) 


where 3, = 
d 

d(t) = 


^ei 


(a) 


SX 


and 0gj(t) s 


t 


I 


I 


i 


4 


where 


A* = 


Jq(x) 


and N (t, ^(t)| 


m 

A cos n Jo<Psi> 

1=1 

Bessel ftxnction of first kind of zero order 
m 

Ad(t) sin p, n J-(p .) cos '^(t) 
a u SI 


N (t) cos - N (t) sin ’9?(t) 
c s 
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( 20 ) 


( 21 ) 

(22) 


From (19)» it can be seen that the first term is the loop phase error 
control signal, the second term is phase distortion due to modulating signal 
d(t), and the third term is the noise. Under certain physically justifiable 
assumptions, the noise N(t,'t(t)j can be shown [2] to.be approximately a 
zero mean white Gaussian noise process with the same spectral density as 
that of the original additive process n.(t). 

The stochastic integro-differential equation that relates the VCO 
phase output to X(t) is given by 

K F(P) 

^(t) = X(t) + ^,(t) (23) 

p 6 


where 


P = ~ = Heaviside operator 

F(P) = loop filter transfer function 
K.y = gain of the VCO (radians /volt-sec) 


’®'-(t) = VCO self -generated phase jitter. 

U 


^d 

Since '®'(t) = '®,pp(t) - ^(t), we rewrite equation (23) to 


’l^(t) = «'^p(t) 


F(P) 


X(t) - ^^{t) 


or 


m) = ^^p(t) - 


A'K 




I 


l(t,’«'(t)) N(t,'ir(t)) 

sin 't(t) + — + ~ 


»(t))l 

— J -*, 


(t) 


(24) 


( 25 ) 
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where K - K K K__ = open loop gain. 

1 m V 

A mathematic equivalent phase error model . can be conveniently 
represented by the block diagram of Figure 10. 

For practical purposes, the carrier tracking loop is usually designed . 
with sufficient care so that the tracking phase error ^'(t) will be small 
( < 7r/6). If this is the case, then the linear model can be applied { 1. e. , 
sin 'i'(t) = ^^(t); cos '|»(t) ^ 1; N (t) >>N (t)'®'(t)), and (25) can be reduced 

C 8 

to 

*(t) ? j*(t) + - *j(t) (26) 

where 

m 

I'(t) = Ad(t) sin J] . (27) 

^ i=l 


A slight rearrangement of equation (26) relates the phase error to the input 
processes: 


^(t) == [1 -H(P)3^^p(t) - [1 -H(P)]^2<^) - H(P)^ 



(28) 


where 


H(P) 


A»KF(P) 

P + A'KF(P) 


(29) 


is the closed-loop transfer function and the two-sided noise loop bandwidth 
is defined by [3], 


•'-CO 


(30) 


From (28), we can see that the total phase error is contributed by 
two in-band processes, I'(t) and N (t), and two out-band processes, (t) 
and 

Generally speaking, the data rate of d(t), R^, is designed well above 
the loop bandwidth Therefore, since the data is bi-4>-L< with no d. c. 

component, the contribution to the phase error by d(t) can be neglected. 

It is also well known that the probability density function of ^ due to 
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noise N (t) for a second -order loop can be expressed as [4]; 
c 


e® cosji» 


PNc‘*’ ' 2.rlg(a) ' 


(31) 


where 


a ~ 


N. 


(32) 


ana 


- variance of phase error due to N (t) 
c c 

l>(3c) - modified Bessel function of zero order. 
0 


From (28), we obtain 

< = -t r 


s (0>) 

2 ^c 


do) 


(33) 


where (to) is the power spectral density of N^(t) and is N^/2 in the loop 
c 

bandwidth B_ . Thus, 

1j 


N, 


2A' 


2 


. ...0 ■ . h 


A* 


m 


Since A* = A cos f] 3Ap .) and A = P we have 

d ,1 », 0 si ’*ec 

1=1 


and 


N, 


0 L 


« m , 

2 « n - 2 


P cos /3, n jN/ 3 .) 
rec *^d 0 si 


, m 

P cos^fi.UsUfi.) 
rec a 1 I 0 ^si 


a = 


i=l 


'N, 


0 L. 


(34) 


(35) 


(36) 
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which is the signal-to->noise ratio (SNR. ) in the loop bandwidth. It should 

Xj 

b< noted here that the numerator is the total power allocated to the carrier 
component. 

For large SNR^^, say greater than 12 dB, we can simplify (31) to 



(37) 


a Gaussian distribution with mean zero and variance . 

Since (t) and ^ (t) can be assumed to be zero-mean Gaussian 
processes, we conclude that the phase error 'i'(t) is also a zero -mean 
Gaussian process if SNRj^ a 12 dB and R^ » B^^. The variance of (t) 
under these conditions can be computed as 


where 


'if 


r ff 


TP 


+ (T- , + ff, 

N i 




TP 


1 

2ir 


P |l-H(jw)|^S (w) dw 

TP 


(38) 


(39) 



_L 

2it 


p |l -H(jw)|' 
•4 CO 


S^ (w) dw . 


(40) 


S-. (w) and S* (oj) are power spectral densities of ^'__(t) and 'i'_(t), 

*TP 2 TP 2 

respectively. 

Now, if we rewrite equation (28) in terms of ^(t), we have 


4r(t) = 


[1 -H(P)]4-^p(t) - [1 -H(P)]4'2(t) - H(P)(^) - H(P)(-^) 
$(t) = H(P)«'^p(t) + [l-H(P)]'*r^(t) +H(P)(^) + H(P)(--^^) . 


(41) 


or 


( 42 ) 
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Notice that ^(t) is the turnaround phase jitter to be imposed on the 
return link. From (42), we can see that ^(t) is due to contributions by the 
"in-band" transmitted jitter, thermal noise and modulation interference, 
and out-band contribution by the VCO phase jitter. The narrower the loop 
bandwidth the higher the phase error, but the less turnaround phase 
jitter and phase error contributed by the thermal noise, and vice versa. 
Therefore, a tradeoff among the loop bandwidth, incoming transmitted 
phase jitter, and thermal noise can be made to equalize the performance 
of both links to the limit of the validity of the linear model in equation (26). 


4. 0 DATA DETECTION AND BIT-ERROR PROBABILITY 


If we pass the output of the PLL multiplier through a divider with 
the gain ^ lowpass filter with the bandwidth just to pass d(t) 

as shown in Figure 11, the output of the lowpass filter can then be expressed 
as 


d'(t) 


A sin 



cos 4f(t) d(t) 


+ A cos (3. 


"m T 

J\ JQ(^g.) ain 'i'(t) + N'(t,4^(t)) (43) 

,i=l ® J 


where N'(t,'»&(t)J is the filtered version of N(t,^(t)j as given by (22). 

The first term is the desired signal, the second term is the interference 
due to loop error, and the third term is the noise. We can approximate 
sin tf(t) “ ’9f(t), Since ’^'(t) is also a lowpass Gaussian process, we can 
rewrite (43) as 


m 

d‘(t) = A sin f[ J(,(pg.) cos ^(t) d(t) + n*(t) (44) 

issl 


where n'(t) is a zero-mean Gaussian process with the power (N^) in the 
data rate (R.) bandwidth to be given by 

Cl 


N = N.R, 
p 0 d 


+ .i CO.* • 




(45) 

















i- .j 


t ■ * 
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We assume that all the power of 'i'(t) is in the data rate bandwidth. 

For Rj B_ y the most practical designs* we can also assume SE'(t) 
d 1 j 

to be constant during the data period (1/R^). In this case the conditional 
bit-error probability by detecting d(t) with a matched filter (integrate-and- 
dump circuit) can be expressed as [5]: 


P ” erfc 

0 ^ 


fP’ cos^^ 
rec 


N 


where erfc (x) = complementary error function = 


P» ^ sin^ p, fl j;^(p .) . 
rec /^d .11 0 ^sx 

1-1 

Substituting (45) into (46), we have 


m 


•'x 




dt and 


(46) 


P (’Sr) = r erfc 
e c 


.2 .2 n T^/n \ 2 

A sm U Jo(Pgi) cos ^ 

ifl 

'JO 0 m 

i=l 


(47) 


or 


P {9f) - — erfc 
e c 


E, cos -Jf 

.^d 

Nq + N' 


(48) 


where 


m 




E, 


i=J_ 
R 


(49) 


and 


= Energy per bit of baseband channel 


% C08^ n Jfl'f'si’ 

i=l 


m 


N» - 


R 


(50) 


= Equivalent noise spectral density due to phase error 


’5 


Since the probability density function of ^^{t) is Gaussian distributed 

.2 




exp 


r _iL‘ 

L 


the bit-error probability can now be computed by 


(51) 


= d^ 




erfc 



d cos '3' 


dil' 


The DSN transponder design requirement [6] 


(52) 


where N = Nq + N' = equivalent total noise spectral density. 

Figure 12 presents a set of bit-error probability curves as a func- 
tion of energy per bit/noise spectral density (Ej^/N) for various values of 
(T^ in degrees. If power spectral densities and S^ 2 (*^) known, 

cr2 can be determined from (37) and (38). 

— r- 

Using the phase noise power spectral density in Figure 6, the rms 

phase error ir_ = 2 , 7 ^, 

^TP o 


is for a maximum of 1 static phase error and = 3 . Therefore, the 

o 2 

total phase error will be less than 6. 7 and the degradation due to phase 

-4 

noise from Figure 12 at P s 10 is less than 0. 1 dB. Considering the 

6 

return link, the turnaround phase noise is described by (42). For a 18 Hz 
loop, it may be assumed that a major portion of the phase jitter is con- 
tributed by At the interrogator receiver, a similar frequency 

2 

synthesizer is used as in the transmitter. Therefore, at the receiver, 
there is a received phase noise with or^ = 3° plus - 2. 7° from the 
receiver frequency synthesizer and <r^ =■ 3° from the receiver phase lock 
loop. With a total phase error o*^ = 8,7°, Figure 12 shows that the degra- 
dation due to phase noise for the two-way link is less than 0. 2 dB for 
-4 

P « 10 . 

e 



Figure 12. Bit-Error Pro 
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As sho'wn in Figure 11, the output of the ith subcarrier bandpass 
filter can be given by 


S' .(t) = A cos 
81 






« cos ^(t) + P(t) + N'(t, ^/(t)| 
where P(t) etnd N'(t, ^(t)) are the filtered versions of 

>in d(t)(z sin(a)gjfc+ e^^Ct)) sin ^a/Ct) 


A sin 




(63) 


and N |t, '5'(t)J , respectively. 

The first term is the desired signal, the second term is the inter- 
ference due to sin^(t), and the third term is the noise. Since sin ^(t) is 
is a lowpass process, the spectrum of sin^(t) spilled over the bandpass 
filter bandwidth is extremely low. Thus, for all practical purposes, the 
second term in equation (53) can be neglected. In this appendix, we only 
consider the cases where the subcarriers are phase -shift -keyed by the 
binary data stream, i. e. , 


0„.(t) ^ d^.(t) = + 1 . 

SI SI 


Equation (53) can now be rewritten as 


S' .(t) = A cos 
si 


pd(2^i<p.t>)[n Jo(Psj^j “si^ • 

jM 

(54) 


A subcarrier tracking loop, usually a squaring loop or Costas 
loop [5], is employed in series with a bit synchronizer (integrate -and-dump 
circuit) for choherent detection of d ,(t). Generally, the subcarrier tracking 
loop is designed to operate in the high signal -to -noise ratio region and 
narrow loop filter range (as compared to carrier tracking loop). The phase 
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error introduced by the subcarrier loop can usually be neglected. Thus> 
the bit-error probability can be obtained the Scune way as that of the base- 
band channel: 


erfc 


where 


'b . 

SI 


■ C 


iii 


b . 
81 

N. 


cos 


d^ • (55) 


a2 2 

A cos 


R . 

SI 


- energy per bit of channel i 


and R . = data rate of channel i . 
si 


Figure 12 can also be utilized to determine the bit-error probability 
except that replaces N . 
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CONCLUSION 


The effect of transmitter phase noise on data transmission for a two- 
way link has been analyzed. It was shown that there is a tradeoff between 
the effect of phase noise on the forward link and the phase noise turnaround 
for the return link. As a final result^ the performance degradation to the 
bit error probability was presented as a function of the total rms phase error 
in degrees. As a typical example of a two-way link, the S-band payload 

interrogator in the DSN mode communicating with the DSN standard tran- 

-4 

sponder was analyzed. It was found that at bit error probability of 10 , 

the performance degradation was less than 0. 1 dB and 0.2 dB for the for- 
ward link and the total two-way link, respectively. 
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APPENDIX E 

COMMAND AUTHENTICATION 
Gaylord K, Huth 

For Air Force misaiona, modifications are proposed to the Shuttle 
command system^ to provide command authentication. These modifications 
are shovm in Figures 1, 2 , and 3« In all figures, there is an authenticator 
variable that can be up to 50 bits long that is known at both transmitter and 
receiver. The general approach taken was to modify the parity portion of 
the command codeword using the authenticator variable. 

In Figure 1, the variable is modulo-2 added to the command informa- 
tion word i_ generated by the command source. This combination word is 
then BCH encoded to produce a parity word £. The codeword x = (i_, £) is 
then transmitted. Note that the code so produced by a fixed authenticator 
variable is not cyclic. In the receiver, the error detector encodes the siun 
of the variable and the received information word to produce a test parity 
p*' . This test parity is then compared to the received parity . This is 
done by simply modulo-E adding and pi” . If the sum is zero, the code- 
word is accepted; otherwise, it is rejected. 

In Figure 2, the information word and the variable are encoded by 
completely different (127, 50) BCH encoders. The two parities produced 
are modulo-2 added to form p. The codeword x = (i_,£) is transmitted. The 
same process is performed on the received word in the error detector; the 
resultant parity is compared with the received parity p’ as in the first 
method. 

In Figure 3, the information word is encoded and the resulting parity 
word is permuted in some fashion depending on the variable. The permuted 
parity word is then transmitted together with the information word. Again, 

^G. K. Huth and B. H. Batson, ”A Command Encoding Scheme for 
A Multiplexed Space Communication Link,” NTC'73 Record, Atlanta, 
Georgia, November 26-28, 1973, 
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Figure 3, Method No, 3 for Command Authentication 
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the detector performs the same encoding process on the received codeword 
and checks the resulting parity with the received parity. 

The first method is easy to jam while the other two methods are 
harder to jam. The third method is probably the hardest to jam, but is the 
most complex. The second method is recommended since it is relatively 
Ixard to jam and is of considerably less complexity than method 3. 

In analyzing the command authentication scheme in Figure 2 against 
jamming, it is felt that the multiplexing scheme used might provide a sus- 
ceptibility to jamming. The multiplexing scheme is shown in Figure 4. The 
synchronization chaimel is sent as one 32-bit segment to each format frame. 
The voice channels labeled VI and VE in Figure 4 are 128-bit segments sent 
five times per format frame. For an encoded command word length of 
128 bits. Figure 4 indicates that a complete command word is transmitted 
in four 32-bit segments during a format frame. The basic unencoded com- 
mand word is 48 bits in length. Hence, with 50 format frames per second, 
all the data rates are achieved with the command data rate (48 bits /frame) 
(50 frames /second) equal to 2,4 kbps. It is possible, then, to add 80 bits 
of redundancy to each 48 -bit command word to give the 128 bits allowed for 
the command word per format frame. 
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A block diagram of the time -division multiplexed system is shown 
in Figure 5. The Viterbi decoder is switched into the system for the weak 
link through the TDRS relay satellite, in order to maintain acceptable bit 
error probabilities on the voice channel. The COMSEC unit may be 
switched into the system for Air Force missions. The strong direct link 
will exhibit such a low bit error probability that, even with the error 
extension of the COMSEC unit, the bit error probability will still be very 
small. As an example of the error extension phenomenon, consider a 
hypothetical COMSEC unit which outputs a burst of length 100 containing 
a 50 percent error rate for each input error. The possibility of jamming 
occurs by forcing an error before the command information in the multi- 
plexed format so that only the first few information bits are changed due 
to the error extension. Similarly, the last few information bits in the 
first 3Z bits of the multiplexed format could be changed without the error 
extension affecting any other bits. The possibility of jamming occurring 
is whether the parity bits corresponding to the changed information bits 
can be changed, causing an incorrect command to be accepted. This possi- 
bility would occur if the parity bits corresponding to the information bits at 
the beginning or end of a 3 2 -bit segment are also at the beginning or end of 
a 32-bit segment so they can be changed in spite of the error extension. To 
eliminate this jamming possibility, the information bits could be placed at 
the center of each 32-bit sequence, that is, 12 bits in the center of each 
32 -bit sequence. In this way, no information could be changed without 
changing many parity bits due to the error extension and thus rejecting the 
command. It should be feasible to arrange the information bits and parity 
bits so that the possibility of jamming is extremely small. 




Figure 5. Multiplexed Synchronization, Voice, and Command Channel Block Diagram. 
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Summary . -The probability of error for a convolutional coded coxn- 
munication system can be predicted using throe approaches- The first 
approach is using Gallager’s [l] exponential bound for random codes. 

The second technique is the use of transfer function union bounds 
developed by Viterbi [2]. Finally, an approximation to the transfer 
function developed by Huth and Weber [3] can be used. This paper com- 
pares the three approaches for predicting performance and presents 
results for Viterbi decoding bit error probability, burst error proba- 
bility, and decoder memory length. 

Introduction . -To design a convolutional coded conxmunication sys- 
tem, techniques are required to predict probability of bit error (I^j) 
versus energy per bit to single-sided noise spectral dcxisity (Eb/No) for 
the modulation/demodulation scheme employed. In addition, it would be 
desirable to predict the nature of the errors out of the decoder. That 
is, do the errors occur in bursts and, if so, -what is the probability of 
an error burst of a given length. 

There are three basic analytical techniqvies available for predicting 
Pb versus Eb/NQ. The first technique is using Gallager's exponential 
bound for random codes [l]. This bound allows quick comparisons of 
various modulation /demodulation and encodir.g/decoding schemes in 
terms of the basic coding parameters of code rule (R) and encoding con- 
straint length (K). Therefore, it can be used as an initial optimization 
technique. 

The second technique is to analyze the distance Ixotween codewords 
using a generating function (transfer function) developed by '’'iterbi [2]. 
This analysis technique has been used to obtain tight upper bounds on the 
error probability for the binary symmetric channel (BSC), the additive 
white Gaussian noise (AWGN) channel, and general memoryless channels. 
These bounds can be used to optimize the communication system param- 
eters for a particular convolutional code and Viterbi decoding. While the 
bounds on error probability using the transfer function approach are very 
precise, the bounds are difficult to calculate for long constraint codes. 
Also, fo.* even moderate constraint lengths, the transfer function bounds 
are evaluated in matrix form and the individual codeword distances are 
not available to determine the effect of decoder memory truncation or to 
determine the statistics of the errors at the ouLjxut of the decoder. 

The third technique uses an approximation ol the transfer function 
developed by Huth and Weber [3], which is easy to calculate even for long 
constraint length codes. The resulting bounds are as tight as Viterbi 
transfer function bounds and the individual codeword distant t's appear in 
the bound so that it is possible to determine decoder meitiory ti'uncation 
effects and output error statistics. 
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Random Coditig Pi*ediction Techniques . -Random >:oding bounds on 
error probability were developed by Gallagcr as easily calculable func- 
tions that did not require detailed knowledge of the code distance proper- 
ties, Gallager's function assumed that the decoder matched its decoding 
rule to the channel statistics. However, when the actual channel statis- 
tics are unknown at a given time, it is desirable: (a) to choose the mod- 
ulation format and the corresponding demodulation technique; and (b) to 
optimize the demcoulation parameters and the decoder metrics to mini- 
mize the probability of error for all possible channel statistics. Con- 
sidering the decoding metrics, let d(}r,x) be a decoding metric function 
where d(y, x) s 0 for all received sequences ^ and all possible transmitted 
sequences x, IfjCm the transmitted codeword for message ni and ^ 
is the received sequence, then the decoder correctly chooses if 

d(y, X ) > d(y, X ,) for all m' n\ . (1) 

m 


Let Pj.[error[m, be the probability of error conditioned on the mes- 

sage m entering the encoder, on the selection of a particular sequence 
2£m mth codeword, and on the received sequence y; . This proba- 

bility can be upper bounded by using a random coding argument similar 
to that of Gallager by letting the codeword be chosen randonily with 
probability Q(Xrn^' The resultant bound Is 


P < 
e, m 


(M-l) 2 
1 






-m 




m- J 


Q(x) ^/dCxT^ 


X 


( 2 ) 


where M is the number of possible codewords and P( 3 j;[x) is the channel 
transition probability. For a memoryless channel, note that 

N N 

= n = n 

n=l n=l 

where yjj takes on the values 0, 1, . . . , J-1 and x^^ takes --n the values 
0, 1, . . . , I-l . Let Q(k) be an arbitra i*y probability assiKuiiiont on the 
channel input alphabet and let each K Lter of each codeword of length N 
be chosen independently so that 


Q(x) 


N 


n 




(-t) 


Upon substituting (3) and (4) into (2), we get 


where 



P 

e, m 


^ ^-N(Rj3-R) 


-log^ 



n-1 


lu=o 


Q(k) 


P(ilk) 
N/d(j,k) . 


‘I-l 

D Q(i) s/d(J,i) 
-i=0 


(5) 

( 6 ) 


and 2^R = M. Note that R is the code information rate and Rp is ilie 
largest rate sucli that reliable communication is possible. Furthermore, 
Rd = ^comp “ ^0^^^ (Gallager's function) when d(j, i^) - P(j|k) which 
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corresponds to a maximum likelihood decoder but requires knowledge of 
the channel statistics* 

Viterbi [2] has derived the bit error probability bound for Viterbi 
decoding using Gallager*s function. For unknown channel statisLics, 

Viterbi* s bound can be rewritten using Rjj as follows: 

where K is the convolutional constraint length and R is the code rat. . 

Using (6) and (7), an initial optimization of the communication system 
can be performed since the modulation and channel characteristics deter- 
mine P(j|k) and the decoder metrics such as soft decision, hard decision, 
or list quantization used in the Viterbi decoder determine d(j, k). 

Transfer Function Union Bounds , -The transfer function of a convo- 
lutional code, introduced- by Viterbi [2] in 1971, is essentially a Hamming 
weight enumerator of certain types of codewords. Specifically, it enum- 
erates all the closed paths in the code's state diagram [3], The general 
form of this function is 


CO 

T(D,N) = a D^f, (N) (8) 

K^Clf K IC 

where is a polynomial of the form 

f, (N) = c, N + c, . + . . . c, .N^ + . . . (9) 

k kl k2 kj 

and d£ is the fi^ec distance of the code [2,3]. The general term in (8) has 
the following interpretation. For the given code, there are code- 

words that diverge from the all -zeros path at some point and r emerge 
with it later, and these codewords have Hamming weight k and are 
caused by input seque^nces with j "1*‘ bits. 

The trausler function has its usefulness in upper bounding the bit 
error probability of the code* To derive this bound, let be the proba- 
bility of incorrectly choosing a weight k codeword over the all-zeros 
codeword at the point of their remergence. If the weight k codeword 
was caused by j input "1** bits, then the bit error probability for this 
event is If there are type, the probability of the 

union of these bit error events is upper bounded by Finally, 

the probability of the uruon of all possible bit err.vr events*; upper 

bounded by the sum of the probabilities of these events, 


P^ < 
B 


CO c) 1, ( N ) 

k=d£ 


N=1 




From (9)> note that the derivative ten i in ( I 0) merely causes the remain- 
ing terms to be multiplied by the exponents of the variable N, and these 
exponents are just the number of bit errors. 

Suppose there exiaied a bound on each Pj^ of the form 

Pk “ ’ 


k e dj. . 


( 11 ) 


where x is the signal-to-noise ratio and a is a function of x and the 
modulation type* Substituting (11) into (10) and reari*anging terms gives 


CO 


exp(ad^x) Y, aj^exp(-akx) 
f k-clf 


yyN) 


0N 


lN=l 


Comparing (12) and the derivative of (8) with respect to N yields 

, , V 9T(D,N) 

3N ' ■ 


N=l, D=exp(-ax) 


( 12 ) 


(13) 


Thus, for the bound in (13), T(D,N) or its derivative need only be evalu- 
ated at some fixed values of D and N, a possible task. 

As an example of evaluating (8), consider the simple rate R=l/2, 
constraint length K=3 encoder in Figure 1. The contents of the last two 
register stages on the right constitute the encoder state, hence there are 
four states defined by ^r=-00, b = 0l, £=10, and d = ll. The state transition 
diagram for this code is shown in Figure 2 where the solid lines mark a 
transition caused by an input "0", the dotted lines mark a transition 
caused by an input "1", and the two binary symbols on each line repre- 
sent the code output for that transition. Looking at this diagram, the 
first few terms in the transfer function series can be written down. The 
path ab££ has weight 5 and was caused by one input "1”. The two paths 
a be b c a and a]jd£a have weight 6 and were caused by input sequences 
containing two input ”1’* bits, so the next term in the series is 2D^N^. 

This counting procedure could be continued, but it is easy to see that 
it is not too efficient, especially for more complex codes. A better way 
is to obtain directly a closed form expression for T(D, N). First, redraw 
Figure 2 so that state a is split into both source and sink nodes. Next, 
label each transition with the dummy variuMe D raised to the weight of 
the code output for that transition and N raised to the number of input ”1" 
bits that caused the transitioii. This construction yields the diagram 
shown in Figure 3. Now T(D,N) is just the transfer function of this dia- 
gram regarded as a signal flow graph, Using standard signal flow graph 
reduction techniques results in 


T(D, N) = 


D^N 
1 - 2DN 


= D^N + 2D^N + . . . + 2dV“^ + . . . (1-1) 


which checks with the two terms gotten by. counting. 

Although the above* method is straightforward, it is not suitable for 
computer computation. For the present example, define the state vari- 
ables x^, xjj, Xj,,x^ corresponding to the states a, £, and d, respectively. 
Then, from Figure 3, the following node equations hold: 


X = D^x ; X. = ND^x +Nx ; x = Dx, +lJx ; x = NDx, l-NDx . (15) 

a c b acebdd b d 

rp 

Next, define the state vector x = (xj^, x^,, x^)-^ , Then (15) can be rewritten 
in compact form as 


X “ Ax + B and = Cx , 


(16) 


where 
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A 


0 

N 

0 

D 

0 

D 

ND 

0 

ND 


is called the state transition matrix, and where 


B = (ND^ 0 0)^ and C - (0 0) . 

Solving for x in (16) gives 

X = (1-A)"^B = (I + A+A^ + . . . )B 


(17) 


(18) 

(19) 


where I is the identity matrix. Equation (19) lends itself very easily to 
computer evaluation. Obviously, the sum in (19) is truncated to a finite 
number of terms such that the truncation error is negligible compared to 
the sum. Now, once (19) and (16) have been evaluated, T(D, N) is equal 
to Xg^. ^ 

These bounds have been shown to be extremely tight for the AWGN 
channel. In this case, the bit error probability is bounded by 


P 


B 


< 



9T(D, N) 
8N 


N=l,D=exp(“E /N ) 
s 0 


( 20 ) 


where Eg/No is the binary coded symbol energy per single-sided noise 
spectral density, and 


(l>(x) 


■f 


exp(-v^/2) 


n/2tt 


dy 


( 21 ) 


For bit error probability less than lO"'^, the bound is within 0. 1 dB of 
the actual measured results for an AWGN channel. 


Transfer Function Approximations . -The use of the slate transition 
matrix to evaluate the transfer function union bounds loses much of the 
information contained in the individual terms of the transfer function 
series notation. An approximation developed by Huth and Weber [3] 
alleviates the problem and allows additional useful results to be obtained. 
It has been shown [4] that the minimum weight Dp of all codewords 
of length i which have not returned to the all -zeros state (i.e. , urrmerged) 
is bounded by 

[d^; +h^] Djf S id^i’ + b^J for all f > 1 , (22) 

where do is the minimum average weight per branch, and are 
biases^ \k] is the jimallest integer greater than or equal to x, and [xj is 
the largest integcj^ less than or equal to x. 

For some length thi“ path of smallest weight is lower bounded by 
do^ + bi* The path with the maximum weight is estimated by where 

the code rate is R==l/n* All intermediate weight paths of length j^,are 
approximated by appropriate combinations of these two extreme weights. 
The generating function for the weights of these unmerged codewords is 


(23) 


do+bj 


0 . 




dob 


A^(D,N) = ND (ND +D ) (ND +D ) 


wAere a 


= M 


and b = 


M 

L 2 J- 


The probability of ei*ror due to truncation, of the Viterbi decodv'ir 
memory at length I has the effect of causing the Viterbi decoder to 
choose among unmergea paths to determine the correct path. Thus, 
the bit error probability of error due to truncation Is 


9A^(D,N) 

9N 


N=1 


i; (j) E(.) (i+i-j+b)p 

i = 0 j:=0 


(24) 


where k = bj + fdo(l + i+j)l - 1 - i - j) . 

Using similar techniques, the probability of error due to the merged 
paths can be estimated:r ’■ 


B(D,N) = £ a D''t (N) + E (d‘'° +Nd“)‘'^^ 

'' k=J 


k=df 


(25) 


where the free distance dr is the minimum distance of the merged [>aths, 

T - '^ ll 

and I is defined as | , 

Thus, the approximation to the merged transfer function consists of 
the first few terms of the transfer function in (8) and an approximation to 
the higher distance codewords similar to the unmerged codewords. 
Therefore, the bit error probability due to the merged paths is: 


9B(D,N) 


9N 


afj^(N) 


= V d P 

N=1 Ic k 9N , . , 

« k/2 , ,, k/2 


k=J i=0 ' j = 0 ^ " 


where r = b^ + + i + j)1 +n(k - i - j) . 

The total bit error probability is appf oximated by 

9A^(D,N) ^ 




9N 


9N JN=1,D*'=Pj^ * 


(26) 


(27) 


Using (27), the bit error probability can be computed for arbitrary Viterbi 
decoder path memory as shown In Figure 4, The convolutional code used 
is the rate 1/3, K=7 co-le found by Odenwalder [5] to have df=14. The 
measured distance values are dQ = ll/17 and bj =2. Figure 4 presents a 
family of curves for various information bit energy per single -sided 
noise spectral density (Ei^/Nq). Note that the results agree very closely 
with the results of simulation. 
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Figure 5 presents the bit error probability versus Ej^/Ng for two 
decoder memory lengths and compares the predicted results with mea- 
sured results of a Viterbi decoder operating at the output of a STDN 
receiver [6]. Also shown are the results for the Viterbi decoder oper- 
ating on digitally generated Gaussian noise. The Viterbi decoder quan- 
tizes the received signal to eight levels while the Huth -Weber approxima- 
tion is for infinite quantization. Allowing for the 0.2-0.Z5 dB loss 
usually associated with 8-level quantization, the Huth-Weber approxima- 
tion is very close to the measured results for the ideal Gaussian noi.Te at 
large Ejj/Ng. At small E^/Ng, the union bound technique used to derive 
the transfer function bound results in a somewhat pessimistic prediction. 
However, for the region of Pj 3 < 10"^, the agreement is very good. Note 
that measured STDN receiver results are within 0. 2 dB of the ideal 
Gaussian noise when the bit synchronizer degradation is accounted for. 
Thus, the Huth-Weber approximation is an excellent prediction technique 
for real channels. 

The undetected errors of the Viterbi decoder tend to occur in buruts. 
This phenomenon has been observed from the operation of Viterbi 
decoders but also can be predicted from (25). Suppose an incorrect path 
is chosen by the Viterbi decoder that first differs from the correct path 
at the i^^ bit and does not have K-1 identical consecutive bits until the 
i + L+K-1 bit. Then, the decoder will make burst error of length L. But 
the sum of undetected errors in (25) is just the sum of the various lengths 
of paths that remerge with the correct path after L + K-1 bits. There- 
fore, the probability of a burst error greater than L is the sum of those 
paths that remerge with the transmitted path in lengths greater than 
Li + K-1, Let P(b>L) be the probability of a burst error greater than oi 
equal to L. Then, from (25), for L :S J -K+1, 


I ni 

P(b>L) < E E 

i=L+K-l k=df 


+ E E(;)p. 


where g|(k) is the number of paths of length i and weight k and 
s = bj + fdg(l+j)l +n(i-j). For L,>J -K+1, 

CO i 

P(b>L) < E E • (29) 

i=L+K-l j = 0 

It is possible to simplify the approximation in (29) for large L by bound- 
ing Pg as in (20) for the AWGN channel. Therefore, for large L, P(b£L) 
is of the form 


P(b>L) 


F E 

i=L+K-l 


< fg^'*'^“Vi-g , 


where F = <p(jz{d-+b ) E /N ) and G = exp ( -d E /N ) + exp(-nE /NJ. 

UsO sO 


Conclusions . -Three analytical techniques have been developed to 
predict the probability of error for a convolutional coded communication 
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system. Quick comparisons of modulation and coding parameters can be 
performed with the Gallager random coding exponential bounds. Finer 
parameter optimization can be obtained by using Viterbi transfer function 
vinion bounds and the Huth-Weber approximations to these bounds. The 
Huth-Weber approximations provide a techniqw to predict decoder mem- 
ory truncation degradation and burst error probability, as well as bit 
error probability. These analytical approximations are important to 
optimize the modulation/coding parameters at very low probability of 
error where measurements are impractical. It was noted that the 
approximations are extremely accurate over tlu; region 10"^>P|^>10"^ 
where measurements have been made. Since analytical techniques are 
used to predict performance, communication system optimization can be 
obtained with little effort. 
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SUMMARY . - Error correcting coding is required for 50 Mbps data link 
from the Shuttle Orbiter through the Tracking and Data Relay Satellite 
System (TDRSS) to the ground because of severe power limitations. Convolu- 
tional coding has been chO'^°n because the decoding algorithms (sequential 
and Viterbi) provide signivicant coding gains at the required bit error 
probability of 10"6 and can be implemented at 50 Mbps with moderate hard- 
ware. While a 50 Mbps sequential decoder has been built, the highest data 
rate achieved for a Viterbi decoder is 10 Mbps. Thus, five multiplexed 
10 Mbps Viterbi decoders must be used to provide. a 50 Mbps data rate. This 
paper discusses the tradeoffs which were considered when selecting the 
multiplexed Viterbi decoder approach for this application. 

INTRODUCTION . - A wideband communications link from the Space Shuttle 
Orbiter to ground via the NASA TDRSS is currently planned. This link will 
be used for On-Orbit transmission of up to 50 Mbps of scientific data and, 
on a time-shared basis, will also accommodate television, analog scientific 
data, experiment or maintenance tape recorder dumps, etc. The operating 
frequency will be in the Ku-band (approximately 15 GHz). 

In order to accommodate digital data rates on the order of 50 Mbps, 
under the constraints which exist due to limitations on Orbiter antenna 
diameter (present planning is for a 20-inch diameter deployable parabolic 
dish) and on Orbiter power amplifier output (approximately 50 watts), the 
use of coded, coherent PSK has been selected. The large performance gains 
associated with convolutional coding, the simple spaceborne encoder design, 
and the relative simplicity and availability of ground station decoders has 
led to the choice of convolutional codes over block codes. Due to a 
restriction in available TDRSS repeater bandwidth to 225 MHz, a code rate 
of 1/2 has been selected. 

This paper addresses the problem of selecting a convolutional encoding 
and decoding approach that provides the best possible performance without 
severe penalties in cost, complexity, reliability, etc, 

COMPARISON OF SEQUENTIAL AND VITERBI DECODING PERFORMANCE . - Convolu- 
tionaV coding with either sequential or Viterbi decoding has been shown to 
provide significant performance gains for communications over the memoryless 
white Gaussian noise channel. Both sequential and Viterbi decoders have 
been constructed for multi-megabit data rates and are currently available 
commercially. Soft-decision Viterbi decoders for short constraint length 
(K = 7) convolutional codes are available at information rates up to 10 Mbps, 
and provide coding gains in the range 4-6 dB. Figure 1 shows the measured 
performance of the Linkabit LV7015 Viterbi decoder, which operates at a 
maximum information rate of 100 Kbps. This measured performance coincides 
with the theoretical performance of the LV7017 Viterbi decoder, which 
operates at up to 10 Mbps. The measured performance data was obtained in , 


the Electronic Systems Test Laboratory at the NASA Lyndon B. Johnson 
Space Center, and reflects the performance achievable for a practical PSK 
receiver (non-ideal carrier reference) with the degradation due to symbol 
synchronizer performance corrected out. 


Figure 1 also shows the theoretical performance of the Linkabit LS4157 
sequential decoder, which is a commercially available hard-decision device 
originally developed for a 40 Mbps application, but marginally capable of 
operation at up to 50 Mbps. 

As indicated by the performance curve, the coding gain provided by 
sequential decoding increases with decreasing information rate. This is 
because the perfonnance of the sequential decoder is primarily determined 
by the buffer memory size and the speed advantage of the decoding logic 
(the ratio of computation rate to data rate), and the speed advantage is 
higher at lower information rates. 

Figure 1 indicates that the soft-decision Viterbi decoder significantly 
outperforms the hard-decis-ion sequential decoder at high error probabilities 
(greater than about 10-5) and is still somewhat superior at an error prob- 
ability of 10 - 6 . The almost 2 dB improvement in sequential decoding per- 
formance offered by soft-decisions would make sequential decoding a much 
more favorable candidate for the application being considered, but use of 
soft-decisions would greatly reduce the maximum information rate at which 
the decoder could operate. Soft-decision sequential decoders are not 
commercially available, but the technique is straightforward and has been 
demonstrated at low information rates. 

POSSIBLE 50 MBPS DESIGN APPROACHES . - Based on the performance consider- 
ation^ described above and on the availability of decoder hardware, the 
following options for the 50 Mbps Shuttle Ku-band link were evaluated: 

a. High-speed (> 10 Mbps) Viterbi decoder(s) 

b. Multiplexed 10 Mbps Viterbi decoders 

c. 50 Mbps hard-decision sequential decoder 

d. Multiplexed low-speed (< 50 Mbps) soft- decision sequential 
decoders 

The results of these evaluations are summarized in subsequent sections 
of this paper. 

HIGH-SPEED VITERBI DECODER(S) FOR 50 MBPS OPERATION . - Currently, a 
Viterbi' decoder capable of 50 Mbps~operation can not be implemented with 
off-the-shelf components. The arithmetic unit of the Viterbi decoder is 
limited to about 35 Mbps using ECL off-the-shelf components and to about 20" 
Mbps using SOS components. Thus, with off-the-shelf components, at least 
two Viterbi decoders would have to be operated in parallel on multiplexed 
data if a 50 Mbps capability is to be provided. If fewer than five parallel 
Viterbi decoders are used, then a potentially expensive decoder development 
effort would be required to provide a decoder for operation at information 
rates greater than 10 Mbps. 

The arithmetic unit could achieve 50 Mbps operation if it is imple- 
mented with a custom LSI chip using ECL or SOS. For a K = 7 decoder, 64 
identical custom arithmetic chips would be required to implement the 
arithmetic unit. These custom LSI chips would significantly decrease the 
complexity of the arithmetic unit as well as allowing 50 Mbps operation. 


However, this custom LSI chip development would be relatively expensive; 
also finding a semiconductor manufacturer to develop the custom LSI chip 
could be extremely difficult without a guaranteed large scale production 
follow-on. 

50 MBPS OPERATION USING MULTIPLEXED 10 MBPS VITERBI DECODERS . - One 
obvi ous approach which would provide a 50 Mbps information rate capabi 1 i ty 
without additional decoder development is to utilize five parallel 10 Mbps 
encoders and five parallel 10 Mbps Viterbi decoders. This implementation 
is illustrated functionally in Figure 2. The performance provided by a 
parallel approach such as this would be the same as would be provided by a 
single 10 Mbps Viterbi decoder. 

It should be noted that the encoder portion of Figt re 2 can be imple- 
mented with a single 31-stage shift register encoder operating at 50 Mbps 
with no mux/demux units. 

HARD-DECISION SEQUENTIAL DECODING AT 50 MBPS . - Sequential decoding of 
convolutional codes ls arpowerful decoding technique. This technique has 
an advantage over Viterbi decoding in that long constraint length convo- 
lutional codes can be used since the code constraint length has only an 
insignificant effect on the decoder hardware complexity. TTierefore, 
undetected errors can be made highly improbable for even small values of 
EbTNoT A? noted previously, the performance of sequential decoding is 
primarily determined by two parameters: the buffer memory size and the 

speed advantage of the decoding logic (the ratio of computation rate to 
data rate). 

In the Linkabit LS4157 hard-decision 50 Mbps sequential decoder, the 
received information bits are delayed by a 65,000 bit MOS shift register. 

The delayed information bits are exclusive-ORed with the output of the 
sequential decoder, resulting in the corrected information sequence. This 
large buffer is necessary because the number of computations to decode an 
information bit is a random variable. 

The decoding logic is implemented using the fastest commercially avail- 
able logic family, MECL III. The operating speed of the decoder logic is 
70 megacomputations per second. Nearly all output errors are the result of 
the sequential decoder not having a correction to the information bit at 
the end of the buffer soon enough. When this event occurs, the information 
bit is output with no error correction. For the sequential decoder to 
"catch up," at least the next 500 bits will be output from the buffer with 
the channel errors uncorrected. Occasionally, 5000 bits will be output 
from the buffer without the channel errors corrected. Thus, the errors 
output from the ouffer will tend to occur in long bursts. 

50 MBPS OPERATION USING MULTIPLEXED SOFT-DECISION SEQUENTIAL DECODERS. - 


From strictly a performance standpoint, this approach is superior to 
the other three approaches described previously. There are several dis- 
advantages associated with soft-decision sequential decoding, however, and 
these will be briefly outlined. 

The basic arithmetic in a soft-decision sequential .decoder is much 
more complicated than in a hard-decision sequential decoder, because of the 
increased branch metrics. In addition, a larger buffer memory is required 
to store the same number of branches as a hard-decision device. The net 
result is that a 3-bit (8 quantization levels) soft-decision sequential 
decoder would be three times as complex as a hard-decision sequential 


decoder, and would operate at a much lower information rate. A good 
estimate is that the information rate for a 3-bit soft-decision sequential 
decoder (using commercially available logic) would be limited to about 8 
Mbps, such that six or seven of these decoders could be paralleled to 
provide a 50 Mbps capability. The performance provided by such an approach 
would be about 2 dB better than that provided by a single 50 Mbps hard- 
decision sequential decoder. 

A 2-bit (4 quantization levels) soft-decision sequential decoder would 
only be twice as complex as a hard-decision device and could be made to 
operate at a somewhat higher information rate (perhaps 12.5 Mbps), such 
that four decoders could be paralleled to provide a 50 Mbps capability. 
However, a 2-bit soft-decision sequential decoder provides only a 1.25 dB 
improvement over hard-decisions, so the performance advantage would not be 
as significant. 

Soft-decision sequential decoding suffers from two additional dis- 
advantages. First, the performance of a soft-decision sequential decoder 
is highly susceptible tOfVariations in receiver A6C voltage, so somewhat 
greater care must be exercised in the design of the modem with which the 
decoder must operate. Secondly, the memory introduced by the costas loop 
of a PSK receiver results in a pronounced degradation in the performance of 
a soft-decision sequential decoder; therefore, the use of interleaving 
(2000 bits or more) may be required in order to realize the full potential 
performance gain. 

DESCRIPTION OF SELECTED APPROACH . - Based on the considerations noted 
above, the following conclusions were made regarding the four options for 
50 Mbps operation which were evaluated: 

a. High-Speed (> 10 Mbps) Viterbi Decoder(s ) - This approach was 
considered unacceptable because of (1 )' the developmental effort required to 
provide high speed (< 50 Mbps) Viterbi decoders using off-the-shelf com- 
ponents or to provide the custom LSI chip required for a 50 Mbps Viterbi 
decoder, and because of (2) the fact that the performance would be no 
better than that achievable using commercially available Viterbi decoders 
operating at 10 Mbps. 

b. Multiplexed 10 Mbps Viterbi Decoders - This approach was 
considered acceptable and was selected for implementation, 

c. 50 Mbps Hard-Decision Sequential Decoder - This approach was 
considered unacceptable because of the poor performance characteristics of 
hard-decision sequential decoders. 

d* Multiplexed Low-Speed (< 50 Mbps) Soft-Decision Sequential 
Decoders - This approach was considered unacceptable because of (1) the 
development costs associated with providing soft-decision sequential decoder 
hardware, and (2) the relatively high susceptibility of soft-decision 
sequential decoders to non-ideal conditions such as receiver AGC variations 
and receiver memory. 

The 50 Mbps design approach selected for the Shuttle Ku-band link is 
summarized in somewhat greater detail in Figures 3 and 4. Figure 3 
illustrates the single shift- register implementation for the encoder and 
also indicates the bit-by-bit interleaving scheme which was selected. 

Using the encoder interleaving scheme shown in Figure 3, the parallel 
decoder configuration is totally self-synchronizing. The demultiplexer 
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which precedes the parallel decoders is a single deconmutator which merely 
routes every fifth incoming symbol to each decoder. Each decoder establishes 
its own branch synchronization and the resulting decoded sequences are com- 
bined with an output multiplexer. Figure 4 illustrates the details of these 
operations. 

One aspect of the selected parallel Viterbi decoder approach which is 
of concern is that of reliability. The reliability of a configuration 
consisting of five decoders operating in parallel is obviously much lower 
than the reliability of a single component decoder. This decreased 
reliability may necessitate some sort of fault detection and sparing 
provisions* possibly even an automatic replacement scheme. 

CONCLUSIONS . - It is felt that the selected 50 Mbps encoding/decoding 
approach offers a reasonable compromise between performance and system 
complexity. The coding gain provided by this approach is relatively high 
(5.3 dB at a probability of error of 10”°), the system degrades gracefuHy 
(a coding gain of 2.7 dB is available at a probability of error of 10 “^), 
the spaceborne encoder is extremely simple to implement, and the ground- 
based decoders are available commercially and can be combined in a straight- 
forward manner. Such an approach is not limited to 50 Mbps data rates. It 
appears reasonable to extend the parallel 10 Mbps Viterbi decoder concept 
to systems requiring operation at data rates v;ell in excess of 100 Mbps, 
although reliability becomes of increasing concern as the number of parallel 
decoders becomes large. 
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Figure 1 

Performance of the Linkabit LS4157 Sequential 
Decoder and Linkabit LV7015 and LV7017 Viterbi 
Decoders. 
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Figure 2. Functional Configuration for 50 Mbps Operation 
Using Multiplexed 10 Mbps Viterbi Decoders 
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SIGNAL DETECTION AND ANGULAR SEARCH PROCEDURE 
FOR SHUTTLE KU-BAND COMMUNICATION STUDY 

Sergei Udalov 


1 . 0 INTRODUCTION 

To establish the Ku-band communication link between the Shuttle 
and the TDRS, the antennas of these two satellites must be pointed at each 
other. Various techniques for pointing the TDRS antenna at the Shuttle 
have been considered and the latest informatioii available to Axiomatix 
indicates that initial pointing can be accomplished without requiring the 
Ku-band Shuttle antenna axis to be in perfect alignment with the axis of 
the TDRS antenna. During this initial pointing, however, the effective 
EIRP available from the TDRS to Shuttle is only 40 dBw instead of the 
nominal 5 0 dBw. But, before the TDRS makes its final pointing correction, 
the Shuttle antenna must search out the random designation errors which 
exist between its computer -estimated linc-of-sight (LOS) to TDRS and the 
actual pointing direction of the Shuttle antenna. In this appendix, we conside 
the procedures for searching out the residual designation errors of the 
Shuttle antenna and we estimate the time requirements for these procedures 

The main argument presented here is that the times required to 
detect the presence of the TDRS signal within an arbitrary beamwidth (1 dB 
beamwidth is assumed) of the Shuttle antenna are short compared to the 
typical time requirements for mechanically scanning the LOS uncertainty 
volume. 

Before proceeding with the quantitative analysis, we would like to 
point out that the overall probability of detecting and acquiring the TDRS 
signal consists of two terms: 

Pj^ = Probability of detecting the presence of the TDRS signal 
when the Shuttle antenna is pointed, within a specified 
offset, at the TDRS, and 

P.y, = Probability that the true LOS (to TDRS) is contained within 
the angular volume searched out by the Shuttle antenna. 

The Pj^ is determined by forward link power budget and by the detection 
procedure used at the Shuttle receiver for determining the presence of the 
signal when the Shuttle antenna is within a predetermined offset from true 
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LOS to TDRS. The is determined by such factors as Shuttle antenna 
beamwidth and beam shape (directivity characteristic), the LOS designation 
uncertainty, and the angular search pattern used. To a first-order approxi- 
mation the overall detection probability is 





( 1 ) 


and with proper selection of detection and search techniques, one can insure 
that P,p S P^(min) where P^(min) is an arbitrarily selected minimum prob- 
ability of signal acquisition. 

Acquisition inrthis appendix is defined as termination of the angular 
search procedure and the initiation of the angular autotrack. 

2. 0 SIGNAL DETECTION 

Before considering the angular scan procedure for the Shuttle 
antenna, we will analyze the procedure for detecting the presence of the 
signal within the antenna beamwidth. For this, we will assume that the 
angular scan over the designation uncertainty volume has been initiated 
and that the purpose of the signal detection unit is to stop the scan when the 
received signal exceeds a preset value. 

The simplest and the fastest technique for determining the presence 
of the signal is to use an energy sensor, such as a square law envelope 
detector followed by an integrathig filter and a threshold crossing indicator. 
Because this technique does not require either the code sync or the carrier 
lock, it can be applied directly to the incoming spread spectrum signal. 

The simplified block diagram for implementing such a detection scheme 
is shown in Figure 1, 

The first step in determining the performance of this detector is 

to determine the power budget for the forward (TDRS to Shuttle) link. 

Table 1 shows the pertinent parameters of tliis budget. From this table, 

we see that the available P /N_ is 64, 8 dB-Hz. 

rec 0 







Table 1, Forward (TDRS-to-Orbiter) Ku-Band Link 

Parameters 


Parameter 


Units 


1. TDRS EIRP dBw 

2. Space Loss (13.775 GHz; 22^ 786 nmi) dB 

3. Orbiter Antenna Gain dB 


4. 

5. 


6 . 


Total Receive Power, (Sum 1 through 3) 

Orbiter System N*bise Temp (10 log 

-23 

Boltzman Constant (10 log 1.38x 10 ) 


dBw 

dBrK 

dBw/K-Hz 


7. Orbiter Noise Spectral Density, 
(Sum 5 and 6) 


dBw/Hz 


8, Total Received Power /Noise Spectral 

Density, P /N„ (7 minus 8) 
rec 0 


dB-Hz 


Applies to TDRS acquisition phase only. 
Based on a 20 -inch diameter parabola. 


Value 
40* 
-207.7 
34, 6 
-133. 1 

30.7 
-228. 6 

-197.9 

64.8 


If we assume a rectangular IF filter*' at the receiver, the loss as 
compared to a matched filter is about 0. 9 dB and the Bt product is !• 37. 

The IF SNR for rhe spread spectrum signal of 14. 1 Mbps (71. 5 dB) is 
therefore 

(SNR)jj, 64.8 dB^Hz - 0.9 dB - 71. 5 dB = ^7.6 dB (2) 

and the corresponding IF bandwidth is 

= 1.37 X 14.1Mbps = 19.3 = 20 MHz. (3) 

It! T 

The -7,6 dB value for the IF SNR corresponds to the condition when the 

Shuttle antenna is pointed directly at the TDRS. 

To insure that signal detection occurs when the beam of the Shuttle 

antenna is off the true LOS within a predetermined amount, we have to set 

the parameters of our detection circuit accordingly. For example, if we 

want to achieve a specified probability of signal detection within the -1 dB 

antentta pattern contour, we set our detection circuitry to declare a "signal 

present" condition when the IF SNR is -8. 6 dB. Similarly, if we decided to 

achieve signal detection within the -3 dB antenna beamwidth contour, we 

have to adjust the parameters of the signal detection circuit so that the 

threshold crossing indicator is actuated at -10,6 dB IF SNR. 

For any given combination of probability of detection, Pj^, and 

the probability of false alarm, the parameter which we can vary is 

the nuniber of samples integrated, N, This number is approximately equal 

to the ratio of IF bandwidth to the bandwidth of the post -detection integrating 

filter, i. e., N = B/b. Figure 2 shows N versus IF SNR for two values of 
-9 

Pj^ and P|^ = 10 . The corresponding integration time, T^, based on the 

relationship N s B/b, is also shown in this figure.’*'* From this figure, we 

Butterworth or Chebichev design. 

Computational details are given in Addendum A, 


Figure 2. Number of Samples Integrated and Approximate 
Integration Time vs. IF Signal-to-Noise Ratio 


N (Samples Integrated) 
i r •• 


10 , 000 

9, 000 


P, = 10 
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see that to achieve specified detection criteria (use = 0.999) with the 
-1 dB antenna pattern (SNR^^ == -8,6 dB), about 4,100 signal samples must 
be integrated and the corresponding integration time is about 220 micro- 
seconds. As shown in the next section, the integration time is short com- 
pared to dwell times achievable with a mechanically scanned antenna. Also 
the false alarm time for the parameters used is 
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fa 



220 X 1 0 ^ sec 


10 
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220 X lO seconds 
or about 6l hours . 


This time is well withfin- any reasonable specification. 

We may also consider a special case for v/hich we do not have a 

paramp in the sum channel of the Shuttle Ku-band receiver. This may be 

a case where only a single 2l6 kbps data stream spread by the 14. 1 Mbps 

PN sequence is required to be transmitted to the Shuttle. The noise figure 

of the sum channel will then bo about 8 dB and the corresponding noise 

temperature is T = (6. 3)(290*^ K) = 1825° K or 32.6 dB-K. This will 
s 

degrade the (SNR) by about 1,9 dB, resulting in (SNR) of -10,5 dB. 

Ir Ir 

From Figure 2 we see that, at this SNR, about 10, 300 samples have to be 
integrated for = 0. 999 and the integration time is increased to about 
51 5 microseconds. This is still a short time compared to the dwell times 
encountered with a mechanically scanned antenna. 


3. 0 ANGULAR SEARCH PROCEDURE 


Having established the criteria for detecting the presence of the 
signal within the 1 dB antenna beamwidth, we now consider the angular 
search procedure which insures that the designation uncertainty volume 
is scanned to provide a predetermined probability of true LOS intercept. 
For this we consider the relationship between the antenna beamwidth and 
the designation uncertainties. The baseline numbers are 

0_ S Z. 8° s Antenna 3 dB beamwidth (20-inch parabol.. 
® dish) 
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and = cr^ = 0. 933° - Designation uncertainties of Az 

and El axes (based on 3 ff = 

= 2. 8°) . A E 

It is reasonable to assume that a. and are independent, have zero 
mean, and are Gaussian-distributed along their respective coordinates. 

The conversion to polar coordinates then results in a probability distribu- 
tion of the true LOS around the designated LOS. This distribution has the 
form: 

G exp (- 0 /2cr ) 

p(6,o) = , 0 S e < ® (5) 

f •• 0 < o < 2 tt 

where o* = cr = 0'„j and 0 represents the antenna boresight deviation 

J\ llr 

along any redial direction away from the designated LOS. The 6-dependent 
portion of the two dimensional distribution given by (5), i, e, , 

p(0) = exp (- e^/ZOg^) , 0 < G < m (b) 

is known as the Rayleigh distribution and its general form is shown in the 
upper part of Figure 3, The angular direction, a, of this distribution, 
however, can be anywhere in the interval 0 £ a ^ 2 tt. In other words, one 
can picture the probability distribution of the true LOS as a Rayleigh- 
shaped "wind vane" which is "pivoted" around the designated LOS and 
whose direction a is uniformly distributed across the interval from 0 to 

2tt, 

Because we do not know the direction along which the true LOS 
distribution lies, we have to search out all possible directions and also 
move the Shuttle antenna beam outwards until the cumulative probability 
of intercepting the true LOS reaches the predetermined value, such as 0.99 
or 0.999, etc. The cumulative probability of intercepting the true LOS as 
the beam completes successive revolutions and moves "outward" in the 
6-direction is 
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Figure 3, Spiral Scan Provides Means for Interception and 

Detection of True LOS 
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P„(0T ) = f “ (- B^/cTq^) de = 1 - 

*'o 0 

where 0^^ is the limit of the "outward" beam displacement. 

To accomplish the required search, one can use cither a stepped 
circular or a spiral scan. With a stepped circular scan the beam is moved 
outward along either the Az or the El axis by a preselected angular distance 
and then a circular scan is initiated. If the true LOS is intercepted during 
this scan the antenna motion is either terminated immediately after signal 
detection, so that the autotrack can lock on to the true LOS, or the detection 
coordinates are storeef in a memory and the beam returned to these coor- 
dinates upon completion of the scan. If detection is not indicated, the beam 
is moved outward and a circular scan along a wider radius is performed. 
This procedure is repeated until the true LOS to TDRS is intercepted and 
detected. 

With a spiral scan, which is more common,*’' the antenna beam is 
moved outward from its initial position in a spiral pattern which eventually 
covers all possible locations of the true LOS up to a predetermined value 
of the cumulative intercept probability. Figure 3 shows how an idealized 
antenna beam is moved in a spiral pattern to "cut through" the probability 
density function of the true LOS. The lower portion of Figure 3 shows how 
each successive "cut" through the density function contributes to the cumu- • 
lative probability of true LOS intercept. 

Two facts must bo noted while we are considering the idealized 
search depicted in Figure 3. One is that the "height" of the idealized 
antenna beam shape (shown as a cylinder in the figure) represents the 
i. e,, the probability of signal detection as determined by the parameters 
of the detector shown in Figure 1 and the available IF SNR, IF < 1, 
then the contribution of the individual "cuts" through the true LOS distri- 
bution must be weighed accordingly. The second fact evident from Figure 3 

Stepped scan requires fast changes in antenna direction and thus 
increases servo loading. 


is that even without scanning there is a finite probability that the true LOS 
is within the beamwidth regardless of a. The value of this probability 
depends on the relationship between the antenna beamwidth and the "spread” 
of the Rayleigh distribution. 

First let us consMer the Shuttle's antenna pattern. For this we 
will employ a commonly used assumption that, within a considerable portion 
of the main lobe, the antenna pattern can be approximated by a Gaussian 
form: 

G(G) = exp (- 2.776 0^/6g) (8) 

where G(6) = normalized boresight gain (maximum gain) 

6_ - antenna 3 dB beamwidth 

B 

6 = angular deviation from boresight . 

The loss, in dB, as a function of 0 is therefore 

L(dB) = 10 log [exp (-2.776 e^/0g)] . (9) 

The plot of L versus 6 for 0^ = 2.8*^ is shown in the upper left corner of 
Figure 4. The cumulative probability of true LOS intercept, Py« versus 0 
is also shown. From plots of Figure 4 it is evident that using -1 dB antenna 
pattern contour the initial value of prior to scan is only 0. 30. Note that 

only a half of the antenna pattern is shown prior to scan. As the spiral scan 
commences the boresight begins to move outward and the -1 dB contour 
begins to "cut” into the p(G), thus contributing to the cumulative probability, 
Py(0)» of true LOS intercept. 

Suppose we decide to make Py - 0.9999. From the Py(9) plot of 
Figure 4 we see that this requires "spiraling out” the -1 dB point from 0.8 
to 4.0 degrees. If we spiral out at the rate of 0.8° per 360° turn, the 
number of turns required, n^, is 

Distance to maximum 6 _ 4.0° - 0.8° _ 3.2° 

" AO per turn 0.8°/turn 0,8°/turn 




V 
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With these parameters, we barely touch the 4° contour and thus 
we are not certain what dwell time we have to stop the scan when we reach 
our selected limit. Thus, we increase the rate of A6 advance. If we let 
A6 be 0.9® per turn, then n^_ becomes 


3, 


= 3.56 turns . 


”t 0. 9®/turn 

But we let the scan continue to an even 4 turns. Then the excess is 


01 ) 


An. = 4-3.56 - 0,44 turns . 
t 


The time spent beyoncf the 4® limit then becomes 


At = An X 


1 ^ . J_ 

turn/ sec o x ^ 

s 


( 12 ) 


(13) 


As explained in the next section, the turn/ sec is the scan frequency f 
which determines the dwell time of any antenna beam contour (-1 dB, 

-3 dB, etc. ) and also the time required to spiral out the boresight to any 
required cone limit. For the example chosen, the -1 dB contour and the 
boresight are spiraled out to 4. 2® and 3. 6®, respectively. Also as stated 
in the next section, if we assume servo loop bandwidth of 1 Hz and we want 
the servo to pull on the first pass over the true LOS, the scan frequency f 

s 

should be 0.425 Hz. The time of overscan, t , is then 

os 


os 


= An X T — 7 
t 0.4 


0,44 

0.425 


= 1 . 04 seconds . 


(14) 


This time is far in excess of the detection time of about 220 microseconds! 

The time to complete the four turns at the scan frequency of 0.425 Hz is 
about 9.4 seconds. 

The combined minimum probability of true LOS intercept and 
detection for the procedure described is 

P,j.{min) = = (0.999)(0. 9999) = 0.9989. (15) 


i . 


J 
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The actual probability may be somewhat larger because of the overscan 
we have allowed for the P^term, Nevertheless, the procedures outlined 
above permit one to select various combinations of Pj^ and P^ to obtain 
the desired value of P^(min). 

When the energy detector commands the Shuttle antenna scan 
mechanism to terminate the search, the antenna is stopped and the track 
mode is initiated. There may be a pull-in transient, but since we have 
assumed that the antenna is stopped within the -1 dB antenna pattern con- 
tour, the servo operation should be in its linear region. A conservative 
estimate is that the pull-in time should bo approximately one third of the 
inverse of the servo loop bandwidth. For a 1 Hz loop, this implies that 
about 0. 3 second will bo required for the servo loop to go from scan to 
track mode. 

At the time when the search is terminated,, the value of P /N„ 

rec 0 

may be 1 dB below the nominal 64. 8 dB-?lz value, i. c. , it may be about 

64 dB-Hz. The recent analysis'*' of the Shuttle monopulse system indicates 

that, with the paramp in the sum channel only and no despreading, the 

maximum RMS error due to thermal noise at about 64 dB-Hz is less than 

0.05 degree. Consequently, we can safely assume that the thermal noise 

is not a limiting factor for establishing a reliable monopulse autotrack at 

P /N- values as low as 64 dB-Hz. 
rec 0 

When the Shuttle antenna begins to track the TDRS signals, com- 
mands to correct the TDRS pointing are sent to the ground via an S-band 

link. As the TDRS pointing is zeroed in at the Shuttle, the P /N„ avail - 

^ rec 0 

able to the Shuttle receiver increases by about 10 dB. The code sync and 
doppler acquisition then take place and the Ku-band link communication is 
established. 


C, L, Weber, "Shuttle Monopulse System for Ku-Band Communi- 
cation Signal from TDRS - Part II, " Axiomatix Report No, R7509-4 (under 
Contract NAS 9-14614), September 26, 1975. 


N. 


^ I 


I 


15 


t 






4. 0 EXTENDED SCAN COVERAGE AND COMPATIBILITY 
WITH RADAR SEARCH MODE SCAN 


In the preceding section, we have considered only the requirements 
for searching over the residual designation uncertainties of the TDRS LOS 
acquisition. In practice, however, various factors may increase these 
uncertainties and a scan over a wider angular area may be required. But 
the search requirements for some of the radar modes of the integrated 
radar /communication system do require scan over relatively wide angles 
and thus the adaptation of the radar scan modes to communication angular 
search appears feasible. In particular, the ±5® (10° cone) radar search 
mode for the cooperative target should be considered for an adaption to 
the communication angular search mode. As shown below, scan over 
wider areas can also be performed for the communication mode at the 
expense of increased search time. For all of the quantitative discussions 
which follow, it is assumed that a spiral scan is utilized. The detailed 
analysis of spiral scan is given in Appendix I. 

There are two basic spiral scans to consider. One involves con- 
stant angular velocity, i.e., the rate at which the boresight is rotated 
around its starting position is constant. This results in a progressively 
increasing tangential velocity as the boresight is spiraled out. As a con- 
sequence, the dwell time decreases and eventually it may fall below the 
value required by the servo to acquire the signal source (TDRS or a radar ‘ 
target) on the first pass. The antenna must then make at least one or 
more revolutions before acquisition. Figure 5 shows the relationship 
between the scanned angle (conical coverage) and the dwell time versus 
scan frequency (rev/sec) and the elapsed scan time. The equations used 
are 


e = A0 f t 

s 


and 


-i. ^ /J-)^ I 

Ztt A6 t 


(16) 


(17) 







where 


0 


- scanned an^lc 

f = scan frequency 

s 

0 - 3 dB antenna beamwidth 

B 

A6 = radial advance in deg/rev . 

For data in Figure 5, wo used 6„ =2. 8° and AO = 1 . 6°/sec which corre- 

JD 

sponds to -1 dB overlap. 

As an exannplc for use of Figure 5, assume that we want to deter- 
mine the time to scan out a ^5° cone, with a 1 Hz second-order servo loop, 

A second -order servo requires an acquisition time of about one third of its 
inverse bandwidth. Therefore, a 1 Hz loop has to be scanned at an 0.3 Hz 

f • _ 

rate to insure that , *• the 5 limit the dwell time is not less than about 0. 3 
seconds. The scan time required to roach this limit is then about 1 0. 5 
seconds. Slower scan frequencies require longer times to reach any given 
scan limit but they also provide longer dwell times,' 

Because the dwell time at the maximum scanned angle is a major 
factor which affects servo acquisition i>ehavior at the scan limit, the rela- 
tionship between the parameters wh. 7 V. determine this dwell time is important. 
As derived in Appendix B, the dwell tinie at ma cimum scan limit is 


t 


dm 



2 



(18) 


The time to scan out to the maximum angular limit is 


t 

m 



A0 f 

s 


(19) 


where A6 is the radial advance in deg/rev. 

Figure 6 shows plots of t and t , versus f for various values of 

m dm s 

0 and A0. The baseline value of 6o = 2,8° is assumed for all plots of t , 
m ri ' dm 

From Figure 6 is is -evident that the 0. 9°/rev overlap chosen as an example 
in Figure 4 is conservative and that a shorter scan time can be obtained, 
i. e. , 5,2 seconds versus 9. 4 seconds if the AG is increased to 1 . 6°/rcv, 


Time - Seconds 
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The 1. value corresponds to -1 dB worst case beam overlap and 

thus does not violate the energy detection criterion. This applies to the 
1 Hz loop and 0.3 second minimum dwell time assumed as the requirement 
for servo acquisition. 

So far we have discussed the scan which utilizes constant rotational 

velocity. To achieve constant tangential velocity and constant dwell time 

over the scan path, the scan frequency has to decrease inversely with time. 

This, however, implies that at the start of the angular search the scan 

frequency must be infinite. To avoid this requirement, the actual scan 

program starts out with a constant angular velocity scan and then switches 

r ' 

over to a constant tangential velocity scan. The switchover takes place 

when the minimum allowable dwell time is reached. For example, from 

Figure 5 we see that, for f - 1 Hz, the 0.3 second dwell time is reached 

s 

after about 1 second of scan time. If the servo bandwidth remains constant 

at i Hz, compatible with the initial f‘ = 1 Hz, the scan frequency must then 

s 

decrease to keep the dwell time at a constant value of approximately 0.3 
second. The program for the time -variable scan frequency is then 

0 £ t < tjj (20) 

to^t (21) 

where is the initial constant scan frequency and is the time after which 
the scan frequency is slowed down.. 

Figure 7 shows the application of the time -variable scan frequency. 
The criterion for the data shown, is to achieve scan limit with dwell 
time no less than 0,3 seconds and a I Hz servo loop. Three cases are 
shown. Case 1 starts with £ = 1 Ha constant scan frequency and then 
switches to constant dwell time scan after = 1 second. Case 2 starts 
with f = 0. 5 Hz and switches over at ” 4 seconds. The Case 3 scan 

S U 

uses a constant f value of 0,31 Hz. The most significant result indicated 
s 




^^-4 M - - 1 . . 4 1 I . 
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by this example is that, for a selected dwell time at a given search limit,* 
the time to reach this limit is independent of the scan method used. But 
because the constant tangential velocity scan is most efficient from the 
standpoint of providing uniform dwell time over the major portion of the 
search area, the switchover to this scan should take place as soon as the 
minimum allowable dwell time is reached. Thus, for a 1 Hz loop and the 
f' = 1 Hz, the switchover should take place at t. = 1 second. 

The data presented so far indicates that times required to search 
the uncertainty angles for TDI^ acquisition are short compared to radar 
search times. This is due to the fact that SNR available from TDRS for 
signal detection during acquisition is much higher than the SNR of radar 
targets. Consequently, the scan time limitation for TDRS acquisition is 
determined by servo considerations rather than signal detection. The 
approach for radar search is the same, however, as that for TDRS acqui- 
sition and only the dwell time requirements are different. We therefore 
recommend that the radar search pattern be adapted for TDRS acquisition 
when the integrated system is in the communication mode, 

5.0 CONCLUSIONS 

The energy detection method for determining the presence of the 
spread spectrum TDRS signal within the beamwidth of the scanning Shuttle 
antenna requires integration times in the order of several orders of mag- 
nitude shorter than the dwell times associated with a mechanically scanned 
antenna. The latter times are typically in the order of seconds or fractions 
of seconds. Consequently, the time required to search out the initial desig- 
nation inaccuracies of the Shuttle antenna is determined by the servo drive 
capabilities and scan pattern rather than the signal detection technique. 

This applies even to the reduced power mode when the effective EIRP from 
TDRS is 40 dBw instead of the nominal 50 dBw. 


In this case, approximately 0. 3 second a*- a 5® limit. 


For a Gaussian -distributed, uncorrelatcd, and zero mean value 
deviation and azimuth designation errors, the transformation to polar 
coordinates presented in this appendix provides a convenient analytical 
technique for estimating the cumulative probability of true LOS intercept 
as the function of spiral (or stepped circular) scan strategy. 

Because the dwell times required for the angular acquisition of 
the radar targets are typically in excess of the times required to detect 
TDRS LOS, we recommend that these scan programs be used, with minor 
modifications, for the establishment of the Ku»band communication link. 


■Ni 
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ADDENDUM A 

SIGNAL INTEGRATION REQUIREMENTS VERSUS 
SIGNAL.-TO-NOISE RATIO 


The energy detector used to determine the presence of the spread 
spectrum TDRS signal within the beamwidth of the scanning Shuttle antenna 
employs a square law device followed by an integrating low pass filter and 
a decision threshold device. The square law device is preceded by an IP 
filter whose bandwidth is approximately matched to the bandwidth of the 

spread spectrum signal. Because the incoming signal has a constant envelope 

r 

there is no need to perform its actual sampling at specific intervals to deter- 
mine the number of pulses integrated required to achieve a given probability 
of detection versus IF signal -to-noise ratio. In this case, it is the ratio of 
IF to low pass filter bandwidth which determines the equivalent number of 
pulses, or more specifically, signal samples integrated. The analysis used 
for discrete sampling can be applied to this problem, however, to determine 
the effective number of samples integrated and the corresponding integration 
time required to achieve a given probability of detection. For this, we use 
the following equation;*^ 


SNR 




IF 




(A-1) 


where 


SNR,„ = the rms S ignal -to -noise ratio in IF 
IF 


N = number of pulses (or samples) integrated 


,-l 


(|> ” function related to false alarm probability 

= function related to probability of detection P^ , 

The functions <)j ^(P, ) and ^(P,) are given by Figures 10,7-1 (p, 368) and 

la a 

10,7-2 (p, 369), respectively, in the reference cited. 


DiFranco and Rubin, Radar Detection , Prentice -Hall, Inc,, p, 3^7, 


For one of the cases, we select P , = 0.999 and P* = lO" , For 
“1 “ 1 ^ 

these values, ^ ~ ^^fa^ ~ 5.95, Substituting in Equation 

A-1; 


SNR 


5.95 - (>3,10)] ^ 9.05 

^fw 


Solving for N, we obtain: 


(SNR)^^ 


(SNR) 


2 * 


For example, for SNRjj, = -10 dB, the numeric value is 0, 1. Substituting 
this value into A -3 we get: 

N = = 8190 samples integrated . 

( 0 . 1 )^ 

Since the sample time is approximately equal to the reciprocal of the IF 
bandwidth, the corresponding integration time will be: 


N _ 8190 

^IF 20x10^ MHz 


s 410x10 sec 


or 410 microseconds . 


Similar procedures are used for other values of (SNR) _ and also for the 

li? 

case of Pj = 0.99. Figure 2 of this report shows the plot of N and corre- 
sponding for (SNR)jp ranging from -7. 5 dB to -10. 5 dB. 
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.SPIRAL SCAN ANALYSIS 
by 

Gaylord K. Huth 
Sergei Udatov 

This appendix analyzes the relationship between the dwell time, the 
total scan time, and the scan rate for the following types of spiral scans: 

(1) spiral scan with constant rotational velocity; 

(2) spiral scan with constant tangential velocity; 

and (3) hybrid scan which starts out with constant rotational velocity 

and then switches over to constant tangential velocity. 

It will be shown that the hybrid spiral scan offers the best compromise 
between the antenna servo power requirements and the dwell time. Several 
numerical examples pertinent to the Shuttle Ku-band radar /communication 
system are provided. 

l.O INTRODUCTORY ANALYSIS - STEPPED CIRCULAR SCAN 

Before considering the spiral scan, it is instructive to develop first 
the relationship between the major parameters for the stepped circular scan 
which is a discrete approximation of the continuous spiral scan. Figure I 
defines the basic parameters of the stepped circular scan. After each com> 
plete rotation about boresight, the antenna ir stepped to the next radius for 
the next rotation. The number of complete rotations N can be found by 
noting that 

= -r +N(?„(1-A) (1) 

m z is 

where $ is maximum scan limit, is the antenna beamwidth, and A is 
m ±s 

the scan overlap. Equation (1) can be rewritten to determine N; 


= beamwidth (3 dB) 

Q = maximimv scan limit 
m 

A = beam overlap 

tj>. s angle of antenna dwell 
bn i^ rotation 


Figure 1. Stepped Circular Scan Increases the Dimensions of 
Bo resight Loci by Discrete Increments 




4.N 




^3(1 -A) 


where fx] is the largest integer greater than x. 

The dwell time on the ith rotation is denoted by t , and is given by 


d. £. 2-7 f, 2n 
11 i 


4l^g(l-A) £.»r 


4i(l - A) 


where £. is the rotational £requency (i. e . , revolutions per second) £or the 
ith revolution and the radius o£ the ith revolution is id— (1 - A)* Note that» 
for reasonable values of A (i. e . , 30% or less)» equation (3) can be rewritten 


*dj ~ fj2n'i(l-A) 

with only a 2. 2% error or less for i-1 and 0. 5% error or less for i=2. 

* 

The total scan time T , i, e, , the time required to reach the scan limit 

, ■ s ■ 

boundary is 

^1 

= ‘d + £ f +“’■ <51 

“ “0 i=l ‘i 

where t^ is the dwell time on the beamwidth at bore sight and V is the time 
0 

to step between scan radii. For constant rotational velocity f • 

c 


T = t 
8 d 




^m“ 2 




(t'f 


For constant tangential velocity, 


^drt ^d. f. 2jri(l-A) 
0 11 


^i t, 2ffi(l-A) 





( 9 ) 


Therefore* T t 
8 d 


Fl +2f^(l-A) £ i] 

oL i=i J 


+ Nt 


= t [l + JT(1 -A) N(N + d] +Nt. 


But* using equation (2) for N* 


6‘ 


N(N + 1) = 




( 10 ) 


Hence, for constant tangential velocity* 


T = 

B 


1 + 




0m- 


0^(1 -A) 


2A)j‘ 


9 




;b 

2 


^g(l -A) 


( 11 ) 


2. 0 C03MSTANT ROTATIONAL VELOCITY SPIRAL SCAN 

With the antenna axis and the monopulse comparator arranged for 
tracking along two orthogonal axes* such as Az and El* the spiral scan 
can be Implemented by driving each iixis with an orthogonal sinusoid of 
progressively increasing amplitude. The boresight displacement along 
each of the orthogonal axes is then 

AAz = 0{t) cos Oit 

AEI = 0(t) sin Oit (12) 

where d(t) is the time -dependent function which determines radial bore- 
sight displacement and ct>- Znf = angular rate of boresight displacement. 

5 

The will be referred to henceforth as the ’’scan frequency." This fre- 
quency may be equal to or may be smaller than the servo bandwidth* but 
it should not exceed the latter. 

In polar coordinates* the boresight locus is defined by a spiral 
0{t)<0(t) shown in Figure 2. Li Figure 2, 0„(1 - A) is the boresight dis- 

. . J3 

placement per one complete turn of 9 b(t), 


0„(1 - A) = deg/rev 

Jd 

6g = beamwidth (3 dB) 

0 = maximum scan limit 

ni 

A = beam overlap 


Figure 2. Spiral Scan Increases the Dimension of the Boresight Locus 

in a Continuous Manner 






Let tf(t) be proportional to ^(t). Then, 9{t) is given by 


^ (1-A) 0„{1-A) 

^(t) = ot = 0g(l-A)f^t. (13) 


From (13), we can determine the time T required to spiral out the bore- 

S 

sight to any given maximum angular (conical) limit, 9 . Thus, 


T 

s 


9 


e - 

m 


B. 

2 


f - A) 

s B 


(14) 


Equation (14) can be used directly only if f (i. e., scan frequency) is con- 

s 

stant. Note that (14) is very close to (6) for the stepped circular scan with 

T equal to 0, Thus, (14) applies to a constant angular velocity scan. For 

a constant tangential velocity, f is time -dependent and (14) has to be modi- 

s 

fied accordingly, as explained in Section 3. 0. 

As the boresight position spirals out toward a given maximum limit, 
the dwell time decreases if the (i> is held constant. The dwell time can be 
computed by considering the tangential^ velocity along the spiral at the 
corresponding Vctlue of 6(t)t Thus, 

= 0(t)a>. (15) 


The dwell time is then considered as the time required by a predetermined 
antenna beamwidth contour, 9^ (-1 dB or -3 dB, etc.) to traverse a direc- 
tion in space as the boresight moves across that direction with velocity V. , 


^B ^B 1 ^B / 1 1 

*d ^ 0g(l -A)fg<ut "" 17 0g(l - A) t * 

Of particular importance is the dwell time available at the position 
where the boresight reaches the maximum scan limit 0 ^ . The expression 


The radial component due to spiraling out is not included here 
because it is generally low compared to the tangential component of the 
total velocity along the spiral. 


I 


* 


for the value of this dwell time is obtained by substituting (14) into (16): 


J 


0. 


yi - A) 


dm 


2»7 0^(1 -A) £2 


d. 


1 

2vf 


e 


B 


e 


m 


0--T 

txi 2 


(17) 


Equations (13) and (16) can be used to plot the angle scanned and the cor- 
responding dwell time as a function of time and scan frequency, f , Cor- 

6 

respondingly, equations (14) and (17) can be used to plot the maximum time 

to reach the scan limit, 0 , and the corresponding t , as the function 

m dm 

of f . 


3. 0 CONSTANT TANGENTIAL VELOCITY SPIRAL SCAN 




For a spiral scan with constant tangential velocity V , from (15)- 

V 


with o - 2irf , 
6 


V^ = 


(18) 


where 0(t) is the instantaneous radius given by (13). Hence, 


V. = 2ir0 (1-A)f''t . 
t i5 s 


(19) 


e 


■ The total scan time T can be found from (5) for the circular scan 

s 

by letting t equal zero and the sum be an integral over the continuously 
changing 6^(t). Therefore, note that maximum value of 0(t) is 


0B 

max 0(t) = 0^ - . 

itl 2 


( 20 ) 


The instantaneous number of turns completed, n, is the total radius the 
spiral has reached divided by the advance per turn. That is, 


9 


n = 


m 


ag(l-A) 

and the incremental increase is 


( 21 ) 


I s 


i: 


f . 


I C 


t;; 




8 


dn - 


dg(t) 


6>g(l-A) * 


(22) 


Thusj the sum in (5) becomes the integral 


"X 


^m" J_ dfl(t) 




(23) 


Neglecting the radial component due to spiraling out, then from (16), 


t 

d 


6 


B 


Znf ^(t) 
s 


or 


e 


f = 


B 


s ■ Znt^eii) ’ 


Substituting (25) into (23) results in 


T = 
s /1 2 


<?3(1-A) 


Jo 


6{t)de(t) , 


Also, note that by substituting (25) into (13), one obtains 


^„(1-A)t 

e\i) = ® 


2zrt . 


Therefore, 


^(t)dd(t) = 


0 g(i -A)dt 


4wt 


and l^y rearranging (27), 


t = 


Znt^e (t) 


^3(1 - A) 


The r.iaximum value of t is 


(24) 


(25) 


(26) 


(27) 


( 28 ) 




9 


max t = 




Substituting (28) and (30) into (26) gives 


-k£ 


max t 


nttd -ej2) 

which is very close to the result for stepped circular sccUi with t equal 
to zero. 


4. 0 HYBRID SPIRAL SCAN 

The design of a scan procedure typically requires the maximization 

of the dwell time for a given total scan time T . Figure 3 compares the 

s 

dwell time for a spiral scan with constant rotational velocity [i. e. * scan c(t)] 
and constant tangential velocity [i. e * , scan a(t)]* It can be seen that con- 
stant tangential velocity maximizes the dwell time but at i9(t) - it 
requires infinite rotational velocity as predicted by (25). Therefore, to 
overcome the infinite rotational velocity at 0(t) and to increase the dwell 
time from that obtained by constant rotational velocity, a hybrid scan tech- 
nique is proposed. The hybrid scan is shown in Figure 3 [i. e. , scan b(t)] 
where initially the scan is performed at a constant rotational velocity and 
then is s^tched over to constant tangential velocity. To compute the total 
scan time of the hybrid scan, let the constant rotational velocity scan switch 
over to the constant tangential velocity scan at Therefore, the radius at 
switchover is 


- T 


(32) 




Dwell Time 


Switchover 


maximum scan limit 



frame time 

Scan a(t); f (0) infinite and then progressively decreases to maintain 
constant t^. 

Scan b(t): f constant until switchover, then decreased to maintain 
constant t^ for the remaining portion of the scan. 

Scan c{t): f^ constant to butt^ decreased continuously. 


Figure 3. Scanned Angle and Dwell Time for Constant Tangential Velocity 
a(t), Hybrid b(t), and Constant Rotational I^ate c(t) Spiral Scans 
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The total scan time is composed of the constant rotational velocity scan 
time T and the constcmt tangential velocity scan time T • From (14) 


sr 

and (26)* 


sv 


T s T +T 
s sr sv 


^r 2 


2nt 

J 0{t) d^(t) . (33) 


0g(i-A) 


Making the substitutions similar to equations (27) to (31), 




T = 


^r“ 2 




[(0^ - - {9^ - 0„/2)^3 . (34) 

2,. * » m Jo to 


For any given initial constant rotational velocity f and any final dwell time 

s 

the following boundary condition from (17) must hold 

e. 


1 


'B 


^d Zni 

® - — 
r 2 


(35) 


or 


e. 


e - 


B 




r ■ 2 ■ 2fff t, * 

8 d 


(36) 


Therefore* for given f and t * the half -cone angle reached at switchover 

8 d 


xs 




(37) 


3 d 


Substituting (36) into (34) 


T = 
a 


4Tr(l-A)f t, 
s d 


e|(i-A) 


( 38 ) 
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Note that the second term in (38) is the same as (31). Thus, for a given 

t^, the hybrid scan requires a total scan time greater than the constant 

tangential velocity scan by the addition of the first term in (38). If £ = <», 

s 

then (38) becomes the constant tangential velocity result in (31). Alter- 
nately, t^ can be found for a given by rewriting (38) as 


rrid^ - d^/Z) 

m £> 

9b(1-4) 


T t^ + 
s d 


= 0 


4;r(l -A)f' 


(■'\‘ 


and solving this quadratic equation for t^ . It is important to point out that 

solution of (39) for t^ depends not only on but also on such parameters 

as $ I On I a-nd f . Because T and 0 arc generally fixed for a given 
m li 6 s m 

mission, maximizing t^ for a particular value of 6^ can be performed by 

analyzing the relationship betweeii t^ and f^ . Such analysis reveals the 

minimum value of £ yielding an acceptable value of t , . In the next section, 

8 a 

several numerical examples of such an analysis are presented, with par- 
ticular emphasis on Shuttle Ku-band radar /communication system parameters 


5. 0 NUMERICAL EXAMPLES 

In the preceding sections, the analytical relationships for various 

spiral scan strategies were developed. It was shown that the hybrid scan 

provides an optimum compromise between excessive servo requirements 

and dwell time. We will now use the analytical results to obtain numerical 

data pertinent to hybrid spiral scan such as may be employed by Shuttle 

Ku-band radar /communication equipment during acquisition. Specifically, 

we will use equation (39) to obtain the values of dwell time t . as a function 

a 

of f for a set of various mission-dependent parameters T and 0 and 
s iT JT s m 

for antenna sizes ranging from 20 inches to 36 inches in diameter. Specific 
mis 8 ion -dependent parameters for Shuttle Ku-band equipment are: 
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» 


» 


t 


4 ^ 





rc 


(1) 

e = 

±30° 

m 


(2) 

= 

±20 

(3) 


±10° 


The first two are characteristic of radar acquisition for a passive target. 
The third (0^ = 10°) scan limit parameter is typical of cooperative target 

acquisition in the radar mode ^d TDBS acquisition in the conununication 
mode: 

Figures 4 through 7 show the solution of (39) for t, as a fxmction of 

a 

the scan frequency f and the miss ion -dependent parameters T and $ . 

s s m 

The most sighific 2 mt res\ilt evident from these figures is that, at the scan 

frequency of 2 Hz, the dwell times for all cases considered are near their 

maxima and exhibit only an insignificant increase at higher values of f^. 

As can be e:mected, the case of ^ = ±30° and T =30 seconds 

m 8 

sets the lower limit on f for all dish sizes considered. This is because, 

s 

in this case, the antennas have to scan the maximum value within the mini- 
mum allotted time. Hence, the f^ has to be relatively "fast" to spiral the 

beams across the scan volume determined by the 6 limits. 

m 

Figure 8 shows the minimum dwell time available from the hybrid 
scan versus antenna dish shize and the sets of mission -dependent param** 
eters used. Essentially, these "minimum" dwell times are taken at 
f =2 Hz, because above this frequency there is no significant gain in the 
dwell time. The 2 Hz frequency, however, is within the operating range 
of typical antenna servo drives. 

6.0 CONCLUSIONS 

Analytical expressions for various types of spiral scans have been 
derived and the hybrid scan appears to be most efficient from the standpoint 
of minimizing servo loop loading required to achieve "best" (i* e. , longest) 
dwell times for any given set of parameters. Specific numerical examples 


indicate that a Z Hz scan frequency should meet the dwell time optimization 
requirements for most missions performed by Shuttle Ku-band radar/com- 


mtmication equipment 
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FALSE LOCK PERFORMANCE OF COSTAS RECEIVERS 
by Marvin K, Simon 

i 

ft 

In. this appendiXf we first demonstrate that a Costas loop can exhibit | 

'jj 

a false lock condition in the sense that the loop can lock up on a data side- i 

band. In particular » for the case of random data« the loop can theoretically 

1 

lock at a frequency v/hich is any integer multiple of half the data symbol ' 

rate* An expression for the dc output of the lock detector is then derived 

i 

in terms of the Fourier transform of the symbol pulse shape, the squared 
magnitude of the lowpass arm filter transfer function and the loop phase 
error. This expression is then evaluated at the true lock and false lock 
conditions and numerical result's are given for the case of Manchester coded 
data and a single pole (RC) arm filter. 

Consider the Costas loop and associalted lock detector illustrated in 
Figure J-1. The input signal s[t, 6(t)] is a bi -phase modulated carrier of 
the form 

®[t,0(t)3 - f\IZS m(t) sin [u>pt + 0(t)] (1) 

where S is the average signal power, m(t) is a binary modulation (a i:l 1 

digital waveform), w is the radian carrier frequency, and e(t) - 6« +n«t | 

is the received carrier phase. The total received signal x(t) is then ] 

x(t) = s[t, e{t)] + n^(t) (2) 

where n.(t) is the additive channel noise which can be expressed in the 
form of a narrowband process about the actual frequency of the observed 
data, viz., 

n.(t) = <s^ |N^(t) cos [wpt + e(t)] - N^(t) sin [w^t +0(t)] J 



( 3 ) 












where N (t) and N (t) are approximately statistically independenti stationary, 
c s 

white Gaussian noise processes with single-sided noise spectral density 

N. w/Hz and single -sided bandwidth B__< w /2 . 

U HO 

For the true lock case, the in-phase and quadrature reference 

signals r (t) and r (t), respectively, would be given by 
s c 

r^(t) = 4T sin [w^t + ^(t)] 

r^(t) = nTz" cos [wgt + 0(t)] {4) 

where d(t) is the VCO's estimate of 0{t) andK^ is the reference signal 
power. In the false lock mode, one would have 

r^(t) = s/F sin [(w^ - «)^)t + ^(t)] 

r^(t) = a/F cos [(Wq - (o^t + S(t)] (5) 

t 

where is the false lock radian frequency as yet to be determined. 

Denoting the in-phase and quadrature phase detector (multiplier) 
gains by K^, then using (3) eind (5), the output of the quadrature 

phase detector is (ignoring second-harmonic terms) 

^^it) = K^x(t) a/Z Kj cos [(«q - o>^)t + 0(t)] 

= K^K^[4s m(t) - N^(t)] sin [w^t +?i(t)] 

while the in-phase arm phase detector output is 

€^(t) = K^ x(t) a/F Kj sin [(Wq - w^)t + ^(t)] 

= K.K [a/F m(t) - N (t)l cos [w,t + •^{t)l 

X IXl 8 1 

- K,K N (t) sin[wt +0(t)] 

1 m c •• f 'J 


( 7 ) 
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where 0(t) - 0(t) - ^(t) is the loop phase error. After lowpass filtering 
with the in-phase and quadrature phase arm filtersi these same signals 
become, respectively, 

H A(t;«,^) H h N (t;w,^) H 

SI 8p X 

z^{t) == KjK^In/S' G( p)[m(t) sin [oj^t + ^(t)]| - G(p)[N^(t) sin [o)^t + 

H 

+ |g(p) [Nc(t> cos [u)^t + </»(t)]]| 

K *n^(t; *\ H — H 

z^(t) = K^K^{^/S G(p)[m(t) cos [w^t + 0(t)]j - G(p)[N^{t) cos [w^t + 

1* H 

- KjKm{G{p)[N^(t)sin[w^t + <^(t)]i} (8) 

4 

Here, the signals A (t?(O.,0) and m (t; u>,, are the in-phase and quadra- 
ture signals emerging from passage of the data modulated false lock sinu- 
soids, m(t) sin (ta)^ + 0{t)] and m(t) cos [«^t + ^(t)], respectively, through 
the lowpass arm filters G(s), and likewise, N (t;w,, ^), N (t; u> , 0), 

A A ® ^ ^ CC X 

N (t; <t>f, 0), and N (t^ w , are the equivalently filtered versions of the 

SC X'.CS'X 

noise components. 

Squaring z (t) and z (t) and differencing these quantities produces 
c s 

a process whose statistical mean (over the data sequence) and time average 
gives the dc voltage at the lock detector output. From (8), the normalized 
signal component of this dc voltage is 


In what follows, we shall write differential equations in compact 
form by introducing the Heaviside operator p ^ d/dt. 

For convenience, we drop in dur notation herein, the dependence 
of on time. 




J 


. , 


€ 


dc 


= lim 


/ T /2 

® - m^(t;u,, ^)] dt 

TJZ ° ^ ® ‘ 


(9) 


where the overbar denotes statistical mean. It now remains to characterize 
in terms of the statistical characterization of the modulation m(t) anrf 
the filter transfer function characteristic G(ju>). 

The digital modulation m(t) xnay be represented as 


00 


r. 




i 


m(t) s a^p(t-nT) 


( 10 ) 


n=-eo 

where = ’^l is the data symbol sequence and p(t) is the symbol pulse 
shape, h.g., for Manchester coding, 


p(t) = 


"1 ; 


1 ; 


T 

fstso 


Tf 

Ost^f 


(11) 


Substituting (10) into the definitions of m (t; u,,0} andm (t:u., given in (8) 

C £ S . t 

results, after considerable manipulation, in 
m^(t;w^,<^) ^ G(p)[m(t) sin (w^t +0)) 

= 2Z ^ |G{jw)‘y[p(t -nT) sin (Wjt + 0)]| 


n=-co 


CO 


J(Wft+«^) 


= E 




where ^ denotes Fourier transform with ^ its inverse and 


( 12 ) 
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i: 


I- ■' 
I- 

f ‘ 

rr 

I' «. r. 


Similarly, 


/ CO 

0[j(“+ “,)]P(Me^"* I; 

ct> 

PjW = / al^w-c-jnPOWe^"* 1^ 

mc(t;u)£,0) - G(p)[m(t) C 08 (o>jt + 0)3 


(13) 


CO 


J(uj£t +0) 




n==-o) 


00 


Xx=^co 


+ E 


-j(«ft+0) 


(14) 


Squaring (12) and (14) and taking the difference of their statistical means 
gives 

—5 ^ — ](2«.ft + Z^) „ , , 

m^(t!u>j,^>) - m^(tjUj,^) = — £ £ ^^(t-nT) g (t-mT) 

tt=-0O m=-a> 

-j(2w£t+20) 00 {o 

j. e ^ aT^^p^(t-nT)p2(t-mT) 

(15) 


nss-oo m=-co 


Finally, taking the time average, we get 


= i i: E iTT {lim ^ f 

dc 2 ^ ^ n m\- T_ I 

n=-oo m=-oo ® *'" 

im ~ / 

-*w *0 ./-' 


1 /" ^0^^ 2j(w£t +0) 

r- I e 


T„/2 


^j^(t-nT)pj(t-mT) dt 


+ lim 
T 


Trt/2 -2j(u)£t+0) 


^^(t-nT) ^^(t-mT) dt} 


(16) 


7 


where 


and 


1 2j(w£t+0)^ 


lim ^ f 


Tq-*® ' 0 *'-Tq/2 


ii^(t-nT)gj(t-mT)dt 


gJ20 


lim 


^ r* /“G[j(» + «p] P( ja>) G[ j(-w'+ w^] 


sin [(«-«’+2wJ T./2 ] , , 

^ (w-««+2w )(T /2) ® 2ir 2ir 


(17) 


^ A 

-*® 0 -1 


►T„/2 -2j(cD£t+^) ^ 


0 *"-Tq/2 


P2(t-nT) dt 


= lim r /“ G[j(<w-<t)p] P(jco) G[j(-a>* -o>^)] 

•'-00 •'-CO 


v|(a,n-u»«m)T d® dw* 

(w-O)' -2<yp(TQ/2) ® 2ir 2ir 


In order to simplify (16) any further^ we must consider a specific statistical 
characterization of the data sequence |a^[ . For the case of random data* 

X* e* * 

* I ; m = n 

(19) 

0; miiin 


a a s: 5 
n m mn 


The double sum of (16) reduces to a single sum. Recognizing further that 


* -j(w-u>*)nT Zv ^ f . Zirk. 

y. - -;;r T (w-w«- -— ) 


( 20 ) 




C 


i t 


V. : . . 

t ' 

\ : 
i f 


S' : 


:; T*- 

■ >> 


V, " lim 


dc 


To-, 


‘ £ r G[j(u.'+,f+^)]p[j(«'+^n 

ks-«> J~C0 

eta[(2tr|-+ 2u,)(T /2J 

X a(j{.„'+,,)]p(.j,') j-5 L_o_' is. 


(2tt^ + 2w^)(Tp/2) 


+ e 


■j20 ^ 


E /“G[j(<*>' -«f + ^)] P[j(«* + ~)1 

-CO •^^co 


. ., •‘“t(2l.|-2,j)(V2)]o„.j 

X G[3(-w’ -wjJ P(-jw*) j- — 

(2tt| -2w^)(T^/2) 


( 21 ) 


ZtrK 


Now suppose that 2«^ = " - - where 0, 1, 2, . . , is a particular 
value of k . Then, the value of namely Vj^ t at this frequency is 

given by • 


dcj 

2v 




(22) 


or equivalently 



k^s 1.2,3,... (23) 


where Re { } denotes the real part of what is contained within the braces. 
Note that, for the case of true lock, i. e. , s 0, we get 


Vq = cos 


f;] ' ““ 2i> \GlMf ^ 


( 24 ) 





1 1 2 

vrhere S^(m) = ~ |P(jo))| is the power spectral density of the modulation 
m(t). For any of the false lock frequencies f^ =: w^/2ir = k^/2T, the corre- 
sponding lock detector output (signal component) is given by (23). 

A particular example will now be given to illustrate the theory just 
developed. Consider the case where m(t) is a Manchester coded data mod- 
ulation and G(s) is a single pole RC filter for which 




1g(Jw)1 


1 + (w/w ) 
c 


where <o = 2-trf is the radian 3 dB cutoff frequency of the filter and is 
o c 

related to its two-sided noise bandwidth by 


<«) = 2B. . 

c 1 


For the Manchester pulse described by (11)« its Fourier transform is 
given by 


P(jw) = jT 


2 wT 


After considerable algebraic manipulation and integral evaluations* we 
arrive at the following results 


V„ = c. [3 - 4. 


-2fffcT 
+ e 1 


Y, = - cos 
k 




-2»rfcT 


; k = l,3*5,... (29) 


= - cos Z(p 




-irfcT -2vfcT 
3 4 4e -I- e 


2rrfcT 


k-2,6,10,.,, (30) 


10 


\ = 


cos 2<i> 







-sf T -277fcT 


1 


3 - 4e + e ® 


‘ *{zt rf 


ZnicT 






ks 4,8,12*... (31) 


At this point, it should be mentioned that, i£ one was to consider 
the error signal into the loop filter F(p), i^ e. , 


ZQ(t) = z^(t) z^(t) 


(32) 


where z (t) and z (t) are given in (8), then the dc value of the normalized 
c 8 

signal component in ^Q(t)t namelyi 


^dc 


^ lim j- f 

T 0 J-T 


TJZ ^ 
0 A 


./2 


m^(t;w^,</>) m^(t;to^,<^) dt 


(33) 


would be given by (28) through (31) with cos 20 replaced by (sin 20)/2. 

In view of this, we may reach the following conclusion. Since, for any 
f^T, the braced quantity in Vq of (28) is positive, then as is well known, 
true lock corresponds to the lock points 0= 0 ,± )t, ±2ir, ±3w, ,,,, Further- 
more, since the remaining braced quantities in (29) through (31) are also 
always positive for any f T, then in the false lock mode, we must have 
cos 20 = -1 or 0 = * ff/Z, db3»r/2, ±5w/2, . . . . Thus, the lock points for any 
false lock condition are interleaved midway between those for the true lock 
state. 

Figure J-2 is a plot of V, (in dB); k = 0, 1,2,3,4, versus f T for 
the example given above. We note that, by continuing to widen the arm 
filter bandwidth beyond the optimum value which yields minimum squaring 
loss (f T = 1.4), the separation between and the false lock values con- 
tinues to increase* This property is desirable from the standpoint of dis- 
tinguishing false-lock from true-lock. However, widening the arm filter 
bandwidth also allows more noise to pass through, thus degrading the 
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signal -to -noise ratio at the input to the threshold detector. Thus, in con- 
clusion, widening the arm filter bandwidth is potentially a solution to the 
false lock detection problem up to the point where the signal -to -noise ratio 
which distinguishes between true lock and out-of-lock (as opposed to false 
lock) becomes too low for a specified level of lock detection performance* 



